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1.  This  publication  documents  the  computer  program  called  OOP  which 
is  a mathematical  tool  used  in  the  study  of  the  doppler  content  of  rever- 
beration. The  program  has  been  used  primarily  with  the  characteristics  of 
the  Torpedo  MK  46,  however,  the  characteristics  of  other  acoustic  systems 
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SECTION  I 

\ 

\ 

INTRODUCTION 

\ 1.1  Background 

Torpedoes  are  the  prime  conventional  antisubmarine  weapons  and 
they  use  acoustic  systems  to  detect  the  submarine.  Some  of  these  systems 
use  frequency  dependent  characteristics  of  the  returned  echo  from  the 
submarine  to  enhance  performance. 

1.2  Purpose 

A digital  computer  program  has  been  developed  at  the  Naval  Undersea 
Center  (NUC)  which  facilitates  the  analysis  of  systems  that  use  frequency 
dependent  characteristics.  The  program  is  called  DOP  and  is  used  in  conjunction 
with  two  other  NUC  developed  programs  called  SONAR  and  RAYSRT  which  are 
documented  in  references^l-afld-2.^  This  report  documents  the  DOP  program. 


1.3  Publications 

See  list  of  references,  page  32. 

SECTION  II 

...  — ~ GENERAL  DESCRIPTION 

2.1  Doppler  in  Reverberation 


When  a single- frequency  pulse  is  emitted  from  an  active  sonar  system 
on  a moving  platform,  the  reverberation  seen  by  the  system  is  spread  due  in 
large  part  to  doppler  effects.  Program  DOP  computes  the  spectrum  of  such 
reverberation  in  a refractive  medium.  The  doppler  content  is  computed  as  a 
function  of  the  speed  of  the  sonar  platform  and  (optionally)  of  circular  turning 
or  motion  of  the  scatterers,  or  both.  In  addition,  the  spectrum  of  the 
original  pulse  can  be  included  as  a spreading  effect,  since  even  a "single- 
frequency"  pulse  has  a harmonic  content  due  to  its  finite  duration.  The  energy 
level  is  computed  in  frequency  bands  of  specified  width  at  specified  times 
relative  to  the  transmitted  pulse.  Four  values  are  computed  for  each  band/ 
time  combination:  surface,  bottom,  volume,  and  total  reverberation.^ 

2.2  Boundaries 

For  boundary  reverberation  (surface  or  horizontal  bottom),  increments 
are  summed  in  random  phase  from  all  areas  of  the  boundary  returning  energy 
in  a given  band  at  a given  time.  Scattering  strength  is  a function  of  grazing 
angle,  and  all  combinations  of  paths  (direct,  refracted,  reflected)  to  the 
scattering  areas  are  included. 
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2.3  Volume 


Volume  reverberation  is,  at  present,  computed  in  an  unbounded 
uniform  medium,  summing  contributions  from  all  volumes  returning  energy  in  a 
given  band  at  a given  time. 

Both  boundary  and  volume  computations  consider  two-way  losses  In 
the  environment  and  average  transmit  and  receive  beam-pattern  losses  to 
each  incremental  scattering  unit.  In  addition,  optional  filtering  can  be 
applied  to  each  band  and  TVG  (Time  Varied  Gain)  action  can  bp  applied  at 
each  time.  Also,  total  energy  in  all  bands  at  each  time  is  computed  for 
surface,  bottom,  volume,  and  total  reverberation. 

SECTION  III 

INSTRUCTIONS  FOR  RUNNING  DOP 

3.1  General 

The  program  is  written  primarily  in  FORTRAN  IV  for  execution  on  a 
UNIVAC  1110.  The  program  is  not  self-contained  in  that  there  are  some 
functions  and  subroutines  that  must  be  supplied  by  the  user  to  match  the 
particular  system  being  studied.  (The  program  has  been  exercised  here  at 
NUC  with  the  characteristics  of  the  Torpedo  MK  46  Mod  1 and  an  experimental 
torpedo.)  A description  of  these  user  supplied  routines,  the  input  data,  and 
ancillary  programs  necessary  or  useful  to  the  execution  of  DOP  follows. 

3.2  User  Supplied  Functions  and  Subroutines 

For  any  specific  application,  one  or  more  of  the  following  FORTRAN 
functions  will  be  required.  Since  they  are  vehicle  dependent,  they  must  be 
supplied  by  the  user. 

3.2.1  Function  OXL  (Off-axis  Losses) 

This  routine  computes  transducer  pattern  attenuation  in  any  direction. 
The  call  sequence  is: 

VALUE  = OXL  (IFLAG,  COSA,  COSB,  COSC) 

where  IFLAG  is  0 for  receive  pattern,  1 for  transmit  and  COSA,  COSB,  and 
COSC  are  the  X,  Y,  and  Z direction  cosines  in  the  direction  of  interest. 

Other  values  necessary  for  the  computation,  e.g.  transducer  type,  frequency, 
sound  velocity,  etc.  may  be  supplied  via  a COMMON  statement.  The  value 
returned  is  a fraction  of  the  on-axis  intensity,  from  0.0  to  1.0. 

3.2.2  Function  RRF  (Reverberation  Rejection  Filter) 


band. 


This  routine  will  interpose  a filter  to  modify  the  energy  in  each 
The  call  sequence  is: 


-2- 


■IHi BBaBH 
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VALUE  = RRF  (FREQ) 
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where  FREQ  is  the  center  frequency  of  the  band  in  question.  The  value 
returned  is  the  fraction  of  energy  which  is  passed  by  the  filter  in  that 
band,  from  0.0  to  1.0. 

3.2.3  Function  TVG  (Time  Variable  Gain) 

This  routine  computes  receiver  gain  as  a function  of  time.  The 
call  sequence  is: 


GAIN  = TVGF  (TIME) 


Examples  of  OXL,  RRF,  and  TVGF  are  provided  in  the  program  listings. 
Appendix  G.  These  are  unrelated  to  any  real  system. 

3.2.4  Subroutine  SPRCMP  (Spreading  Computation) 

In  addition  to  the  above,  SPRCMP  may  be  provided  by  the  user  to 
generate  one  or  more  of  the  spreading  function  tables  when  the  program  is  run. 
The  subroutine  has  no  arguments. 

3.3  Input  Data 

Input  data  to  OOP  is  from  two  sources:  the  output  file  of  sorted 

ray  data  from  programs  SONAR  and  RAYSRT  (references  1 & 2)  and  input  data 
cards.  Besides  the  sorted  ray  data  returned  from  insonified  portions  of 
the  boundaries,  the  file  contains  eleven  parameters  passed  from  the  SONAR 
program. 

The  card  input  format  is  free- form  with  blanks  ignored.  Variables 
are  punched  in  fields  of  arbitrary  length,  separated  by  commas.  Data  may 
be  integer,  real  (including  a decimal  point)  or  alpha-numeric  (appearing 
between  single  quotes).  Note  that  no  check  is  made  for  the  appropriateness 
of  any  piece  of  data  to  any  name.  Neither  names  nor  numbers  may  be  split 
between  cards.  All  data  cards  appear  literally  in  the  printed  output.  Two 
kinds  of  data  may  appear  on  the  cards:  option  fields  and  data  fields.  Option 

fields  contain  only  the  name  of  the  option  to  be  invoked.  These  options  may 
modify  the  form  or  content  of  input  data,  computation,  and  output  data. 

Data  fields  consist  of  a variable  name  followed  by  an  equal  sign  and  one  or 
more  values,  as  appropriate,  separated  by  commas. 


( 


■ x 


For  option  or  data  names  longer  than  six  characters,  only  the 
first  six  characters  are  interpreted.  Additional  characters  may  be 
used  to  improve  readability,  but  are  ignored  by  the  input  routine.  All 
data  variables  are  initially  zero. 

3.3.1  Options 

Options  for  DOP  are  listed  below  for  input,  for  output,  and  for 
computation. 


l — T'  T'WM— KTflT  . 1 35^*1  ■ 
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meanings: 


below: 


3. 3. 1.1  For  Input 

For  input,  there  are  two  options,  which  have  the  following 

• GO  - Stop  reading  data  and  begin  computation. 

•NO  TAPE  - No  input  data  file  is  provided  and  no 
boundary  reverberation  will  be  computed. 

3.3. 1.2  For  Computation 

For  computation,  the  options  and  their  meanings  are  listed 


•CENTER  - Doppler  bands  are  computed  such  that  the 
transmitted  frequency  is  centered  in  one  of  the  bands 
instead  of  appearing  at  the  edge  of  a band. 

• END  - Stop  all  program  activity  and  exit. 

•FILTER  (Used  only  with  SPREAD  option  ) - Apply  filter 
to  the  spread  output  data. 

• KNOTS  - Compute  the  intensity  in  bands  of  equal 
apparent  range  rate  or  "knots  of  doppler”  instead 
of  equal  frequency  range. 

•NO  BOTTOM  - Do  not  compute  bottom  reverberation. 

•NO  SURFACE  - Do  not  compute  surface  reverberation. 

•NO  VOLUME  - Do  not  compute  volume  reverberation. 

•SPREAD  - Apply  spreading  function  to  output  data  and 
change  the  format  of  the  output  data  listing.  (See 
note  under  data  variable  BSPRED.) 

•TIME  COMPUTATION  - Compute  additional  values  of  time 
for  which  reverberation  is  to  be  determined.  The 
additional  values  of  time  are: 

•Ping  interval  (See  PING  below.) 

•1/2  At  (See  DELT  below.) 

•For  every  path  or  combination  of  paths  to  surface 
or  bottom  whose  earliest  arrival  time  is  t,  the 
additional  values  are:  t,  t+l/2At,  and  t-l/2At. 

•Seventeen  fixed  values  ranging  from  .01  to  2.0  seconds. 


-4- 
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• Additional  values  every  1/2  second  from  2 seconds 
to  the  value  of  the  PING  Interval,  NOTE:  The  total 
number  of  time  entries  read  from  cards  and  computed 
as  a result  of  the  use  of  this  option  Is  limited  to 
400. 

; 

•TVG  (used  only  with  SPREAD  option  ) - Apply  a time- 
varying  function  to  the  spread  output  data. 

3.3. 1.3  For  Output 

For  output,  the  options  and  their  meanings  are: 

•NO  PRINT  - Suppress  the  printed  output  of  doppler  data. 
Printing  of  Input  data  cards  Is  not  affected. 

•PLOT  - Write  a tape  of  doppler  data  for  use  by 
subsequent  programs. 

• RELATIVE  BANDS  - If  the  KNOTS  computation  option  has  not 
been  specified,  print  the  band  limits  In  kilohertz 
relative  to  the  source  frequency;  with  the  KNOTS 
computation  option,  print  band  limits  In  knots  of 
doppler  relative  to  vehicle  speed;  l.e.  zero  doppler 
(unspread)  is  returned  from  dead  ahead. 

•TOTALS  - Print  only  the  totals  of  reverberation  at 
each  requested  time. 

•(PRINT  EVERY  - Under  Data  Variables  also  modifies  output.) 

3.3.2  Data  Variables 

In  the  following  listing  of  data  variables  for  DOP,  an  asterisk 
(which  Is  NOT  part  of  the  name)  denotes  variables  whose  values  are 
normally  taken  from  the  input  data  tape.  If  the  same  variables  are  supplied 
on  cards,  however,  the  data  from  the  cards  would  be  used.  It  should  be 
noted  that  some  values  which  could  be  changed  are  Implicit  In  the  ray 
data  supplied  by  the  SONAR  program  and  that  changing  these  values  by  card 
would  be  meaningless  and  misleading.  Examples  are:  source  depth  (DO), 
sound  velocity  of  medium  at  the  source  depth  (CO)  and  bottom  depth  (DBTTM). 
(Please  note  that  the  second  character  of  CO,  DO,  and  FO  Is  the  numeral 
zero. ) 

ALPHC* 

The  product  of  attenuation  coefficient  and  sound  velocity  at 
source  depth.  Units  are  dB  per  second. 
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BSPRED 

Spreading  function  table  for  bottom.  Since  the  bottom  scatterers 
would  usually  be  considered  stationary,  this  table  would  usually  represent  the 
spectrum  of  the  transmitted  pulse.  NOTE:  The  appearance  in  the  data  deck 

of  BSPRED,  SSPRED,  or  VSPRED  makes  the  use  of  the  SPREAD  option  redundant 
and  unnecessary.  Each  table  is  entered  as  half  a symmetrical  table  of  an  odd 
number  of  entries.  First  value  is  the  proportion  of  energy  in  a band 
remaining  after  spreading.  Next  value  is  the  proportion  spread  to  the 
two  neighboring  bands,  etc.  Therefore,  twice  the  total  of  all  values 
entered  should  equal  one  plus  the  first  value.  Each  of  the  tables  has  a 
maximum  of  150  values.  Failure  to  include  all  three  tables  with  the 
SPREAD  option  causes  a call  to  SPRCMP  in  an  attempt  to  generate  the  missing 
table(s). 


BW1DTH 

Band  width  in  hertz  or  knots  as  appropriate,  based  on  the 
computation  option  selected. 

CO* 

Sound  velocity  in  yards  per  second  at  source  depth. 

DATE* 

Alpha-numeric  date,  maximum  of  two  machine  words. 

DBTTM* 

Bottom  depth  in  feet. 

DELT 

Effective  pulse  length.  At,  in  seconds.  Effective  pulse 
length  is  the  length  of  a square  pulse  with  the  same  energy  content 
as  the  pulse  of  interest,  which  may  not  be  a square  pulse. 

DO* 

Source  depth  in  feet. 

FO 

Transmit  frequency  in  kilohertz.  If  no  value  is  supplied, 

1 kilohertz  is  used.  Although  this  parameter  is  not  supplied  from  the 
input  tape,  it  is  implicit  in  the  ray  data  from  the  attenuation  values 
(and  the  spreading  loss  correction,  if  used). 


IT*.. 


* » ' ■'%*■*••• 


-6- 


TT_:_  - - . 


1 


r 


OD  52258 

IDC* 

Alpha-numeric  identification  as  on  the  "constant  card"  of 
the  SONAR  program. 

IDV* 

Alpha-numeric  identification  as  on  the  "semi -variable  card"  of 
the  SONAR  program. 

LOGMV* 

Volume  scattering  coefficient,  dB. 

NBEAM 

Coded  description  of  transducer  patterns,  if  more  than  one 
can  be  generated;  intended  for  use  by  subroutine  OXL. 

OMEGA 

Platform  turn  rate  in  degrees  per  second. 

PING* 

Interval  between  successive  transmits  in  seconds. 

PRINT  EVERY 

For  a value  n,  print  only  every  nth  doppler  band  at  each  time. 

Bands  to  be  printed  are  chosen  so  as  to  include  the  band  containing  the 
transmit  frequency  as  its  lower  bound  or  its  center,  if  CENTER  option  is 
used.  If  no  value  is  supplied,  1 is  used. 

PULSE 

Coded  description  of  pulse  shape,  if  several  options  exist. 

Intended  for  use  by  subroutine  SPRCMP. 

S* 

Source  level  in  dB  relative  to  1 yard. 

SSPRED 

Spreading  function  table  for  surface.  (See  note  under  BSPRED.) 

THTMAX 

Approximate  value  in  degrees  of  largest  angle  (between  velocity 
vector  and  sound  rays)  to  be  considered  in  computing  unspread  doppler 
bands.  If  no  value  is  supplied,  90  is  used. 


Values  of  elapsed  time,  In  seconds,  (measured  from  the  midpoint  of 
the  transmitted  pulse)  at  which  reverberation  Is  to  be  computed.  Maximum 
of  400  values. 

VS 

Vehicle  speed.  In  knots. 

VSPRED 

Spreading  function  table  for  volume.  (See  note  under  BSPRED). 

3.3.3  Sample  Data  Deck 

Figure  3-1  Illustrates  an  Input  card  deck  for  program  DOP.  It 
assumes  that  at  least  four  files  (or  complete  sets)  of  data  are  on  the 
Input  tape. 
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In  such  a case,  each  file  on  the  input  tape  (corresponding  to 
different  "semi -variable"  cards  in  the  original  SONAR  program  deck)  is 
processed  in  turn.  Options  and  data  variables  from  cards  are  retained 
until  over  written  from  the  tape  or  from  new  cards.  (Please  note, 
however,  that  options  once  invoked  or  turned  "on"  can  not  be  turned 
"off"  again.)  In  general,  the  program  proceeds  until  an  END  option  is 
read  from  a card,  or  until  the  last  data  file  on  the  tape  has  been 
processed,  whichever  occurs  first. 


The  indicated  cards  in  Figure  3-1  illustrate  certain  features  of 
the  program.  Cards  at  1 exemplify  the  continuation  of  a string  of  values 
on  several  cards.  Also,  note  that  of  the  six  values  of  time,  two  are 
repeated.  The  program  sorts  values  of  time  in  ascending  order  and 
eliminates  duplicates.  The  card  at  2 redefines  band  width  as  5.  knots 
although  it  had  been  set  to  .25  knots  on  the  second  card  of  the  deck. 

Card  at  3 marks  the  end  of  card  data  for  file  1 on  the  tape.  The  card  at 
4 indicates  that  all  the  same  data  is  to  be  used  for  File  2.  At  5 NBEAM 
is  used  because  the  hypothetical  torpedo  has  two  transmit  beams.  Some 
sonar  problems  might  never  use  this  feature.  At  6 invocation  of  FILTER 
will  cause  a call  to  require  RRF  to  modify  the  level  of  each  band. 

Without  this  option,  no  RRF  would  be  required.  Similarly,  if  one  of  the 
variables  at  7 was  not  supplied,  a subroutine  SPRCMP  would  have  been 
called  by  the  program.  The  card  at  8 replaces  aVj^  of  the  values  of  TIME 
previously  read.  Note  throughout,  the  free,  even  arbitrary  use  of  blanks 
which,  we  repeat,  are  ignored.  Also  note  that  values  used  on  these  cards 
are  arbitrary  and  unrelated  to  any  real  system. 


3.4  Program  Output 

Figures  3-2,  3-3,  and  3-4  illustrate  the  output  from  Files  1,  3,  and 
4 of  an  input  tape  used  with  the  sample  input  deck.  In  addition,  the  input 
deck,  up  to  the  GO  option,  or  between  successive  GO  options,  is  listed  on  a 
separate  page  before  each  output  data  page.  No  illustration  of  this  page  is 
provided. 

Figures  3-2,  3-3,  and  3-4  show  the  three  principal  formats  for 
the  output  data.  The  pages  are  mostly  self-explanatory.  One  point,  however, 
should  be  mentioned.  The  appearance  of  .00  dB  in  a column  seems  ambiguous, 
since  it  may  mean  zero  dB  or  zero  energy  (-  dB).  In  practice,  however, 
this  ambiguity  should  rarely  present  any  difficulty.  The  truly  empty  bands 
will  always  lie  at  the  top  or  the  bottom  of  a column  of  data.  Unless  the 
first  or  last  nonvoid  band  has  exactly  zero  dB,  which  can  usually  be  deter- 
mined from  the  band  and  column  totals,  leading  and  trailing  zeros  in  a column 
represent  zero  energy.  This  form  of  printout  was  chosen  in  lieu  of  printing 
an  arbitrary  large  negative  number,  because  It  made  the  page  less  cluttered. 

It  has  to  date  presented  no  problems  to  users  of  the  program. 

Three  data  decks,  to  programs  SONAR,  RAYSRT,  and  DOP  were  used  to 
generate  Figures  3-2,  3-3  and  3-4.  These  are  listed  in  Appendix  A.  It  will 
be  noticed  that  the  DOP  data  deck  is  not  quite  identical  to  that  in  Figure  3-1. 

3.5  Ancillary  Programs 

Besides  the  SONAR  and  RAYSRT  programs  which  produce  the  input  file  of 
sorted  ray  data,  two  other  programs  exist  which  may  be  of  help  to  the  user 
of  DOP.  These  are: 
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•DENFSP  - This  program  may  be  used  to  prepare  spreading 
function  tables  for  DOP.  It  will  combine  (convolve)  a 
spectrum  and  a density  function  and/or  convert  a spectrum 
Into  a density  function.  Print,  punch,  and  plot  options 
are  available.  Program  DENSFP  Is  documented  In  reference  (3). 

•SRNBT4  - This  uses  the  ‘plot1  output  tape  from  DOP.  The 
program  prints  or  plots  reverberation  in  doppler  bands. 

It  can  translate  the  data  Into  broader  bands  than  were 
computed.  It  can  also  plot  selected  bands  vs.  time,  and 
plot  or  punch  cards  for  enabling  level  for  use  by  the  SONAR 
program.  Program  SRNBT4  is  documented  in  reference  (3). 

3.6  Other  Considerations 

Although  standard  FORTRAN  IV  is  the  primary  language  of  DOP,  there 
are  two  UNIVAC  FORTRAN  V features  which  have  been  employed: 

•The  INCLUDE  Statement  - This  is  a convenient  way  of  adding 
blocks  of  cards  to  several  routines,  insuring  that  such 
blocks  will  be  identical  in  each  routine.  For  compilers 
without  this  or  an  analogous  feature,  it  is  a comparatively 
simple  matter  to  duplicate  the  blocks  of  cards  and  to  add 
them  physically  in  the  proper  places. 

•The  FLD  Function  - This  is  a bit-manipulating  function  which 
is  used  only  in  subroutine  INPUT  for  assembling  individual 
characters  into  variable  names.  For  compilers  without  a 
comparable  feature,  an  assembly-language  routine  must  be 
written. 

As  used  In  subroutine  INPUT,  the  two  constants  NWORD  and  NCHAR 
contain  the  number  of  bits  per  machine  word  and  the  number  of  bits  per  alpha- 
numeric character.  For  the  UNIVAC  1110,  these  numbers  are  36  and  6 respectively. 
They  can  be  changed  via  a data  statement  in  the  BLOCK  DATA  subroutine,  INBLK. 

SECTION  IV 

THEORETICAL  BASIS  AND  PROCEDURE 

4.1  General 

The  thinking  which  led  through  various  models  of  reverberation 
to  the  one  implemented  in  DOP  has  been  covered  extensively  in  earlier  Informal 
works  by  A.B.  Poynter.  Five  of  these  works  are  inlcuded  as  Appendices  B through  F. 
The  conventions,  assumptions,  and  simplifications  employed  In  the  program  are 
listed  below: 


•Both  the  surface  and  bottom  are  treated  as 
horizontal  planes. 

•The  platform's  velocity  vector  Is  horizontal. 
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•The  medium  varies  only  with  depth.  There  is  no 
f current. 


•All  times  are  measured  from  the  midpoint  of  the 
transmitted  pulse  (see  Appendices  C and  D). 

•No  allowance  is  made  for  differences  in  transmit  and 
receive  paths  caused  by  translation  of  the  platform 
(Appendix  E). 

•When  the  transmit  and  receive  paths  are  of  different 
types,  e.g.  one  direct  and  one  reflected,  the  boundary 
scattering  strength  between  paths  is  taken  as  the 
average  of  the  back  scattering  strengths  for  the  two 
single  paths.  This  is  one  of  the  least  defensible 
assumptions,  and  investigation  in  this  area  is  needed. 

4.2  Description  of  the  Models 

4.2.1  Boundary  Reverberation 

Refer  now  to  Figure  4-1.  A ray  is  considered  as  emanating  from 
a source  at  some  depth,  d0>  located  on  the  Z axis.  Since  sound  velocity 
is  a function  of  depth  only,  a ray  lies  entirely  in  a vertical  plane,  and 
the  locus  of  all  rays  leaving  the  source  at  an  initial  angle  0O  from  the 
horizontal  is  a surface  of  revolution.  If  these  rays  reach  the  surface. 
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FIGURE  4-2 

Two  different  ray-bundles,  called  "paths"  (direct  path  and  singly 
bottom- reflected  path)  insonifying  a common  portion  of  the  surface. 
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the  XY  plane,  their  points  of  intersection  describe  a circle.  For  con- 
venience, we  consider  only  rays  in  the  first  quadrant  of  the  XY  plane. 

The  SONAR  program  is  used  to  run  a fan  of  rays.  To  cover  all 
frequencies  seen  by  the  receiver,  these  would  range  from  vertical  upward 
to  vertical  downward  (In  any  case  they  should  cover  a range  equal  to  ± 
THTMAX).  The  points  where  these  rays  intersect  the  surface  and  bottom 
are  examined  by  program  RAYSRT.  The  points  are  separated  into  strings 
defining  annular  areas  of  the  boundaries,  which  are  insonified  over 
independent  ray  paths,  i.e.  direct  or  one  or  more  times  refracted 
or  reflected.  Each  such  string  of  points  is  hereinafter  referred  to  as 
a path. 


Figure  4-2  illustrates  portions  of  two  such  independent  paths 
at  the  surface  and  two  others  at  the  bottom.  In  practice,  for  completeness, 
all  four  paths  would  include  the  vertical  rays,  and  overlap  would  be 
complete  unless  interupted  by  the  refraction  of  some  ray  away  from  a 
boundary.  For  the  purposes  of  the  explanation  following,  however,  we  will 
assume  that  paths  A and  B are  illustrated  in  their  entireties. 

To  understand  the  operation  of  DOP  program,  refer  now  to  Figure  4-3. 
Let  the  dots  along  the  horizontal  axis  represent  points  at  which  rays  have 
reached  the  surface  directly  from  the  source,  called  path  A.  The  region 
so  insonified  lies  between  the  light  solid  lines.  Let  the  X's  represent 
rays  arriving  via  a single  reflection  from  the  bottom,  called  path  B. 

The  insonified  region  is  delimited  by  the  light  dashed  lines. 

Let  the  heavy  solid  line  represent  the  range  at  which  the  two-way 
travel  time  is  equal  to  a requested  data  time,  t,  having  been  transmitted, 
say,  over  path  A and  received  over  path  B.  The  heavy  dashed  line  represents 
the  ranges  at  which  the  two-way  times  are  t ;At  and  t + At,  the  time  limits 
due  to  the  pulse  length  At.  2 2 

Now  the  area  insonified  over  both  paths  and  lying  within  the  time 
limits  will  return  sound  to  the  source  at  time  t.  Figure  4-4  shows  this 
annular  area.  The  annulus  is  now  divided  into  the  incremental  areas  bounded 
on  two  sides  by  equal -time  lines  and  on  the  other  two  by  equal -frequency 
lines.  (Sometimes  one  of  these  lines  degenerates  into  a point).  Such 
areas  are  labeled  _a  and  la  in  the  figure.  The  co-ordinates  of  the  corners 
(numbered  1 to  4)  are  computed  and  the  areas  are  approximated  by  polygons, 
assuming  straight  lines  between  the  corners.  From  the  ray  data  are 
computed  the  average  transmission  loss  and  scattering  strength  over  the 
polygon.  These,  combined  with  average  transmit  and  receive  pattern  losses, 
permit  the  approximation  of  the  back-scattered  energy  from  each  area.  These 
contributions,  expressed  in  intensities,  are  summed  in  random  phase  within 
the  requested  time  and  band  limits. 
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FIGURE  4-4 

The  part  of  the  common  insonified  region  which  also  lies  between 
the  time  limits  returns  energy  to  the  source  at  time  t.  This  region 
has  been  broken  up  into  incremental  areas,  bounded  on  two  sides 
by  equal -time  lines  and  on  two  sides  by  equal -frequency  lines. 
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4.2.2  Volume  Reverberation 

Figure  4-5  Illustrates  the  computation  of  volume  reverberation.  For 
the  straight- running  case,  the  return  In  a band  at  a given  time,  t.  Is 
from  the  frustum  of  a hollow  cone  with  spherical  bases.  The  cone  1$ 
divided  radially  about  Its  face  forming  small  areas,  and  the  pattern 
losses  to  each  area  are  combined  to  yield  an  average  pattern  loss  for  the 
Iso-frequency  volume.  Volume  reverberation  Is  computed  using  this  pattern 
loss,  average  two  way  transmission  losses,  volume  of  the  hollow  frustum,  and 
volume  scattering  coefficient.  In  the  turning  case  the  Insonlfled  volume  Is 
not  a cone,  and  the  transmit  and  receive  directions  are  different.  However, 
the  basic  procedure  Is  the  same  as  for  the  non- turning  case. 


FIGURE  4-5 

Volume  Insonified  at  time  t between  two  equal -frequency  cones 
and  two  equal -time  spheres.  Pattern  losses  are  measured  to 
each  area  marked  on  the  end. 
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4.3  Equations 
4.3.1  Preliminary 

The  following  symbols  and  definitions  are  used  in  this  section: 

C = the  velocity  of  sound  in  water,  a function  of  water  depth 
f = the  frequency  of  a sound  wave 


t = the  elapsed  time  at  which  reverberation  is  to 
be  computed,  measured  from  the  midpoint  of  the 
transmitted  pulse 


V = the  speed  of  a platform 


Y=  the  angle  between  the  sound  direction  and  a platform 
velocity  vector 


\=  the  wave  length  in  water,  and 
uJ  = the  turn  rate  of  the  source 
Let  the  subscripts 

0 (zero)  = source, 

T = target,  and 
R = receiver 


Further,  let  the  superscript  designate  terms  associated  with 
sound  being  re-radiated  from  a target  as  opposed  to  the  unsuperscripted 
terms  associated  with  sound  approaching  a target. 


The  X axis  is  defined  as  the  platform  velocity  vector  at  time 
t = 0 for  volume  reverberation,  or  its  projection  on  the  surface  or  bottom 
for  boundary  reverberation. 


Using  this  symbology 


Ao  = 


Cn  ” v o Cos 


Approaching  a target 


-x0  = 


C r (C0  - V0  Cos  <0) 


Cn  fn 
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The  frequency  observed  by  a receiver  on  the  target,  and  that 
re- radiated  Into  the  water  will  be 


VT  Cos  yt 


Cp  (^T  ~ ^0S  Tj) 

CT  (Co  - Vo  Cos  Y0) 


The  wave  length  of  re-radiated  sound  Is 


v*  = c*  - vt  Cos  y»  = CT  (C0  - V0  Cos  Yn)  (C*  - VT  Cos  y*) 


approaching  the  receiver, 

= CRCT  (Co 


foCo  (Ct  ~ Vj  Cos  Y y) 


Cos  Y 0)  (Ct  - VT  Cos  Yt) 


C?  f0C0C?  (CT  - VT  Cos  yt) 

and  the  frequency  seen  by  the  receiver  is 

^ 3 Cr  ~ Vr  Cos  Yr  ^ ^ CqCt  (CT  - VT  Cos  yt)  (Cr  - Vr  Cos  Yr) 

^r  * 0 CRCT  (Co  - Vo  Cos  Yo)  (C?  - VT  Cos  y£) 

Simplifying,  in  the  vicinity  of  the  target,  C*T  * CT,  and 
considering  reverberation  as  being  from  motionless  scatterers  (their  motion 
is  handled  statistically  as  "spreading"  in  the  returned  sound)  VT  - o. 

Also,  we  are  interested  only  in  the  monostatic  case  (source  and  receiver 
at  the  same  point),  so  CR  = C0  and  VR  = V0.  Applying  these  identities  we 

have  r r „ 

f _ fo  C0  - V o Cos  Yr  (1) 

C0  - V0  Cos  y 0 
4.3.2  Boundary  Reverberation 

For  calculations  of  boundary  reverberation,  we  will  resolve  the 
angle  Y into  its  spherical  components,  the  horizontal  angle  <!>  and  the  vertical 
angle,  8.  The  equation  now  becomes 

^ Cq  — Vo  Cos  6r  Cos  $r 

C0  - V0  Cos  0q  Cos  $o 

Since  C is  a function  of  depth  only,  the  transmitted  and  received 
rays  are  in  the  same  vertical  plane.  For  the  straight-running  case, 

♦r  = w * ♦o  but  in  general  $R  ■ ir  + 60  - ut.  Finally,  then,  the  frequency 
of  reverberation  received  from  a point  on  a boundary, 

f * f Cq  + Vp  Cos  6r  Cos  ($p  - ut)  (2) 

C0  - V0  Cos  e0  Cos  *0 
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Refer  again  to  Figure  4-4.  In  order  to  be  certain  that  all 
incremental  areas  are  assigned  to  the  correct  frequency  bands.  It  Is 
Important  to  know  what  Iso-frequency  line  Is  tangent  to  an  Iso-tlme  line. 

This  Is  equivalent  to  finding  the  maximum  or  minimum  frequency  occurring 
at  a given  time.  It  is  clear  from  Figure  4-4  that  all  the  variables  of 
equation  (2)  except  t0  have  been  given  fixed  values.  Therefore,  we  find 
at  what  azimuth  angle  received  frequency  is  a maximum  or  minimum.  Dif- 
ferentiating equation  (2), 

df  R _ _f  V0  Cos  eR  Sin  (t0  - at) 
dt0  C0  - V0  Cos  e0  Cos  *0 

-f0  Lo*  v<>  Cos  e»  Cos  v0  Cos  6p  Sin  t0  = 0 

(C0  - V0  Cos  e p Cos  $p  )2 

Clearing  fractions  and  rearranging  terms, 

Vp  Cos  0p  Cos  0R  [sin  tp  Cos  (tp  - wt)  - Cos  tp  Sin  (tp  - wt)J 

+ C0  V0  [cos  0R  Sin  ( $p  - uit)  + Cos  et  Sin  to]  * 0 

Reducing  the  first  term  and  expanding  the  second  term  using  the  sin  (a-b) 
identity, 

V02  Cos  0O  Cos  0R  Sin  wt 

+ C0  V0  (Cos  eR  Sin  $p  Cos  wt  - Cos  eR  Cos  to  Sin  wt  + Cos  e0  Sin  to)  = 0 
Again  rearranging  terms 

(Cos  0R  Cos  wt  + Cos  0p ) Sin  t0  - Cos  eR  Sin  wt  Cos  t0 

+ Cos  0 p Cos  0 R Sin  wt  = 0 
Co 

Substituting 


i n 

r < 

r 

r 


we  have 


a = Cos  0R  Cos  wt  + Cos  0C 
b = Cos  eR  Sin  wt 


c = Cos  ft  Cos  0 Sin  wt 

Co 


a Sin  t«  - bCos  tn  + c = 0 


Squaring  and  substituting  for  sin2  t0  we  arrive  at  the  quadratic 
(a2  + b2)  Cos2  t0  - 2bc  Cos  t + c2  - a2  = o 


-22 


Set  into  the  quadratic  formula  and  simplified,  this  becomes  finally 

be  ± a /a2  + b2  - c2 


Cos  *o  = 


a2  + b5 
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(4) 


To  understand  the  significance  of  the  sign  of  the  radical, 
consider  that  when  u = 0,  sin  wt  3 0,  and  cos  <f>o*  ± That  Is,  the 
locus  of  the  points  of  tangency  between  iso- time  and  Iso-frequency  lines 
(the  min-max  locus)  is  the  X axis.  Also,  in  this  case  the  iso- frequency 
line  for  f = f0  is  the  Y axis,  with  greater  frequencies  to  the 
right,  lesser  frequencies  to  the  left.  Thus  it  is  seen  that  the  posi- 
tive and  negative  signs  define  4>0for  frequency  maximum  and  minima, 
respectively. 

We  will  now  find  the  co-ordinates  of  the  corners  of  the  Insoni- 
fied  polygon  labeled  1,  2,  4,  3 in  Figure  4-6.  At  time  t and  frequency 
f , from  equation  (2), 


f RC o - fRV0  Cos  e0  Cos  t0  = f0C0  + f0V0  Cos  eR  Cos  (t0  - wt) 
f0  Cos  eR  Cos  (to  - cut ) + fR  Cos  e0  Cos  t0  - f"Co  ~ focQ  = 0 


or 


or 


f0  Cos  0R  (Cos  ut  Cos  to  + Sin  ut  Sin  to)  + fR  Cos  eo  Cost0 

(f,  - U)  * o 
Vo 

(f0  Cos  eR  Cos  ut  + fR  Cos  e0)  Cos  t0  + (fo  Cos  eR  Sin  wt)  Sin  t0 

(f,  - fo)  - o 
Vo 


or 


Setting 


a = f0  Cos  eR  Cos  ut  + fR  Cos  e0  , 
b 3 f 0 Cos  eR  Sin  ut  , and 

c (f.  - fj 


we  have 

a Cos  t0  + b Sin  t0  - c = 0 

Squaring  and  substituting  for  sin2t0  as  previously,  we  obtain  the  quadratic 

(a2  + b2)  Cos  t - 2ac  Cos  t0  + c2  - b2  3 0, 
and  by  quadratic  formula 

ac  ± b >/a2  + b2  - c2 


Cos  t. 


a + b' 
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Similarly  it  can  be  shown  that 

be  1 a /a2  + b2  - c2 


Sin 


1rT~bT 


opposite  polarity  of  the  radicals  being  required  by  the  Sin  t2+  Cos  <P2  = 1 
identity. 


Now  from  Figure  4-6,  it  is  seen  that  for  point  1 

Cos  t0  = — and  Sin  *o  = — 

X,  X, 


where  x,  is  the  x value  associated  with  time  t,  i.e.  xf  is  the  radius  of 
the  iso-time  circle.  Therefore  for  point  1 


X,  = X, 


Y,  = X, 


ac 


+ b y/a2  + b2  - c2 


TTTP" 


and 


(5) 


be 


+ a ' Ja 2 + b2  - c2 

T + p 


From  considering  the  geometry  when  w = 0,  it  can  be  seen  that  the  upper  signs 
are  used  for  insonified  polygons  which  are  clockwise  of  the  min-max  locus 
when  f„>f0  . counter  clockwise  of  the  locus  when  fR<  f p , and  conversely 
for  the  lower  signs. 

The  area  of  the  polygon  in  Figure  4-6  is  clearly  the  sum  of 
trapezoids  2,  4,  8,  6 and  4,  3,  7,  8 less  the  sum  of  2,  1,  5,  6 and  1,  3, 

7,  5.  That  is. 


A = «*(X4  - X2)(V4  + Y2)  + MX3  - X4)(Y3  + Y4) 

- MX,  - Xa)(Y,  + Y2) - MX 3 - X,)(Y3  + Y,) 


or 


A = >s(X4Y4  + X4Y2  - X2Y4  - X2Y2  + X3Y3  + X3Y4  - X4Y3  - X4Y4 
- X j Y , - X,Y2  + X2Y,  + X2Y2  - X3Y3  - X3Y,  + X,Y3  + X,Y,) 
Consolidating,  we  arrive  at  the  computational  form  used  in  the  program: 

A s \h  [(X4Y2  + X3Y4  + X2Y,  + X,Y3)  - (X2Y4  + X4Y3  + X,Y2  + X3Y,)]|  (6) 
The  taking  of  absolute  value  avoids  negative  areas  in  some  quadrants. 

Equation  6 is  quite  general  and  computes  correct  areas  even  for 

triangles. 
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4.3.3  Volume  Reverberation 

Now  in  the  straight-running  case,  as  has  been  mentioned  before, 
the  iso-frequency  surface  is  a cone.  In  the  turning  case,  this  surface 
becomes  quite  complex,  even  discontinuous  for  certain  combinations  of  fre- 
quency and  turn  angle.  We  will  now  derive  a method  of  expressing  the  volume 
within  an  iso-frequency  surface  in  the  general  case. 

To  begin  with,  equation  (2)  is  equally  applicable  to  volume 
reverberation.  However,  in  the  simplified  model,  using  a uniform  unbounded 
medium,  the  transmit  and  receive  paths  will  necessarily  be  the  same  for 
the  straight  running  case.  In  any  case,  e R will  equal  e0.  Thus  the 
equation  becomes 

f _ f C0  + V0  Cos  e0  Cos  (<p0  - ut) 

C0  - V0  Cos  0Q  Cos  <P0 

clearing  the  fraction  and  rearranging  terms 

Cos  a0  [f0  Cos  (<P0  - ut)  + fR  Cos  *ol  = — 0 (fR  - fo) 

Vo 


Substituting 


Cos  (<p o - wt)  = Cos  <f>0  Cos  ut  + Sin  <t>0  Sin  wt. 
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Equation  (7)  describes  the  Iso-frequency  surface  in  cartesian  coordinates. 
For  simplicity  let  us  substitute 

a = f0  Cos  u»t  + fR  , 

b = f0  Sin  u»t  , and 


yielding 


aX  + bY  = c vX2  + Y 2 + Z2 


By  rotating  axes  through  some  angle  6 f the  Y term  of  the  above  equation 
can  be  eliminated.  Using  the  superscript  "prime"  to  denote  measurements 
in  a rotated  system  and  the  transformations 

x’  = X Cos  6 + Y Sin  6 and 


Y*  = Y Cos  6 - X Sin  6 


it  can  be  seen  that  to  eliminate  the  Y term  In  the  left  half  of  equation  (8),  6 
must  be  chosen  such  that 

a Sin  6 + b Cos  6=0  , or 

a2  Sin26  = b2  Cos 2 6 , or 


2 2 
Sin26  = b and  Cos2  5 = — -5 — . , 0r 
a2  + b2  a + b 
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/(X')2  + (Y')*  + Z2 


>/aJ  + b2 


Resubstituting, 


/(X')z  + (Y1)2  + Z2 


V°  (f.  - ft) 


v^f0  Cos  ut  + fR)2  + f02  Sin2  u.<t 


Vo(f»  - f0) 


✓(X1)2  + Of')2  + Z2 


/f02  + fB2  + 2fRf0  Cos  ut 


remembering  that  v/(X')2  + (Y1)2  + Z2  = C0t,  the  distance  that  the  sound 
has  traveled  in  time  t. 


But  the  left  hand  member  of  this  equation  is  cos  y .the  X direction 
cosine.  For  a particular  set  of  values  of  the  parameters  in  the  right-hand 
member  of  the  equation,  cosy1,  and  hence  Y ' are  constant.  Thus  it  can  be 
seen  that  for  each  time,  t,  the  intersection  of  the  iso-frequency  surface 
of  frequency  fR  and  the  sphere  of  radius  C0t  is  a circle  centered  on  the 
X'  axis,  which  is  at  angle  6 from  the  X axis.  This  circle  defines  a spherical 
segment. 


Let  us  define  our  differential  of  the  volume  enclosed  by  an 
iso- frequency  surface  as  the  area  of  this  segment  times  the  differential 
of  range  or 


dv  5 A dr 


or 


Now  the  area  of  a segment  with  thickness  of  h of  a sphere  of  radious  r is 
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The  volume  returning  energy  from  within  the  iso-frequency  surface 
at  time  t will  be 


V = 2*C3 


t + At 
2 


t - At 

T 


Ld_ 

Vo  (fn  - fo) 

>/f  o2  + fR2  + 2fRf0  Cos  wt 


dt 


(9) 


The  integral  in  equation  (9)  has  no  analytic  solution,  but  is 
evaluated  digitally.  The  volume  returning  energy  within  any  band,  obviously, 
would  be  the  difference  between  two  such  integrals  at  the  two  values  of  fR 
representing  the  band  limits. 


4.4  Verification 


Appendix  F gives  the  results  of  an  extensive  comparison  between 
various  models  for  computing  the  boundary  reverberation.  This  compares 
results  only  for  straight-running  cases.  In  addition,  comparisons  were  made 
in  1973  with  an  independent  model  developed  by  Dr.  J.  H.  Slaton  of  the 
Naval  Undersea  Center  for  a limited  number  of  turning  cases.  These  revealed 
several  difficulties  in  the  turning  case  which  have  been  corrected.  Agreement 
was  very  good,  considering  the  difference  in  mathematical  approach. 

For  the  straight-running  case,  comparisons  of  volume  reverberation 
with  the  method  described  in  reference  1 have  been  excellent.  This  was 
expected  because  of  the  similarity  in  assumptions  in  the  two  very  simplified 
models. 

SECTION  V 

AREAS  FOR  POSSIBLE  IMPROVEMENTS  AND  REFINEMENTS 

5.1  Volume  Reverberation 

The  simulation  of  volume  reverberation  includes  two  gross 
simplifications;  the  medium  is  homogeneous,  and  it  is  unbounded.  Thus, 
several  known  features  of  particular  environments  cannot  be  realistically 
approximated.  In  particular,  returns  from  scattering  layers,  returns  from 
regions  of  greater  or  lesser  focusing,  e.g.  caustics,  and  returns  over  multiple 
paths  are  ignored.  Still  it  is  hoped  that  the  model  can  supply  numbers  of 
the  right  order  of  magnitude  in  some  useful  cases. 

Almost  from  the  beginning,  the  three-dimensional  analogue  of  the 
boundary  calculations  was  recognized  as  more  realistic.  A preliminary  study 
was  made  in  1967  by  M.M.  Jacoby  of  NUC,  which  outlined  the  geometric  problems 
of  computing  incremental  volumes.  Although  precise  formulae  have  still  not 
been  developed,  the  problems  seem  mainly  those  of  bookkeeping;  that  is,  of 
ensuring  that  the  entire  insonified  volume  is  accounted  for.  The  calculations 
would  also  be  very  costly  in  computer  time. 

Implementations  of  this  model:  tedious  and  time-consuming  to 
write  and  check-out,  but  require  only  re-writing  of  one  DOP  subroutine, 
plus  the  changes  in  SONAR  mentioned  in  5.2. 
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5.2  Multi-path  Reverberation 

The  current  model  approximates  the  ray-path  losses  when  transmit 
and  receive  paths  are  not  the  same.  The  approximation  is  simply  the  average 
of  the  two-way  losses  over  each  single  path.  Now  the  transmitted  and 
reflected  rays  have  different  spreading  loss  even  over  the  same  (reciprocal) 
path.  The  difference  is  small,  however,  rarely  amounting  to  so  much  as  0.2  dB. 

Another  discrepancy  is  in  the  implicit  averaging  of  the  scattering 
strengths.  The  tacit  assumption  (not  to  be  dignified  as  an  approximation) 
is  that  scattering  strength  from  one  path  into  another  is  the  average  of  the 
back  scattering  strengths  for  the  two  paths.  In  the  absence  of  a good 
theoretical  or  empirical  alternative,  this  assumption  was  allowed  to  stand. 
There  is  a need  for  a single  expression  of  scattering  strength  in  all  dir- 
ections from  an  incident  ray. 

Incorporating  such  an  expression  in  DOP  would  require  that 
scattering  strength  be  removed  from  total  losses  passed  by  program  SONAR. 

Assesment  of  these  changes:  given  the  scattering  expression, 

the  latter  change  would  be  an  easy  matter;  the  former  would  not  be  much 
more  difficult. 

5.3  Vehicle  Translation  During  Ping 

Appendix  E discusses  the  consequences  of  using  the  same  ray-path 
for  transmit  and  receive,  ignoring  the  translation  of  the  vehicle.  This 
change  is  of  secondary  importance,  but  should  probably  be  included  for 
completeness  if  the  changes  in  5.2  were  implemented.  It  would  be  easy  under 
these  circumstances. 


5.4  “Spreading"  of  Reverberation 

There  are  theoretical  reasons  to  believe  that  doppler  shifts 
caused  by  motions  of  scatterers,  particularily  at  the  surface,  are  a function 
of  grazing  angle.  No  such  provision  is  made  in  this  model.  Incorporation  of 
this  feature  would  require  a fundamental  change  in  the  handling  of  boundary 
"spreading"  and  probably  minor  changes  in  SONAR  and  RAYSRT. 

Difficulty  of  these  changes:  moderately  lengthy,  but  not  complicated. 


5.5  Preferential  Orientation  of  Vehicle  with  Respect  to  the 
Environment. 

Real  environments  are  non-uniform  in  every  direction,  not  merely 
in  depth  as  currently  modeled. 

The  following  three  additions  to  the  model  would  take  account 
of  some  non-uniformities.  All  would  come  under  the  heading  of  major  revisions. 
In  addition,  it  is  clear  that,  conditions  being  different  at  all  points  and 
times  in  the  environment,  a single  ping  cycle  would  no  longer  give  a general 
description  of  reverberation.  The  changes  would,  however,  make  possible  the 
calculation  of  ranges  of  values  that  might  be  expected. 
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5.5.1  Current  and/or  Wind  Direction 

Only  OOP  would  need  revision.  Doppler  content  due  to  both  vehicle 
motion  and  scatterer  motion  Is  affected. 

5.5.2  Sloping  or  Broken  Bottom 

The  SONAR  program  can  compute  rays  reflected  from  any  bottom  made 
of  plane  segments  which  reflect  rays  in  the  vertical  plane  of  Incidence. 

Such  bottom  segments  must  be  perpendicular  to  the  plane  containing  the  rays, 
and  such  bottoms  are  of  limited  usefulness  In  modeling  real  environments. 

Spreading  losses  and  directions  of  reflections  can  be  computed 
for  rays  striking  any  analytically  describable  bottom.  However  a properly 
general  treatment  would  require  considerable  changes  not  only  In  SONAR  and 
DOP,  but  also  In  RAYSRT.  These  changes  would  seem  of  little  value  without 
also  the  ones  in  5.5.3. 

5.5.3  Three-dimensional  Velocity  Profiles 

The  use  of  velocity  profiles  varying  in  two  or  three  dimensions 
requires  completely  new  SONAR,  RAYSRT  and  DOP  programs;  no  mere  adaptation 
would  suffice.  While  ray  tracing  programs  already  exist  and  could  be  adapted, 
the  running  times  would  be  increased  by  an  order  of  magnitude,  at  least. 
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MEMORANDUM  P80203/ABP:pas 

18  January  1965 

From:  Code  P80203 
To:  Code  P802 

Subject:  Error  from  Using  10  log1Qr  as  a Measure  of  Effective  Train  Length 

When  Computing  Boundary  Reverberation. 

In  mathematical  models  commonly  In  use  for  computing  boundary 
reverberation  levels  In  dB,  a factor  of  10  log10T  is  added  In  to  account 
for  the  fact  that  an  Insonlfled  annulus  of  finite  width  returns  energy 
to  the  transducer  at  the  same  Instant  of  time.  In  this  context  t is  taken  to 
be  V(At)/2  where  V is  the  velocity  of  sound  and  At  Is  the  ping  length  In 
seconds.  It  is  rather  obvious  that  t is  an  accurate  measure  of  annulus  width 
only  when  the  transmission  path  Is  essentially  horizontal.  Since  analyses 
are  contemplated  for  applications  in  which  long  pings  and  steep  paths 
are  Involved,  It  seemed  desirable  to  Investigate  the  magnitude  of  errors 
which  might  accrue. 

Some  simple  calculations  were  made  after  the  fashion  indicated 
In  Figure  1,  assuming  an  isotropic  medium  (V  = 5000  ft/sec).  Starting  with 
a source  depth  and  an  initial  ray  angle  relative  to  the  horizontal,  the 
slant  ranae  (R  ) to  surface  intercept  was  computed  along  with  the  horizontal 
range  (X,)  to  this  point  and  the  two-way  travel  time.  If  time  is  measured 
from  the  end  of  transmission,  this  two-way  travel  time  (2t,)  would  identify 
the  instant  when  the  trailing  edge  of  the  wave  train  will  return  energy 
from  point  P,.  At  this  same  Instant  the  leading  edge  of  the  wave  train  will 
be  returning  energy  from  point  (P2)  for  which  the  horizontal  range  (X2)  can 
be  found  on  the  basis  of  a slant  range  (R2)  consistent  with  a two-way  travel 
time  of  2t2  - 2ti  + At.  The  width  of  the  insonifled  annulus  is  t*  = X2  - Xt. 

A measure  of  the  dB  error  involved  by  using  t instead  of  t*  is  given  by 
10  log10  f*. 

Such  calculations  were  carried  out  over  a range  of  values  for  Q<  of  20* 
through  80*  for  combinations  of  conditions  that  can  be  obtained  from  source 
depths  of  500,  1000,  2000,  and  5000  feet,  and  ping  lengths  of  .040,  .100, 
and  .250  seconds.  The  results  are  plotted  in  Figure  2.  The  error  as  plotted 
is  the  number  of  dB  that  should  be  added  to  values  computed  when  10  log  t is 
used.  It  appears  that  the  error  from  this  source  is  negligible  when  the 
significant  paths  are  less  than  20°  from  the  horizontal  as  is  the  case  for 
the  directional  systems  we  have  previously  analyzed  in  circular  search.  Even 
when  a refractive  environment  is  considered,  the  uncertainty  in  using  an 
average  velocity  in  computing  •*  should  not  increase  the  error  appreciably. 

As  expected,  the  Inaccuracy  mounts  rapidly  with  increasing  path  angles  above 
20”  and  with  increasing  source  depth.  However,  the  decrease  in  error  as  the 
ping  length  becomes  longer  was  not  anticipated.  This  trend  occurs  because 
in  the  ratio  4*  the  numerator  increases  less  rapidly  than  does  the  denominator 
with  increasing  ping  length. 
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If  the  annulus  width  itself  were  the  only  potential  source  of  error, 
it  would  be  possible  to  devise  a correction  procedure.  Unfortunately,  the 
contribution  of  each  unit  width  of  the  insonified  area  to  the  reverberation 
level  is  not  necessarily  the  same  since  the  transmission  and  vertical  pattern 
losses  may  vary  substantially  for  the  different  paths  involved.  Nor  in 
general  can  a total  error  be  computed  by  summing  the  independent  errors  for 
each  of  the  factors  handled  separately.  A method  for  assessing  the  composite 
error  will  be  presented  in  a subsequent  memorandum. 


A.  B.  POYNTER 
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APPENDIX  C 

MEMORANDUM  P80203/ABP:pas 

20  January  1965 

From:  Code  P80203 

To:  Code  P802 

Subject:  Review  of  Mathematical  Model  for  Computing  Boundary  Reverberation. 

References:  A.  NAVORD  Report  5606  (NOTS  1818)  - Analytical  Methods  for 
Predicting  Acoustic  Performance  of  Homing  Torpedoes  In 
Circular  Search. 

B.  TPR  246  (NOTS  TP  2498)  - Analytical  Studies  of  Sonar  in 
Refractive  Water. 

C.  P80203  Memo,  "Error  from  Using  10  log^T  as  a Measure  of 
Effective  Train  Length  when  Computing  Boundary  Reverberation," 
of  18  January  1965. 

In  the  analytical  work,  it  is  necessary  to  develop  mathematical  models 
of  the  phenomenon  to  be  studied  so  that  quantitative  results  can  be  computed. 
These  models  are  seldom  exact,  either  because  the  process  is  not  completely 
understood  or  because  of  a desire  to  avoid  complexities  which  do  not 
materially  affect  the  results.  The  models  may  be  adequate  for  the  purpose 
originally  intended,  but  it  is  dangerous  to  extend  their  use  to  new  situations 
without  reassessing  their  adequacy.  The  surface  and  bottom  reverberation 
equations  used  in  our  computer  programs  are  cases  in  point.  The  equations 
in  general  use  are  those  given  In  References  A and  B. 

For  boundary  reverberation  a ray  is  traced  to  the  surface  (or 
bottom)  and  the  expected  reverberation  level  at  the  time  in  the  ping  cycle 
corresponding  to  the  two-way  travel  time  is  computed,  using:  the  transmission 

loss  over  this  path;  vertical  pattern  losses  based  on  the  initial  ray  angle, 
and  a scattering  coefficient  per  unit  area  based  on  the  grazing  angle  at 
which  the  ray  strikes  the  surface.  Since  the  annular  ring  on  the  surface 
which  returns  energy  at  the  time  in  question  is  unlikely  to  be  unit  distance  wide, 
an  additional  term  (101og,oT)  is  added  ( T = half  the  train  length  = V(At)/2 
where  V is  the  velocity  of  sound  and  At  is  the  ping  duration  in  seconds). 

The  equation  was  formulated  in  the  context  of  fairly  shallow  systems  using 
relatively  short  pings  and  highly  directional  transducers  oriented  to  emphasize 
horizontal  coverage.  For  such  systems  the  equation  could  be  expected  to  yield 
acceptable  accuracy.  Future  applications  may  involve  one  or  more  of  the 
following: 

•Long  ping  lengths  for  correlation  detection 

•Broad  vertical  patterns  to  provide  greater  depth  coverage 

•Steep  paths  (deep  source  or  bottom  bounce) 

In  preparation  for  such  work  it  seems  advisable  to  re-examine  the  mathematical 
model  to  ascertain  whether  or  not  it  is  still  satisfactory. 
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In  Reference  C it  was  shown  that  10  log10r  may  not  adequately  account 
for  the  linear  extent  of  the  insonified  area  of  the  boundary  which  returns  en- 
ergy at  the  same  instant  of  time.  Moreover,  when  this  annulus  has  considerable 
width,  as  is  the  case  for  the  longer  pings,  the  contribution  of  various  parts  of 
the  insonified  area  to  the  total  energy  return  may  be  substantially  different 
because  of  variations  in  transmission  loss,  pattern  loss,  and  even  in  the 
scattering  coefficient  for  the  different  paths  involved.  To  better  account 
for  these  things  and  to  provide  a reference  for  assessing  the  accuracy  of  the 
model  now  in  use,  modifications  were  made  in  the  computational  procedures. 

The  key  to  the  method  is  the  temporary  assumption  that  10  log10T  = 0. 
This  is  equivalent  to  saying  that  a ping  length  is  chosen  so  that  the  insoni- 
fied annulus  at  the  boundary  is  1-yard  wide.  This  js  consistent  with  the  way 
in  which  the  scattering  coefficient  is  given  as  a function  of  grazing  angle, 
and,  over  such  a distance,  no  change  in  transmission  or  pattern  loss  need  be 
considered.  Otherwise,  the  desired  sets  of  environmental  conditions  and 
equipment  parameters  are  used.  A whole  family  of  rays  are  run  from  the 
selected  source  depth  to  the  surface,  and  reverberation  levels  are  computed  as 
before.  However,  instead  of  plotting  reverberation  level  versus  two-way 
travel  time  as  is  usual,  both  reverberation  and  travel  time  are  plotted 
against  the  horizontal  range  at  which  each  ray  strikes  the  surface.  If  the 
beginning  of  transmission  is  taken  as  zero  time,  any  selected  time  (2ti)  on 
the  two-way  time  curve  will  identify  the  horizontal  range  at  the  surface  from 
which  the  leading  edge  of  the  wave  train  is  returning  reverberation.  It 
follows  that  2t i - At  identifies  the  range  from  which  the  trailing  edge  is 
returning  energy  at  the  same  instant.  Integrating  the  reverberation  curve 
between  these  two  limits  should  give  the  actual  reverberation  level  quite 
accurately.  It  should  be  noted  that  the  values  in  the  one-yard  increments 
must  be  converted  from  dB  to  intensity  before  summing  and  then  reconverted. 
Reasonable  accuracy  can  be  achieved  with  much  less  work  by  using  average 
values  for  range  increments  several  yards  wide  and  multiplying  each  average 
value  by  the  number  of  yards  in  the  increment  before  summing. 

After  sufficient  points  have  been  determined  in  this  manner,  the 
values  of  reverberation  can  be  replotted  against  elapsed  time  to  yield  the  type 
of  presentation  generally  desired.  To  obtain  comparable  results  from  the 
original  model,  one  can  rerun  the  computer  program  using  the  appropriate 
value  for  10  log10Tas  determined  by  the  ping  length  of  the  transmission. 

It  is  more  efficient,  however,  to  return  to  the  working  curves  and  add  10 log  10 T 
to  the  values  of  the  reverberation  curve  for  t = l at  the  range  determined  by 
the  travel  time  curve  at  the  elapsed  times  for  which  data  are  desired. 

Surface  reverberation  levels  were  computed  by  the  two  methods 
using  parameters  for  the  Torpedo  MK  46  Mod  0 in  circular  search.  Attenuation 
by  the  RRF  was  not  considered.  The  source  depth  was  assumed  to  be  750  feet. 

The  velocity  profile  for  the  environment  is  shown  in  Figure  1A.  The  new 
computer  program  which  permits  use  of  a continuous  gradient  profile  was  used. 

The  comparative  results  for  a 40-ms  ping  are  shown  in  Figure  2.  The  vertical 
patterns,  of  course,  dominate  the  shape  of  the  curves  at  times  less  than 
one  second.  The  integration  procedure  would  normally  be  expected  to  lower  the 
peaks  and  fill  in  the  troughs  as  well  as  to  shift  them  slightly.  The  fact 
that  the  peaks  in  this  case  are  lower  for  the  curve  obtained  via  the  old 
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model  suggests  that  the  10  log  t correction  is  insufficient,  particularly 
at  the  shorter  times  which  would  be  identified  with  steeper  paths.  It  might 
be  kept  in  mind  that  the  curve  obtained  by  means  of  the  original  equation 
has  the  same  shape  as  would  the  10  log  t = o curve.  The  initial  ray  angles 
associated  with  a few  points  on  this  curve  are  shown  as  a matter  of  interest. 

Figure  3 shows  the  comparative  reverberation  levels  that  would 
obtain  if  a ping  length  of  250-ms  were  used.  It  is  clear  that  one  must  be 
very  careful  in  computing  reverberation  levels  when  long  pings  are  used; 
the  two  curves  are  substantially  different.  The  integrated  curve  was 
not  determined  at  times  less  than  .55  seconds  because  before  that  time 
all  of  the  ping  has  not  reached  the  surface. 

For  both  ping  lengths,  the  curves  would  match  more  closely  if  tl 
curves  obtained  by  means  of  the  original  model  were  advanced  in  time  by 
one-half  the  ping  length.  This  follows  since  the  computed  ray  path  would 
then,  in  effect,  intercept  the  surface  near  the  middle  of  the  insonified 
area  and  the  transmission  loss,  vertical  pattern  loss,  and  scattering 
coefficient  would  be  more  representative  of  the  area  as  a whole  than  is 
the  case  when  the  ray  falls  at  one  end  of  the  insonified  area.  Figures  4 
and  5 compare  the  results  when  this  subterfuge  is  used  for  the  .040  and 
.250  ping  lengths,  respectively. 

As  a further  check,  additional  calculations  were  made  for  a non- 
turning vehicle  at  500  feet  in  an  environment  characterized  by  velocity 
Profile  B shown  in  Figure  1.  Our  old  computer  program  based  on  linear- 

f rad lent  layers  was  used  in  this  instance.  Other  parameters  were  unchanged, 
he  results  are  similar  to  those  from  the  first  set  of  calculations. 

Figure  6 compares  the  40-ms  ping  data  from  the  integrated  method  with  that 
from  the  original  model  after  the  latter  had  been  advanced  in  time  by  half 
the  ping  length. 


Figure  7 presents  similar  data  for  the  250-ms  ping.  Part  A,  on  the 
left  shows  the  results  out  to  1.4  seconds.  Calculations  for  this  case  were 
carried  out  to  about  7.5  seconds,  and  the  curves  (after  translation)  remain 
in  close  agreement  until  about  6.7  seconds.  For  this  elapsed  time  the  paths 
involved  are  from  rays  that  start  out  below  the  horizontal  and  are  refracted 
back  toward  the  surface  by  the  positive  velocity  gradients  which  characterize 
this  environment  below  170  feet.  The  ray  having  a two-way  travel  time  of 
6.705  seconds  and  subsequent  ones  having  slightly  steeper  initial  angles, 
reverse  at  depths  where  the  velocity  gradient  becomes  less  positive.  This 
results  in  an  abrupt  increase  in  transmission  loss  and  a resulting  drop  in 
reverberation  level  as  indicated  in  Figure  7B  (recall  that  the  curve  based 
on  the  original  model  has  been  advanced  by  0.125  seconds.)  In  the  vicinity 
of  such  anomalies  the  integrated  method  offers  the  only  approach  for 
computing  expected  values  of  boundary  reverberation  levels  with  reasonable 
accuracy. 

Some  general  observations  appear  warranted: 

•Within  the  limits  imposed  by  the  input  data,  accurate 
expected  levels  of  surface  reverberation  as  a function 
of  elapsed  time  can  be  computed  by  means  of  ray  tracing 
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programs  by  summing  the  contributions  of  1-yard  increments 
of  the  insonified  annulus  which  contribute  at  the 
same  instant. 

•At  present,  this  can  be  accomplished  manually,  but 
it  may  be  possible  to  program  the  computer  to  this 
end. 

•The  present  mathematical  model  can  yield  results  which 
are  reasonably  accurate  over  the  elapsed  time  interval 
during  which  the  sound  paths  are  such  that  10  logioT 
is  a good  measure  of  the  width  of  the  insonified  annulus 
and  the  change  in  transmission  loss,  pattern  losses,  etc. 
over  this  annulus  is  not  appreciable. 

•The  accuracy  attainable  and/or  the  time  interval  over 
which  a given  accuracy  can  be  maintained  will  be  in- 
creased if,  in  effect,  the  elapsed  time  is  measured  from 
the  middle  of  the  transmittal  ping.  If  time  is  measured 
from  the  beginning  of  transmission  the  computed  two-way 
travel  time  to  the  surface  for  each  ray  should  be  in- 
creased by  At/2  to  obtain  the  time  consistent  with  the 
computed  reverberation  level. 

•Even  if  the  technique  in  the  above  paragraph  is  used,  errors 
may  become  substantial  when: 

•The  paths  are  very  steep  so  that  10  logtoT  is 
not  a good  measure  of  annulus  width; 

•The  insonified  annulus  is  wide  enough  to  encompass 
a spread  of  paths  through  the  minor  lobes  of  the 
transducer; 

•In  the  vicinity  of  transmission  loss  anomalies 
(either  substantial  focusing  or  defocusing) 

The  above  observations  should  also  apply  to  bottom  reverberation, 
except  that  additional  complications  may  arise  if  the  bottom  is  irregular. 


A.  B.  POYNTER 
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APPENDIX  D 

MEMORANDUM  P80203/ABP : pas 

25  January  1965 

From:  Code  P80203 
To:  Code  P802 

Subject:  Review  of  Mathematical  Model  for  Computing  the  Expected  Level 

of  Volume  Reverberation. 

Reference:  A.  P80203  Memo  "Review  of  Mathematical  Model  for  Computing 
Boundary  Reverberation,"  of  20  January  1965. 

Reference  A discussed  the  problems  associated  with  the  mathematical 
model  for  computing  boundary  reverberation.  It  was  pointed  out  that  the  model 
currently  In  use  In  our  computer  programs  generally  gives  reasonable  accuracy 
over  a large  portion  of  the  ping  cycle  provided  time  Is  measured  from  the 
middle  of  the  transmittal  ping.  The  errors  may  still  be  appreciable  over  the 
time  intervals  during  which  the  paths  to  the  surface  are  quite  steep  and  involve 
transmission  through  the  minor  lobe  structure  of  the  transducer,  and/or  exhibit 
anomalous  transmission  loss.  A method  is  delineated  by  which  errors  can  be 
reduced  by  summing  the  contributions  from  narrow  Increments  of  the  insonified 
annulus. 


The  model  currently  in  use  for  calculating  volume  reverberation  in 
dB  again  involves  the  addition  of  a factor  10  log10T=  10  log  V (At/2)  to 
account  for  the  lateral  extent  of  the  insonified  region  which  returns  energy 
to  the  transducer  at  the  same  instant  of  time.  (V  is  the  velocity  of  sound 
and  At  is  the  ping  duration  in  seconds.)  Since  for  volume  reverberation  the 
insonified  region  is  a spherical  shell,  the  lateral  extent  is  measured 
radially.  Thus  10  log  lo^is  a reasonably  accurate  measure  of  the  thickness  of 
this  shell  providing  the  value  chosen  for  V is  a good  average  value  for  the 
environment  being  considered.  Since  in  general  the  velocity  is  not  constant 
throughout  the  volume  of  interest,  there  are  resulting  errors.  However,  the 
percentage  variation  in  velocity  is  so  slight  that  the  errors  from  this  source 
can  be  expected  to  be  less  than  0.2  dB,  an  insignificant  amount  in  view  of 
our  inaccurate  knowledge  of  the  scattering  coefficients  which  obtain. 

A somewhat  more  significant  error  can  result  from  the  variation  in 
transmission  loss  over  the  shell  width,  particularly  for  long  pings.  To 
obtain  a measure  of  the  Inaccuracy,  calculations  of  volume  reverberation  were 
made  for  two  postulated  ping  lengths  in  an  environment  with  an  average 
velocity  taken  to  be  1660  yd/sec.  The  levels  computed  with  values  of  101og10r 
appropriate  for  the  selected  ping  lengths  were  compared  with  those  obtained 
by  a modification  of  the  integration  method  described  in  Reference  A with 
10  log,qTset  at  zero.  The  levels  were  plotted  as  functions  of  elapsed  time 
measured  from  the  beginning  of  the  transmission.  The  10  log10T  3 0 curve 
can  be  integrated  over  the  ping  length  In  seconds  except  that  the  resulting 
intensity  is  multiplied  by  a factor  of  0.83  before  reconverting  to  dB  to 
account  for  the  fact  that  one  millisecond  is  equivalent  to  an  0.83-yard  shell 
thickness  when  V has  the  value  previously  selected. 
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The  results  for  ping  lengths  of  0.040  and  0.250  seconds  are  shown 
In  Figure  1.  As  was  the  case  for  boundary  reverberation,  the  corresponding 
curves  from  the  two  methods  would  be  in  much  better  agreement  if  the  one 
for  the  integrated  method  were  shifted  back  in  time  by  half  the  ping  length. 
This  is  equivalent  to  measuring  elapsed  time  from  the  middle  of  the  trans- 
mission. When  this  is  done,  there  is  still  a residual  error  at  the  very 
short  times  due  to  the  fact  that  the  spreading  loss  at  short  ranges  is 
changing  rapidly.  The  integration  method  accounts  for  the  fact  that  the 
latter  portion  of  the  ping  is  contributing  a considerably  greater  portion 
of  energy  to  the  reverberation  level.  The  errors  for  short  pings  do  not 
seem  large  enough  to  justify  the  additional  work  of  using  the  integration 
method. 

Our  analyses  in  the  past  have  involved  systems  with  ping  lengths  of 
.040  seconds  or  shorter.  The  elapsed  time  for  both  interference  and  echo 
levels  was  taken  to  be  the  two-way  travel  time  for  the  paths  being  computed. 
For  the  short  pings  this  is  roughly  equivalent  to  measuring  time  from  the 
transmission  mid-point.  Consequently,  the  results  of  the  analysis  can  be 
considered  valid  insofar  as  the  factors  discussed  here  are  concerned. 

There  is,  however,  a basic  weakness  in  our  volume  reverberation 
model  which  may  be  really  serious,  especially  if  performance  analyses  of 
long-range  systems  are  attempted.  Moreover,  there  is  no  means  at  hand  for 
ascertaining  what  errors  might  be  involved.  This  situation  arises  because 
the  present  model  for  volume  reverberation  is  based  on  the  assumption  of 
an  Insonified  shell  over  expanding  In  time  in  an  infinite,  uniform  medium. 

It  is  well  known  that  an  ocean  is  neither  infinite  in  extent  nor  uniform  in 
its  acoustic  properties.  From  whatever  transducer  position  the  sound  paths 
will,  at  one  time  or  another,  intercept  the  surface  or  bottom,  and  variation 
in  the  physical  properties  of  the  water  results  in  differences  in  trans- 
mission loss  depending  on  which  path  is  considered.  Furthermore,  the 
volume  scattering  coefficient  in  general  is  not  constant  over  any  sizable 
volume  of  water,  and  is  often  observed  to  be  a rather  strong  function  of 
depth  and  time  of  day  (deep  scattering  layer  migration).  To  further 
complicate  the  situation,  the  paths  in  various  directions  are  weighted 
in  accordance  with  the  three-dimensional  transducer  patterns. 

The  difficulty  of  realistically  delineating  the  distribution  of 
acoustic  properties  over  a large  volume  coupled  with  the  even  greater  dif- 
ficulty of  integrating  the  return  from  the  entire  spherical  shell  has  no 
doubt  led  to  the  common  use  of  the  simplified  model.  Its  use  can  be  ration- 
alized as  follows.  While  it  is  true  that  the  transmission  loss  to  various 
portions  of  the  Insonified  shell  may  be  substantially  different  due  to 
refraction,  the  integration  process  will  average  out  such  effects.  Similarly 
the  variation  in  scattering  coefficient  will  average  out  in  the  integration 
so  that  selection  of  a constant  value  consistent  with  the  volume  being 
considered  should  suffice.  When  the  insonified  spherical  shell  is  truncated 
by  a boundary,  energy  aside  from  the  backscattered  as  boundary  reverberation 
or  lost  in  the  reflection  mechanism  will  be  returned  to  the  water  to  insoni- 
fy  volume  scatterers.  If  one  can  assume  that  boundary  losses  are  offset  by 
the  multiple  paths  that  are  set  up,  the  truncating  of  the  sphere  can  be 
Ignored. 
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There  is  little  doubt  that  in  a qualitative  sense  these  trends 
can  be  expected.  There  is  no  available  evidence,  however,  that  they  quan- 
titatively balance  out  within  acceptable  limits.  Faith  that  they  will  do 
so  deteriorates  further  when  one  considers  that  directional  transducer 
patterns  may  weight  the  return  from  anomalous  regions  so  that  their  con- 
tribution is  grossly  out  of  proportion  to  their  share  of  the  total  volume. 

The  depth  and  orientation  of  the  transducer  as  well  as  the  patterns  them- 
selves can  be  significant  factors.  For  example,  if  there  is  a pronounced 
deep  scattering  layer,  the  behavior  of  the  volume  reverberation  level 
with  elapsed  time  can  be  expected  to  be  quite  different  when  a transducer 
is  oriented  vertically  than  when  it  is  oriented  horizontally. 

One  can  also  argue  that  the  need  for  a more  complex  and  realistic 
model  for  estimating  volume  reverberation  is  not  great  on  the  basis  that  most 
systems  have  TVG  thresholds  which  protect  them  against  volume  reverberation 
at  short  range,  boundary  reverberation  is  likely  to  exceed  volume  reverber- 
ation at  intermediate  ranges,  and  noise  will  be  the  dominant  interference  at 
long  ranges.  While  this  reasoning  may  be  valid  for  many  (if  not  most) 
situations  analyzed  to  date,  there  is  no  assurance  that  it  will  always 
hold.  Volume  scattering  coefficients  above  average  when  coupled  with 
long  pings  may  raise  this  type  of  interference  to  new  prominence  in  some 
important  situations.  The  point  is,  no  one  can  be  sure  that  the  simple 
model  is  adequate  in  new  applications  unless  there  is  a comprehensive  model 
available  with  which  to  compare  results.  Efforts  by  able  mathematical 
physicists  to  develop  a more  realistic  approach  to  the  problem  of  estimat- 
ing volume  reverberation  should  be  solicited.  A companion  problem  of  great 
difficulty  is  the  development  of  a practical  method  for  estimating  the 
doppler  spectrum  of  both  boundary  and  volume  reverberation.  Such  infor- 
mation is  needed  for  assessing  performance  of  systems  which  employ  doppler 
discrimination. 

In  the  meantime,  the  best  that  one  can  do  is  to  continue  using 
the  present  method,  taking  care  that  the  elapsed  time  is  measured  from  the 
middle  of  the  transmitted  ping  and  then  adjusted,  if  desired,  to  the  same 
time  base  used  in  computing  boundary  reverberation  and  echo  level.  It  is 
clear  that  the  expected  values  of  both  the  echo  level  and  threshold  must  be 
figured  on  the  same  time  base  for  valid  results.  It  is  worth  noting  that 
for  elongated  echoes  (long  pings  or  other  reasons)  the  applicable  threshold 
level  may  vary  appreciably  depending  on  which  portion  of  the  echo  leads  to 
the  detection. 


A.  B.  POYNTER 
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MEMORANDUM 

From:  Code  P80203 
To:  Code  P802 


P80203/ABP:slh 
24  November  1965 


Subject:  Assessment  of  the  Effect  of  Platform  Translation  on  the  Area  Returning 

Boundary  Reverberation  at  a Given  Time. 

Math  models  for  computing  expected  levels  of  reverberation  nearly 
always  include  the  concept  that  the  return  at  any  given  time  comes  from 
scatterers  contained  in  a sphere-like  shell.  The  outer  and  inner  surfaces 
are  determined  by  the  sound  velocity  multiplied  by  one-way  travel  times 
(equal  to  one-half  the  elapsed  time  in  the  ping  cycle)  plus  and  minus 
one-fourth  of  the  ping  duration.  This  concept  fits  volume  reverberation 
at  least  for  times  less  that  that  required  for  the  expanding  shell  to  be 
truncated  by  the  surface  or  bottom.  The  area  returning  boundary  reverberation 
is  considered  to  be  bounded  by  the  traces  of  the  shell  boundaries  on  the 
surface  or  bottom.  These  traces  will  be  concentric  circles  if  the  shell  is 
truly  spherical. 

The  approach  is  clearly  realistic  if  the  transducer  is  stationary 
and  the  sound  velocity  is  a function  only  of  depth.  It  is  desirable  to  have 
some  quantitative  indication  of  the  error  introduced  if  the  transducer  is 
mounted  on  a moving  platform.  For  practical  cases  it  can  be  assumed  that  the 
velocity  of  the  platform  will  be  small  in  comparison  with  the  sound  velocity. 

Assume  that  the  platform  is  moving  at  a constant  velocity  (V)  directed 
horizontally  along  the  X-axis  in  an  oceanographic  coordinate  system.  The 
X-axis  is  taken  to  be  at  platform  depth.  Since  the  velocity  of  sound  in 
the  whole  ocean  really  varies  over  a range  of  only  about  10%  it  is  reasonable 
to  assume,  as  a first  approximation  in  determining  range,  that  the  sound 
velocity  (C)  is  a constant  in  a limited  volume  of  water.  Consider  the  X-Z 
plane  in  the  following  working  figure: 


Let  q be  the  position  of  the  platform  at  transmission  of  a particular  part 
of  the  ping  and  r be  the  platform  position  when  return  is  received  from  a 
point  (p)  after  some  elapsed  (t).  We  know  that  the  distance  traveled 
over  qpr  = tC  = qp'r  3 qp"r.  The  locus  of  all  p in  this  plane  is  an  ellipse 
with  foci  at  q and  r.  If  the  origin  is  placed  half-way  between  q and  r,  the 
equation  of  the  ellipse  is  x2/a2  + z2/b2  3 1. 
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When  the  point  is  on  the  X-axis  (at  p'),  z = 0 and  a = op'  = Ct/2. 
It  is  also  known  that  qo  = or  = Vt/2.  When  the  point  is  at  p",  X * 0,  and 
qp"  = p"r  = Ct/2.  Therefore, 

b = (Ct/2)2  - (Vt/2)2  = (t/2)  v/c2  - V2 

The  distance  from  either  focus  to  the  origin  should  be: 

>/C  2t  2/4  - (C2t2/4  - V2 12/4)  = /V7t2/4  = Vt/2 

which  checks  with  what  was  known  directly  from  the  movement  of  the  platform. 


The  equation  of  this  ellipse  is 


cV/4 


t /4(C2  - V2) 


The  eccentricity  of  such  an  ellipse  is  e = /(a2  - b2)/a2  = V/C.  Now  C is 

of  the  order  of  5000  ft/sec  while  a V of  15  knots  is  approximately  25  ft/sec. 
In  this  case  e = 25/5000  = .005.  For  a 45-knot  V,  e would  be  three  times 
larger  or  about  .015.  Therefore,  circular  assumption  is  not  bad. 

If  excursions  in  the  Y direction  are  allowed,  the  insonified  area 
would  be  the  surface  of  an  ellipsoid 

- x»  + v2  + z2  = l 

C2t 2/4  (t2/4) (C2  - V2)  (t2/4)(C 2 - V2) 

The  trace  of  this  ellipsoid  intersecting  a horizontal  plane  (e.g.  surface  or 
bottom)  would  be  an  ellipse  with  the  same  eccentricity  previously  computed. 
The  equation  of  such  an  ellipse  may  be  found  by  substituting  ^d),  the 
depth  differential  between  the  platform  and  the  boundary  in  question,  for  Z 
in  the  equation  above.  The  new  X - Y plane  now  coincides  with  the  surface  or 
bottom,  and  the  equation  for  the  ellipse  can  be  expressed  in  the  form 

(C2  - V2)X2  + C2Y2  = »s(C2t2)(C2  - V2)  - C2  (Ad)2.  (1) 

Let  us  consider  a finite  length  ping  of  duration  (At).  Since 
the  ping  is  not  emitted  all  at  once,  some  event  must  be  chosen  as  time 
zero.  This  is  taken  to  be  midpoint  of  transmission.  The  boundary  rever- 
beration at  some  time  (T)  would  come  from  an  annulus  bounded  by  two  non- 
concentric  ellipses  (if  the  platform  is  moving). 

A plot  of  vehicle  positions  at  initial  times  would  show: 


where  A is  the  position  at  transmission  of  the  leading  edge  of  ping;  B is  the 

position  at  t = 0;  C is  the  position  at  transmission  of  the  trailing  edge; 

and  D is  the  position  when  reverberation  is  returned  at  time  T.  The 
distances  A to  C = (At)V;  A to  D V(T  +At/2);  and  C to  D = V(T  - At/2).  By 

alternately  substituting  T + At/2  and  T - At/2  for  the  t in  Eq.  (1),  one 

finds  the  equation  for  the  two  ellipses  when  the  origin  for  each  is  half- 
way between  their  respective  foci. 
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The  origin  of  the  ellipse  stemming  from  the  trailing  edge  of  the 
ping  will  be  located  V(At/2)  further  along  the  X-axis  than  the  one  associated 
with  the  leading  edge.  Choosing  an  origin  appropriate  to  the  midpoint  of  the 
ping  and  substituting  in  Eq.  (1),  one  obtains  the  follov'ing  equations: 


Leading  edge:  (C2  - V2)  [x  + V(At/4)]2  + cV  = [c2(T  + At/2)2/4](C2  - V2)  - C2(Ad)2 

Trailing  edge:  (C2  - V2)  [x  - V(At/4)]2  + C2Y2  = [c2(T  - t/2)2/4]  (C2  - V2)  - C2(Ad)2 

When  y is  0 and  letting  subscripts  L and  T designate  the  leading  and  trailing 
edges,  respectively: 

XL  = ± vtc2/4) (T  + At/2)2  - C 2 (Ad ) 2/ ( C 2 - V2)  - V(At/4) 

XT  = ± >/(C2/4)(T  - At/2)z  - C2(Ad) 2(C  2 - V2)  + V(At/4) 

Using  primes  to  differentiate  between  concepts,  the  similar  equations  for 
the  concentric  circle  approximation  with  the  origin  at  the  same  place  (half- 
way between  the  midpoint  of  transmission  and  the  point  of  reception)  would 
be: 

XJ^  = + /(C2/4)(T  + At/2)2  - (Ad)2 
and  x;  = ± /(C2/4)(T  - t/2)2  - (Ad)2 


It  is  clear  that  XL  < X’Lfirst,  by  the  quantity  V(At/4),  and 
then  by  the  amount  the  evaluation  of  the  square  root  term  is  reduced  because 
(Ad)2  is  multiplied  by  C2/(C2  - V2)  rather  than  by  C2.  Now  V( At/4)  will  be 
small  in  practical  situations.  For  a platform  velocity  of  45  knots  and  at 
At  of  0.25  seconds,  V(At/4)«=  5 feet.  The  other  term  is  more  difficult 
to  evaluate  since  it  depends  on  so  many  factors.  At  V = 45  knots,  C2/(C2  - V2) 
is  the  order  of  1.00025.  However,  the  extent  to  which  this  factor  modifies 
XL  depends  not  only  on  Ad,  but  also  on  At  and  the  value  of  T for  which  the 
evaluation  is  desired.  At  the  times  for  which  the  evaluation  is  important, 
the  first  term  under  the  square  root  will  be  much  larger  than  the  second,  so 
the  difference  in  XLand  X'L  will  be  small. 

In  the  forward  direction  the  differences  in  XT  and  X'T  (when  Y * 0) 
is  less  than  in  the  case  of  the  leading  edge  since  the  correction  for 
translating  the  origin  from  the  center  of  the  ellipse  to  the  center  of 
the  circle  changes  sign.  The  two  factors  causing  the  difference  tend  to 
compensate  rather  than  add.  In  the  -X  direction  the  opposite  trend  is  seen 
and  -X T is  affected  more  than  -XL. 

Another  interesting  case  to  evaluate  is  the  one  where  X = 0.  In 

this  case 

\ - ± /(C2/4)(T  + At/2)2  - (Ad)2  - (V2/4)(T  + At/2)2  - [(C2  - V2)/C2](V  At/4)2 
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Yt  = ± y(C2/4)(T  + At/2)2  - (Ad)2  - (V2/4)(T  - At/2)2  ♦ [(C2  - V2)/C2]  (V  At/4)2 
For  the  circular  approximation 


Y|_  = ± J{ C2/4)(T  + At/2)2  - (Ad)2 

Y|_  = ± Ac5/4)(T  - At/2)2  - (Ad)2 

The  difference  between  the  corresponding  Y and  Y'  expressions  comes  from  the 
two  additional  terms  under  the  square  root  in  the  equations  for  Y.  As  long 

as  V remains  very  small  in  comparison  with  C,  the  Y values  will  not  be 

materially  smaller  than  the  Y1  values. 

In  a JP  training  assignment,  Lee  Sheldon  numerically  evaluated 
the  error  to  be  expected  in  a variety  of  cases  by  computing  the  X-axis 
intercepts  of  the  ellipses  and  circle  approximations,  with  the  center  of 
the  circles  taken  as  the  origin.  It  is  in  this  dimension  that  the  errors 
are  maximum.  Vehicle  speeds  of  45  and  15  knots  were  assumed.  Results 
are  shown  in  Table  1 when  the  sound  velocity  is  considered  to  be  a constant 
5,000  ft/sec  and  the  ping  duration  is  0.25  sec.  Three  elapsed  times  were 
examined  in  combination  with  two  depth  differentials  between  vehicle  and 
boundary.  Table  2 shows  some  results  obtained  by  using  ray  tracing 
data  in  a refracting  medium. 

On  the  basis  of  these  data,  it  would  appear  that  for  vehicle 
speeds  less  than  say  50  knots  the  circular  approximation  of  the  boundaries 
of  the  annulus  returning  reverberation  at  a given  time  introduces  only 
very  minor  errors.  These  are  certainly  negligible  in  the  light  of  other 
sources  of  error,  such  as  in  the  determination  of  the  scattering  coefficients. 
The  maximum  error  in  X occurs  to  the  rear  of  the  vehicle  where  pattern 
discrimination  generally  is  high. 

Another  study  by  Charles  Williams,  a summer  employee,  showed 
that,  for  vehicle  speeds  up  to  50  knots,  the  difference  between  the  initial 
path  angles  from  the  respective  vehicle  positions  at  transmission  and 
reception  of  reverberation  from  selected  points  on  the  surface  amounted  to 
less  than  1*  for  a wide  variety  of  situations.  The  errors  were  largest 
when  geometries  and  elapsed  times  were  such  that  the  effective  paths 
were  steep.  Since,  in  computing  boundary  reverberation,  the  circular 
approximation  uses  what  is  essentially  a median  of  these  two  angles  for 
both  transmission  and  reception  in  computing  pattern  losses,  the  circular 
assumption  should  provide  an  adequate  model  for  most  practical  cases. 

All  of  the  above  discussion  assumes  that  the  boundaries  and 
vehicle  velocity  vector  are  horizontal.  The  situation  is  much  more 
complicated  when  one  or  both  of  these  conditions  are  not  applicable. 

Work  in  this  area  is  needed. 


A.  B.  P0YNTER 
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APPENDIX  F 

MEMORANDUM  P3502/ABP:dd 

Reg  No.  P3502-136 

From:  Code  P35021  15  April  1968 

To:  Code  P3502 

Subject:  Comparison  of  Computer  Models  for  Obtaining  the  Expected  Level 

of  Boundary  Reverberation  as  a Function  of  Elapsed  Time. 

References:  A.  NAVORD  Conf  Report  5606,  "Analytical  Methods  for  Predicting 
the  Acoustic  Performance  of  Homing  Torpedoes  in  Circular 
Search"  (NOTS  1818)  of  26  July  1957. 


B.  NOTS  P80203  memo,  "Review  of  Mathematical  Model  for 
Computing  Boundary  Reverberation",  of  20  January  1965. 


C.  NOTS  P80203  memo,  "Error  from  Using  10  log10Tas  a measure 
of  Effective  Train  Length  when  Computing  Boundary  Rever- 
beration", of  18  January  1965. 

D.  NUWC  P35021  memo,  "Integration  of  Power  under  a Curve 
Plotted  in  Decibels",  of  10  January  1968. 

E.  NAVORD  Report  4962,  "A  Study  of  the  Effects  of  Refraction 
on  Reverberation  (NOTS  1284)  of  7 November  1955. 

i 

F.  NOTS  P80203  memo  "Doppler  Content  of  Boundary  Reverberation 
Due  to  Vehicle  Translation-Refractive  Environment  Case"  of 
21  December  1965. 

t r 

Over  the  last  few  years  various  aspects  of  the  problem  of  computing 
the  expected  level  of  boundary  reverberation  in  a refractive  medium  have 
been  re-examined.  The  purpose  of  this  memorandum  is  to  compare  the  results 
obtained  from  four  models  selected  as  being  practical  for  computer  applica- 
tion and  to  assess  their  relative  accuracy.  In  general,  the  accuracy  to 
be  anticipated  increases  with  the  complexity.  The  study  is  made  in  terms 
of  surface  reverberation.  In  the  case  of  a flat  bottom,  bottom  reverberation 
is  computed  in  an  identical  fashion. 

f 

! c 

Model  I has  been  used  since  the  inception  of  our  first  ray-tracing 
program  and  was  based  on  Reference  A.  The  form  of  the  equation  used  in  the 
program  is 

r * 

r 

* 

. 

Rs  = S - 2H  + 10  log,0  mf  - Js  + 10  log10T  + 10  log10 
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•R,  is  the  expected  level  of  surface  reverberation  at  a given  time 
•S  is  the  on-axis  source  level  in  dB  re  one  microbar  at  one  yard 
•2H  is  the  two-way  transmission  loss  along  the  ray  path 
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•10  log10ms  is  the  surface  scattering  coefficient  in  dB 
as  a function  of  grazing  angle  of  the  ray  path  at  the  surface 

•Js  is  the  boundary  reverberation  index  of  the  transducer 

•10  log10  t is  the  effective  annulus  width  in  dB 

• 10  log  1Q  allows  for  the  fact  that  horizontal  spreading 

C°s  e0 

is  compensated  for  in  the  outgoing  direction  by  the  fact 
that  the  whole  annulus  back-scatters  toward  the  source.  Here 
X is  the  horizontal  range  traversed  by  the  ray  to  boundary 
intercept,  and  e0  is  the  ray  angle  at  the  source. 

In  this  model  10  log,QT  = where  V is  the  nominal  velocity  of  sound 

(1667  yd. /sec.)  and  At  is  the  ping  duration  in  seconds. 

It  was  shown  in  Reference  B that  the  model  produces  the  best  results 
if  the  elapsed  time  T is  considered  to  be  measured  from  the  midpoint  of 
transmission  so  that,  in  effect,  T = 2t,  the  two-way  travel  time  of  the 
ray  being  run.  This  ray  then  reaches  the  surface  near  the  mid-range  of  the 
insonified  annulus  so  that  the  associated  value  of  transmission  loss,  J,, 
and  grazing  angle  are  reasonable  representative  for  the  annulus  as  a whole. 

Reference  C showed  that  t was  not  a good  measure  of  annulus 
width  where  steep  paths  are  involved.  Model  II  is  a simple  attempt  to 
improve  this  situation.  The  only  change  from  Model  I is  to  compute 
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where  Vsis  the  velocity  at  the  surface  and  es  is  the  grazing  angle  of  the  path 
at  the  surface.  This  is  largely  an  intuitive  "fix",  and  it  is  being  tried 
here  for  the  first  time.  Obviously,  both  models  I and  II  get  into  trouble 
when  T = 2t6c<2t90,  + At/2  because  the  full  ping  is  not  insonifying 
the  surface.  ' In  Model  II  a test  was  made  to  ascertain  whether  or  not  this 
criterion  was  being  met.  If  not,  it  was  assumed  that  the  paths  would  be 
sufficiently  steep  to  approximate  straight  lines.  Since  t would  equal  the 
horizontal  range  to  surface  intercept  of  the  leading  edge, 

T ■ j-  12\*y)  Cose°- 

V is  the  average  velocity  over  these  straight  paths  and  can  be  found  by  dividing 
the  source  depth  by  the  one-way  travel  time  of  the  -90°  ray. 

As  pointed  out  in  Reference  B,  the  above  models  in  some  situations 
could  be  expected  to  give  inaccurate  results  for  long  pings  even  if  the  annulus 
width  were  modeled  adequately.  Values  of  transmission  loss,  vertical  pattern 
losses,  and  grazing  angle  found  by  tracing  a ray  to  any  single  point  in  a 
wide  annulus  cannot  be  expected  to  represent  those  over  the  full  annulus, 
particularly  in  regions  where  one  or  more  of  the  parameters  is  varying 
rapidly.  Moreover,  at  elapsed  times  sufficiently  short  :o  that  all  of  the 
ping  has  not  insonified  the  surface,  the  ray  path  for  which  the  two-way 


F-2 


travel  time  is  computed  no  longer  reaches  the  surface  near  the  mid-point  of 
the  annulus.  These  difficulties  are  overcome  in  the  next  model.  Model  III 
in  this  study  follows  the  scheme  delineated  in  Reference  B and  computes  boun- 
dary reverberation  levels  as  in  Model  I with  10  log  set  to  zero.  This  is 
equivalent  to  a one-yard  annulus  width.  From  the  data  produced  by  the  ray-tracing 
program,  both  the  two-way  travel  time  (2t)  and  computed  reverberation  level 
(R0)  are  plotted  as  functions  of  the  horizontal  range  (x)  to  surface  intercept 
as  determined  for  each  ray  path.  For  each  desired  elapsed  time  (T),  which 
can  be  selected  at  will,  the  limits  of  the  true  annulus  corresponding  to  At 
is  found  by  evaluating  x on  the  time  curve  for  T - At  and  T + jAt  . 

2 2 

The  actual  reverberation  level  at  time  T is  then  found  by  summing  in  random 
phase  the  contributions  of  each  one-yard  increment  between  the  above  limits 
as  shown  by  the  R curve.  The  summation  was  a hand  operation  in  Reference  B 
but  it  was  computerized  in  this  study  in  the  manner  delineated  in  Reference  D. 

These  first  three  models  all  require  a table  input  to  account 
for  the  weighting  effect  of  the  horizontal  patterns  on  the  reverberation 
return.  These  data  are  combined  with  vertical  pattern  losses  and  a 
geometric  correction  for  transducer  pitch  to  give  J,  the  boundary 
reverberation  index  of  the  transducer.  NUWC  computer  program  819001  is 
available  for  generating  the  data  for  the  table  showing  the  average 
weighting  effects  of  the  horizontal  patterns.  For  a non-turning  transducer, 
this  is  a single  value.  The  equations  for  evaluation  J,  were  developed 
in  Reference  E.  Because  of  certain  simplifications  introduced  in  the  concept, 
the  accuracy  of  the  results  is  suspect  when  the  transducer  axis  is  tilted 
considerably  and/or  when  the  steep  sound  paths  are  producing  the  reverberation 
at  time  T.  The  basic  simplifications  are  the  assumptions  that  pattern  effects  can 
be  obtained  when  only  the  patterns  in  the  cardinal  planes  are  known  and  that 
the  variables  can  be  separated  for  integration. 

Model  IV  uses  a different  approach.  The  concept  was  developed 
in  Reference  F for  the  purpose  of  computing  the  boundary  reverberation 
returned  in  doppler  bands  when  the  doppler  is  introduced  by  own  vehicle 
speed.  Briefly,  the  insonified  annulus  is  divided  into  incremental  areas 
bounded  on  two  sides  by  equal  doppler  lines  and  on  the  other  two  sides 
by  equal  travel-time  circles.  Rays  are  traced  to  the  corners  of  these  little 
areas  yielding  the  coordinates  of  the  corners  in  space  so  that  the  areas  of 
each  can  be  computed.  Other  ray  data  permit  evaluation  of  the  transmission 
loss,  scattering  strength,  and  the  transmit  and  receive  pattern  losses  needed 
to  ascertain  the  back-scatter  received  from  each  incremental  area.  When  these 
contributions  are  summed  in  random  phase  (expressed  in  intensities)  for  each 
doppler  band,  one  gets  the  reverberation  levels  as  might  be  observed  by  a system 
which  processes  signals  in  narrow  frequency  bands.  The  summation  of  contributions 
from  all  bands  gives  the  total  surface  reverberation  received.  (The  first 
three  models  yield  only  the  total  reverberation  as  might  be  observed  in 
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,v,  —ox  . Later  this  was  changed  to  35  dB  down  as 
noted  in  Figure  1.  These  data  are  useful  in  interpreting  the  reverberation 
levels  which  will  be  shown  later.  Transducer  pitch  angles  of  0*  and  10*  up 
were  used.  The  same  table  of  scattering  coefficients  versus  grazing  angle 
were  used  in  Models  I through  III.  These  values  were  each  reduced  by  101og102ir 
for  Model  IV  so  that  they  would  reflect  scattering  strength  as  required 
in  that  model.  Computations  were  made  for  all  four  models  using  each  of 
two  quite  different  sound-velocity  profiles  in  order  to  see  that  the 
conclusions  arrived  at  from  comparing  the  results  obtained  from  the  four 
models  were  not  prejudiced  by  some  over-riding  characteristic  of  the  profile. 
The  velocity  profiles  are  shown  in  Figure  2.  A transmit  frequency  of 
20  kHz  was  chosen,  and  an  appropriate  table  giving  the  attenuation  losses 
as  a function  of  depth  was  input  with  each  profile.  Ping  durations  of  40  and 
250  milliseconds  were  considered.  Model  IV  assumed  a horizontally-directed, 
straight-running  platform  at  a 40-kt.  speed,  and  reverberation  was  computed 
In  1/4-kt.  doppler  bands.  The  vehicle  speed  Is  an  artlfax  here  since  we 
are  Interested  only  In  the  total  return.  It  was  selected  to  Insure  suf- 
ficiently small  area  Increments  to  yield  good  accuracy.  Actually,  the 
results  can  apoly  to  a stationary  system  since  speed  Is  not  considered 
in  the  other  three  models.  A 120-dB  source  level  was  used  In  all  cases. 

The  transducer  Is  assumed  to  be  at  a 1000-ft.  depth. 

It  proved  to  be  Impracticable  to  show  the  results  for  all'  four 
models  on  a single  graph  because  of  the  clutter.  Consequently,  for  each 
combination  of  environment,  ping  length,  and  pitch  angle,  two  figures  will 
be  used.  The  first  will  compare  the  data  from  Models  I,  II,  and  III.  On 
the  second  graph  the  data  from  Model  III  will  be  repeated  and  compared  with 
that  from  Model  IV.  This  Is  a natural  grouping  for  our  purpose  In  that 
the  first  three  models  assume  comnon  values  for  J,.  Then  by  comparing 
the  most  accurate  of  these  models  with  the  doppler-band  method,  one  should 
be  able  to  obtain  an  Idea  of  the  accuracy  with  which  we  new  evaluate  J.. 
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Figure  3 shows  the  surface  reverberation  levels  computed  in 
Environment  A for  a 40-ms  ping  and  a 0°  pitch  angle.  With  the  source  at 
a 1000-ft.  depth  the  middle  of  the  ping  has  a two-way  travel  time  to  the 
surface  of  about  0.405  seconds  via  a vertical  path  (e0  = -90*).  The 
velocity  profile  is  dominated  by  negative  gradients  and  the  -15*  ray 
is  very  nearly  the  last  to  reach  the  surface.  The  3.8-sec.  ping  interval 
was  selected  with  this  in  mind.  The  combined  vertical  patterns  largely 
determine  the  peaks  and  valleys  in  the  respective  reverberation  curves. 

Of  course,  the  trailing  off  of  the  reverberation  at  times  greater  than 
about  1.7  seconds  is  caused  by  the  transmission  loss  increasing  more  rapidly 
than  the  pattern  losses  are  decreasing.  In  addition,  the  scattering  coefficient 
tends  to  be  smaller  at  the  lower  grazing  angles.  All  three  models  tend  to  be 
in  good  agreement  beyond  about  1.05  seconds  when  the  major  lobes  of  the 
vertical  patterns  govern.  The  integrated  method  (Model  III) tends  to  smooth  out 
the  narrow  peaks  and  valleys  generated  by  the  other  two  models.  It  is 
undoubtedly  the  more  accurate  since  it  does  not  consider  the  particular 
values  of  pattern  and  transmission  loss  associated  with  a data  point  as  being 
necessarily  representative  of  the  entire  annulus  returning  reverberation  at 
that  elapsed  time.  The  cosine  correction  to  the  annulus  width  used  in 
Model  II  seems  to  yield  good  agreement  with  the  integrated  model  over  the 
broad  minor  lobe.  (The  surface  velocity  used  in  Model  II  is  not  signifi- 
cantly different  than  the  nominal  velocity  used  in  the  first  model).  However, 
it  tends  to  drastically  over-compensate  on  the  other  minor  lobe  as  can  be 
expected;  the  cosine  is  rapidly  approaching  zero  as  the  path  angle  approaches 
-90°. 

The  results  for  the  250-ms  ping  under  otherwise  similar  conditions 
are  presented  in  Figure  4.  The  curves  for  Models  I and  II  have  the  same 
shape  as  in  Figure  3;  they  are  merely  at  a higher  level  because  of  the 
longer  ping.  The  agreement  between  the  three  models  is  still  good  beyond 
about  1.2  seconds,  but  the  importance  of  the  integrated  method  when  long 
pings  are  used  is  quite  evident  at  shorter  elapsed  times.  The  results  of 
its  use  are  most  striking  in  modeling  the  onset  of  boundary  reverberation. 

The  flat  portion  of  Model  III  curve  centered  at  about  T = 0.35  is  due  to  the 
narrow  lobe  being  fully  covered  by  some  portion  of  the  effective  train 
length  while  the  rest  is  contributing  negligible  return.  With  the  quarter- 
second  ping  the  broad  minor  lobe  dominates  the  return  until  such  time  that 
the  contribution  from  the  major  lobe  begins  to  build  up. 

Figure  5 shows  the  corresponding  data  for  the  40  ms-ping  in 
Environment  A when  the  transducer  is  pitched  up  10*.  The  gross  effect  is 
to  make  the  peaks  and  valleys  In  the  expected  level  of  reverberation  more 
narrow  than  was  the  case  when  the  transducer  was  directed  horizontally.  This 
is  to  be  expected  since  each  path  making  a particular  angle  with  the  transducer 
axis  is  10*  steeper  in  the  oceanographic  coordinate  system  and  returns 
reverberation  at  a shorter  elapsed  time.  The  difference  In  results  for 
the  three  models  are  of  the  same  order  as  were  found  for  the  horizontally- 
directed  transducer. 
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Figure  6 shows  the  results  under  the  same  conditions  as  in  the 
preceding  paragraph  except  that  the  ping  duration  is  increased  to  250  ms.  As 
anticipated,  the  integration  method  (Model  III)  gives  a decidedly  different 
curve  than  the  other  two  models  at  times  less  than  one  second.  This  model  permits 
development  of  the  onset  of  the  surface  reverberation  starting  where 

1 - *w  - 

The  flat  portion  starting  at  0.3  seconds  covers  the  time  region  where  the 
narrow  minor  lobe  alone  is  contributing  a substantial  return.  The  level 
then  rises  as  the  next  lobe  also  contributes.  The  drop-off  starting 
at  about  0.52  seconds  is  caused  by  the  narrow  minor  lobe  dropping  out  of 
the  picture.  The  next  flat  portion  occurs  when  only  the  second  lobe  Is 
contributing  materially.  The  shape  of  the  remaining  portion  of  the 
curve  is  rather  obvious. 


The  next  group  of  figures  (Figure  7-10,  inclusive)  compare  the 
results  from  the  same  three  models  when  Environment  B applies.  Note  that 
the  dB  scale  has  been  changed.  The  two-way  travel  time  from  1000  feet  to 
the  surface  via  a vertical  path  is  about  0.418  seconds.  Surface  rever- 
beration is  continuous  thereafter  to  an  elapsed  time  greater  than  the  ping 
interval  which  was  taken  as  10  seconds.  The  path  yielding  this  two-way 
travel  time  to  the  surface  has  an  initial  path  angle  of  approximately 
+3.7°  with  respect  to  the  horizontal.  At  the  scale  plotted  there  is  no 
discernible  difference  in  the  results  from  Models  I,  II,  and  III  at  elapsed 
times  beyond  1.6  seconds.  Therefore,  only  the  first  portion  of  the  ping 
interval  is  shown  in  these  figures. 

For  Environment  B the  coincidence  of  the  three  curves  begins 
shortly  after  the  major  lobe  of  the  patterns  assume  dominance.  This 
occurs  at  an  earlier  time  than  in  Environment  A for  a given  pitch  angle 
since  refraction  is  such  that  a ray  with  a given  initial  angle  will  tend 
to  reach  the  surface  much  sooner.  Also,  in  the  case  of  Environment  B, 
the  reverberation  declines  at  a lower  rate  after  the  peak  is  reached  so 
that  the  integrated  curve  does  not  begin  to  rise  above  the  levels  shown 
for  the  other  two  models  later  in  the  ping  cycle  as  occurred  to  some  degree 
in  the  case  of  the  first  environment. 


At  the  shorter  elapsed  times,  the  differences  in  the  three  curves 
tend  to  be  quite  similar  (for  each  of  the  four  combinations  of  ping  length 
and  pitch  angle)  to  those  prevailing  for  the  same  combination  when  Environ- 
ment A applied.  Model  III,  the  integration  method, tends  to  smooth  out  the 
peaks  and  valleys  as  shown  for  the  other  two  models,  particularly  when 
they  are  narrow  with  respect  to  the  ping  duration.  This  method  also  is 
better  in  modeling  the  onset  of  boundary  reverberation.  For  the  shorter  ping, 
the  annulus  width  correction  applied  In  Model  II  appears  to  compensate 
adequately  in  the  region  of  the  broader  minor  lobe,  but  it  tends  to  over- 
compensate when  the  paths  become  steep  enough  for  Cos  0O  to  become  very 
small.  The  introduction  of  a 10°  up  pitch  produces  about  the  same  effect 
in  both  environments. 
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Now  we  turn  to  the  interesting  task  of  comparing  the  results  of 
the  integrated  method  with  those  of  the  doppler-band  model.  Figure  11 
shows  the  data  obtained  in  Environment  A for  a zero  pitch  angle  and  a 
pulse  length  of  40  ms.  The  larger  scale  again  is  used  for  clarity. 

Earlier,  the  hypothesis  was  advanced  that  agreement  between  data  from 
the  two  models  at  the  longer  elapsed  times  (where  flat  paths  obtain)  would 
be  a good  indication  that  the  doppler-band  program  (Model  IV)  was  func- 
tioning as  intended.  The  agreement  is  reasonably  close  over  the  region 
dominated  by  the  major  lobes  of  the  patterns,  but  there  remains  a slight 
difference  at  3.8  seconds.  However,  in  this  environment  the  path  angles 
yielding  this  elapsed  time  is  still  over  15°  off  the  transducer  axis.  As 
anticipated,  some  rather  substantial  differences  are  observed  in  the  region 
dominated  by  the  minor  lobes  of  the  patterns.  Since  there  is  little  reason 
to  suspect  that  Model  IV  is  materially  less  accurate  in  this  region  than 
in  any  other,  the  evidence  tends  to  substantiate  fears  that  the  current 
method  for  evaluating  J,  is  in  error  when  large  off-axis  path-angles  are 
involved.  In  this  example  (as  in  subsequent  ones  where  the  transducer  is 
oriented  horizontally)  the  differences  in  reverberation  levels  at  the  peaks 
are  rather  moderate  (no  more  than  2 dB).  The  wide  troughs  are  something 
else  again.  The  higher  minimums  observed  in  Model  IV  results  are  believed 
due  to  the  complexity  of  the  minor  lobe  structure.  For  the  results  shown 
in  Figure  11,  all  of  the  incremental  areas  contributing  a return  at  an 
elapsed  time  of  about  one  second  are  not  subjected  to  maximum  pattern  loss. 

Of  course,  in  practical  applications  the  actual  values  in  the  valleys  are 
not  likely  to  be  important  since  the  interference  level  probably  will  be 
dominated  by  volume  reverberation  or  noise  at  corresponding  times. 

Figure  12  shows  comparable  data  when  the  ping  duration  is  increased 
to  250  ms.  The  reverberation  levels  agree  to  about  the  same  degree  as  for 
the  shorter  pulse  over  the  peaks,  and  the  differences  in  the  valleys  have 
been  reduced.  The  narrow  peaks  and  the  valleys  have  been  smoothed  consider- 
ably. 

Figures  13  and  14  compare  the  results  from  Models  III  and  IV  for 
ping  lengths  of  40  and  250  ms,  respectively,  when  the  transducer  is  pitched 
10*  up.  Environment  A still  applies.  The  extent  of  the  agreement  between 
the  respective  curves  is  substantially  the  same  as  for  the  comparable 
cases  when  zero  pitch  was  used  except  at  the  shortest  times  shown. 

It  is  felt  that  when  path  angles  are  involved  which  approach-90*  the 
geometric  correction  factor, 

-10  log )0  [Cos  (e0  - s)/Cos  e„] 

(where  5 is  the  transducer  pitch  angle) 

in  J,  over-compensates  and  makes  the  levels  a few  dB  higher  than  they 
should  be  for  Model  III.  At  the  longest  times  (where  the  initial  path 
angles  are  approaching  -15“  ) the  angles  with  respect  to  the  transducer 
axis  are  approaching  -5*.  One  might  expect  the  two  curves  to  coincide. 

The  one  for  the  integrated  method  is  still  a few  tenths  of  a dB  above 
the  other  one  as  it  was  for  the  zero  pitch  cases.  This  leads  one  to 
suspect  that  the  Cos  (0O  - C)  / Cos  e0  correction  is  slightly  over- 
compensating  even  at  these  angles. 
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Figures  15  through  18  compare  data  from  Model  III  and  IV  for  . j 

the  various  combinations  of  pulse  lengths  and  pitch  angles  when  Environment 
B applies.  In  all  four  figures  the  two  curves  come  into  coincidence  near 
or  before  an  elapsed  time  of  about  1.9  seconds  which,  in  Environment  B, 
corresponds  to  an  ititial  path  angle  of  approximately  -12  degrees.  For 
zero  pitch.  Figure  15  and  16,  the  agreement  persists  to  the  10-second  ping 
interval  which  corresponds  to  a path  angle  of  approximately  +4  degrees. 

This  is  the  most  convincing  evidence  that  Model  IV  is  properly  programmed 
since  it  is  under  these  conditions  that  J,  should  be  most  accurate. 


For  pitch  angles  of  10  degrees  up.  Figures  17  and  18,  the  near 
perfect  coincidence  begins  near  1.7-sec.  elapsed  time,  corresponding  to 
a path  angle  of  about  -13.5  degrees.  Note,  however,  that  this  corresponds 
to  -3.5  degrees  from  the  transducer  axis.  On  the  basis  of  path  angle  alone 
one  might  think  that  coincidence  should  occur  earlier.  It  is  suspected 
that  the  delay  is  caused  by  the  Cos  (90  - £) /Cos  90  term  in  J,  over-com- 
pensating for  the  pitch.  The  agreement  persists  to  about  an  elapsed  time 
of  7 seconds,  at  which  time  the  reverberation  level  computed  by  means 
of  the  integration  method  begins  to  fall  slightly  below  that  computed  by 
Model  IV.  At  7 seconds  the  initial  path  angle  is  about  +0.6  degrees  while 
the  angle  with  respect  to  the  transducer  axis  is  approximately  +10.6  degrees. 
It  seems  possible  that  the  slight  divergence  from  7 to  10  seconds  results 
from  the  Cos  (90  - S)/Cos  e0  correction  to  J,  tending  to  under-compensate 
at  appreciable  off-axis  angles  when  the  sign  of  the  pitch  angle  is  opposite 
that  of  the  path  angle.  However,  the  evidence  is  by  no  means  conclusive  since 
the  divergence  does  not  increase  consistently  as  the  elapsed  time  approaches 
10  seconds.  The  differences  are  so  slight  that  they  well  may  be  due  to  errors 
in  the  numerical  integration  process  in  one  or  both  models.  Although  data 
points  are  taken  at  relatively  short  intervals,  they  are  still  at  finite 
intervals  apart,  and  linear  interpolation  is  used. 
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For  the  relatively  short  elapsed  times  (less  than  1.5  sec.),  the 
differences  in  results  from  the  two  models  in  Environment  B are  essentially 
the  same  as  they  were  when  Environment  A was  used  with  one  exception.  For 
both  ping  lengths  with  zero  pitch,  the  minor  lobe  at  minimum  time  for  the 
integrated  method  (Model  III)  is  a little  over  one  dB  lower  than  for  the 
doppler  method  for  Environment  B computations.  For  Environment  A the  dif- 
ference is  about  the  same  amount  but  in  the  opposite  direction.  This  is 
attributed  to  the  fact  that  we  changed  the  constant  pattern  value  assumed 
for  the  back  portion  of  the  pattern  as  was  announced  earlier.  The  transition 
was  most  severe  when  Environment  A was  used,  and  these  results  are  considered 
to  be  unreliable.  For  the  10  degree  up  pitch,  this  portion  of  the  vertical 
pattern  does  not  come  into  play  In  Model  III. 

It  was  decided  to  make  one  further  effort  to  clarify  the  effect 
of  transducer  pitch  by  computing  a case  involving  a 40  degree  up  attitude. 

As  other  conditions  it  was  decided  that  Environment  B and  a 40-ms  ping 

would  combine  to  produce  the  most  useful  results.  Remaining  parameters 
have  the  same  values  used  throughout  the  study.  Figure  19  compares  the 
results  produced  by  means  of  Models  III  and  IV.  The  on-axis  ray  (eQ  * -40°) 

reaches  the  surface  at  T = .647593  sec.  At  this  pitch  angle,  at  least 

in  Model  III,  the  major  lobe  of  the  pattern  controls  the  reverberation 
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level  over  a relatively  small  portion  of  the  total  ping  interval.  The 
first  minor  lobe  above  the  transducer  axis  controls  the  level  of  the 
peak  following  the  onset  of  surface  reverberation.  The  first  minor  lobe 
below  the  axis  governs  the  reverberation  level  over  the  last  75%  of  the 
10-sec.  ping  interval.  It  is  difficult  to  visualize  physically  just  how 
the  patterns  interact  in  Model  IV. 


In  comparing  the  results  from  the  two  models,  it  is  noted  first 
that  the  two  curves  coincide  only  briefly  at  a time  centered  about  T = 0.85  sec. 
This  corresponds  to  an  initial  path  angle  in  the  vertical  plane  of  -29.1  degrees 
or  +10.9  degrees  with  respect  to  the  transducer  axis.  At  times  shorter  than 
this  the  curves  diverge  with  the  integrated  method  giving  the  higher  values. 

At  about  0.43  sec.  the  separation  is  nearly  8 dB.  At  times  greater  than 
0.85  seconds,  the  curves  cross  and  there  again  is  an  increasing  difference 
but  with  the  integrated  method  yielding  the  lower  values.  The  large  difference 
in  the  vicinity  of  1.9  seconds  can  be  attributed  to  the  fact  that  the  null 
in  the  vertical  plane  patterns  is  not  characteristic  of  the  whole  annulus 
which  returns  reverberation  at  times  of  this  order.  If  this  region  is 
discounted,  then  the  remainder  of  the  ping  cycle  shows  a very  gradual 
increase  in  the  difference  between  the  two  curves  from  about  1.7  dB  at  2.5  sec. 
to  2.8  dB  at  10  sec. 


The  above  data  would  seem  to  support  certain  conclusions.  In 
the  first  place,  the  doppler-band  method  (Model  IV)  is  working  properly. 
Although  it  is  judged  to  be  the  most  accurate  way  of  computing  the  expected 
level  of  boundary  reverberation,  it  increases  the  computing  time  by  a factor 
of  ten  over  that  required  by  any  of  the  other  methods.  Therefore,  its  general 
use  as  part  of  the  ray  tracing  program  is  not  recommended.  It  should  be  a 
separate  program  reserved  for  special  cases.  Even  then  it  requires  insertion 
of  appropriate  equations  for  the  transducer  patterns  for  obtaining  pattern 
losses  in  any  directions.  Alternately,  one  could  use  a matrix  of  measured 
patterns  in  such  a large  number  of  planes  that  accurate  results  could  be 
obtained  by  interpolation. 

The  other  three  models  have  in  common  the  inaccuracy  resulting 
from  errors  in  computing  values  of  the  boundary  reverberation  index,  for 
various  geometries  which  obtain,  as  functions  of  transducer  depth  and  elapsed 
time  in  the  ping  cycle.  The  evidence  is  that  the  errors  inherent  in  the 
present  method  of  computing  J,  are  small  for  path  angles  falling  in  the 
major  lobe  of  the  transducer  when  the  pitch  angle  is  not  much  greater  than 
10  degrees  from  the  horizontal.  For  small  pitch  angles,  the  error  in  the 
first  minor  lobe  may  be  tolerable  in  veiw  of  the  usual  uncertainty  as  to 
the  proper  values  to  assign  for  the  scattering  coefficient  per  unit  area. 

Paths  falling  in  the  nulls  between  pattern  lobes  lead  to  values  for  J,  which 
are  too  nigh,  but  in  general  the  boundary  reverberation  at  corresponding  times 
will  be  below  other  types  of  interference  and  is  of  no  practical  consequence. 
Model  I,  the  method  presently  incorporated  in  our  ray  tracing  program,  also 
suffers  from  underestimating  the  width  of  the  insonified  annulus  as  the 
paths  become  steeper.  Model  II  will  improve  this  situation  over  a con- 
siderable range  of  path  angles,  but  this  improvement  is  unimportant  relative 
to  another  source  of  error  common  to  both  models  I and  II.  Here  we  refer 
to  the  assumption  that  values  of  J,  and  transmission  loss  computed  to  a 
single  point  in  the  annulus  are  characteristic  of  the  entire  insonified  area. 
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Except  for  very  short  pings,  this  assumption  can  introduce  substantial 
errors  when  the  rate  of  change  of  one  or  more  of  these  parameters  is  large. 

On  this  basis  it  seems  logical  to  adopt  the  integration  method  (Model  III) 
as  the  regular  routine  in  the  ray  tracing  programs.  Surprisingly  enough, 
this  can  be  accomplished  without  a significant  increase  in  run  time  on  the 
computer.  With  a suitable  selection  of  rays,  accurate  modeling  of  the  expected 
level  of  boundary  reverberation  should  be  achieved  within  the  accuracy 
inherent  in  J, . If  an  improved  method  of  evaluating  this  factor  is  found, 
the  benefits  will  automatically  accrue  in  the  reverberation  computation. 


Because  of  the  other  variables  involved,  it  is  difficult  to  pin 
down  the  actual  errors  in  J,  by  comparing  the  reverberation  levels  computed 
by  means  of  Models  III  and  IV.  During  the  course  of  this  study,  a technique 
was  envisioned  whereby  J,  could  be  investigated  directly. 
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Following  are  listings  for  the  OOP  program  as  used  for  the  test  probleai 
whose  output  Is  Illustrated  In  Figures  3-1  to  3-3.  There  Is  no  subroutine, 
SPRCMP  and  the  sample  OXL,  RRF,  and  TVGF  were  created  for  Illustrative  pur- 
poses only. 

Patterns  from  the  listed  OXL  are  Illustrated  In  Figures  G-l  thru  6-3. 
They  are  all  radially  symmetrical  about  their  axes. 
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PROCEDURES  FOR  OOF 


p»uc 

COMMON  /(HaNk 


COMMON  /CPANP 


lOMMCf!  / CC'JUk 


COMMON  /CFCNS 


COMMON  /CICNS 
COMMON  /CfCNS 
COMMON  /CMCNS 
COMMON  /CMCNf. 
COMMON  /CIMPA 
COMMON  /CINDA 

COMMON  f C I N 0 f 
COMMON  /CINDt 
COMMON  / C I N J f 
COMMON  /CJNPE 

COMMON  / C T N P U 
COMMON  /IINPU 
COMMON  /CINPU 
COMMON  /C1MPU 
COMMON  /C!MPU 


DCOMMM  1 

/ lMBAND,lMBNDlvNBAND,MBAND.MNDMD,rNDAMD,VFTTRM  DCOMMM  2 

900  , 801  DCOMMM  3 

DCOMMM  4 

/ O*Nt>C801>,RRFS(800>#FGAM<800>.8MDOUT<80l>  DCOMMM  S 

< l MB  NO  1 ) , (l MB AND > , (LMBAND > » < iMBMD 1 > DCOMMM  6 

DCOMMM  7 

U LMRBvLMrRStiMKS9(.ffCS2««SVfff  vMBTRvffT«CVT  DCOMMM  0 

961?,  17  , 400  » 800  DCOMMM  9 

DCOMMM 10 

T/  TCS  <6),fPT5,I«S?C4>,M,T«S3(4>,f2,FA,Fl0.F?0  DC0MMM11 

< IMTtS  > ICDMIN1? 

DCOMMM 13 

1/  FL0610VINFNTVF1?3VPI  «Tt«OFl  ,DECRAD  tSMlFT,V»KT»MC3  DCOMMM  14 


1/  LObApi ,F1mim,F3#F90,F180 

T/  HBANDC3l9MUMIT(?l9MSFRD<S>tMSFftM<5>*NFBOMC2> 
T/  HEADS«4l1*NT0Tl(4>tMT0T2<3>9NT0T3(3>»MC8«f I 
T/  ONf  6 A ,DELT  ,CSK  SI ,SMKSI «KS1D ,DR  ,FC OVS ,C OKI 
M EIPS ,F2SO,F7?vFC03»8W1NT,L04MV1 


DCOMMM 15 
DCOMMM 14 
0C0MMM17 
DCOMMM 18 
DC0MMM19 
DCOMMM20 
DC0MMM21 
DCOMMM22 


COMMON  /CPR1N 
COMMON  /CSPRE 
COMMON  / C T A P t 
COMMON  /CT81H. 
CORRON  /CXCNS 
0 IRE  US  ION 

0 IK!  NS  I ON 
D IRf NS I ON 
t IRFNSION 

EOUIVAL f Nit 
EQUIVALENCE 
equivalence 
EQUIVALENCE 
EQUIVALENCE 

1 QUIVALENCF. 

I Q1II  VALENCE 
EQUI  V Al ENCE 


u nnarcs,nanoatuq,3).narcnt<24>  BC0RRN21 

1/  l EN  I E»  ,C0,  ENO.E  HTEP  .XTSBNO  .NOBTTR  , NOPINT  BCONNN22 

>/  NOSURF , NOTAPE. NOVOLR, PLOT, SPREAD, TIN CNP, TOTALS .TVCBCONNN23 
r / RELBND  DC0NNN24 

DC9NNN2S 

1/  IIIC,1>ATEC2>,!(VI»0  DCONNN20 

1/  C f . 1LPHC ,PlR4,tRTTR,l04RV ,S ,«S 1 DCONNN27 

1/  VS,RH1DTH,DELT2,FIR0,NBEAN .ORE  4AD.TMTNAB  DC0NNR28 

T/  PULSE .1PEVRV  DC0NNN29 

W TIRE (A00I.SPBEBC130.33  DCONNN30 

ILNTINI  CLNSPRD 1 DC0NNN31 

BCORRB32 

1/  KIR ir.rTRAI ,PA6E ,NPA4F .NPSTBT  BCORRB  33 

0/  NSPPl.NSPBH.NSPBMl  DC0NNN34 

/ AB1.BB1.1NPT.IPBT.IRLT  DC0NNN3S 

PC  PEL IND.OELOEP.f ACTOR, ITABLE  DE0NNNJ4 

1/  kO,a1,k2,K3,BS,K6,k8.k10.KAO  BC0RRN37 

NP0UNBC2I  DC0NNN38 

N0UTPT1141  DCONNN39 

1PATA11)  DCONNN40 

LCLACSC16)  DC0NNN41 

BC0RRN42 

IKC.FO.QECVl.lKl.RRlT), CBS  .INN  I)  BC0RPN43 

(RSURT .NBOUNBI  BC0RRN44 

< IBt  ANK.HFAOSlIt) 1 , (KUN1T.H0UTPT1  DC0NNN43 

< tFlRO,ElBO),<lBC.lBATA>  DC0NNN46 

ILRTIR.NARQAT121 ,3)1 .INTINE.NANCNTI211)  DC0NNN47 

tlNSPRD,NANDATt22,3)),<LF  LASS,  CENTER)  BC0NNN48 

INSSPBO.NARCNT ( ? 2 > > , (NB SPBB .N ANCN T < 2 3> > DC0NNN49 

<NVSPBB,NARCNT<24)>  DCONNN50 

BC0NNN31 

INERT, KSI.KSID .L04RV ,L  044P 1 .L04RV I BCONNN52 

API ,BR 1 .PAGE  BCORNN33 

CENTER.60.ENB, FILTER, KTSBNB .NOBTTR .NOPBNT  BCONRN34 

NO  SUB  E .NOTAPE .NOVOLR .PLOT .SPBE AB.T IN CNP. TOTALS .TVODCONNN55 
PFLPNO.VPTTPN.LELABS  BC0NNN56 

BCONNN37 

BC0NNN38 

»VS(3,e01) ,PVB 15,801 1.BVV <3,80 11.BVT 13,801 > BC0NNN3* 

1 1LRNT.LRBN01)  > BCONNNtO 


E QUI VALENCE 


PEVFBB19AT2),RV(3, 801,31 

ILNPB),(LRNT,LNBNB1,3> 

IPE VERB  1 11.RVS1 1,1 >,RV< 1,1, 111 


ARINtAOOl.TNINlAOO) 


BC0NNN61 

BC0NNN62 

BC0NNN63 

9C0NNN64 

BCONNNAS 

BCONNNtO 

DC0NNN67 


68  CORRON  RP« 1400) .PBT1A00) ,BBTH I4C01 ,BBHf 4001  BC0NNN68 
1.9  C I LRAS  A RAX  ■ NO.  OF  BATS  SORTED  BV  'BATS*)')  BE0NNN69 
70  END  BC0NNN70 


E 
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7 

7; 

?: 

7t 


DC0MN4  PROC 
COMMON 
C 

END 

DCOMNS  PROC 
COMMON 
COMMON 
COMMON 
C 

c 

COMMON 

DIMENSION 

C 

C 

EQUIVALENCE 

EQUIVALENCE 

C 

C 

END 

DCOMN6  PROC 
DATA 
DATA 
C 
C 
C 

c 

c 

END 


rCSCAM<801> 

(LM8ND1) 


X (fi '■'0),  T (80P)  fTHA  (800),  IMP  (8  00  >.  A <$00>,DOP  18  CO) 

CO$TMA<$00),COSTMB<8CO> 

OMT(800),COSONTC800>vSlNCMT(8P0) 

( LMKS2  ) 

RBIAC400) ,RBTA(400) ,RBTMA(400)  ,RBHA<400) 
RBIB(400>9RBTB(400)9RBTMB<4C0> VRBMB(40C) 

( (LMKS)  > 

(RBIBfDOP) 

(RBTB,T(400))t(RBTHBvTMB(400))»<RBM6,R(400)) 

( (LMKS)  ) 


LMBAND /800/, LMBND 1/801/ 9LMK$/4CP/ ,LPkS2/80 0/ 
LMRfi/9612/fLMTlM/400/fLMTRS/17/,LMSPRD/150/ 

LMRB  * 4 • LMNT  * LMBND 1 
LMBNDl  * LMBAND  ♦ 1 
LMKS2  * LMKS  • 2 
LMNT  * 3 


OCOMPN  71 
DCOPPN72 
DCOMMN73 
DC0PPN74 
OCOPMN  75 
DCOMMN76 
OCOMPN77 
DCOMMN78 
DC0MMN79 
DC0PMN80 
DCOMMN81 
DCOMMN8? 
OCOMMN83 
DC0MMN84 
DCOMMN85 
DC0MMN86 
DC0MMN87 
OCOMMN88 
DC0MMN89 
DCOMMN90 
DCOMMN91 
DCOPMN92 
DCOPPN93 
DC0MMN94 
DCOMMN95 
DC0MMN96 
DC0MMN97 
0C0MMN98 
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I N » E > 


rati 


PROCEDURES  EOR  OOP 


t , 

♦ 


SYMBOL 

* B * * 1 

* * * « 

ALBMC 

• 

27CO 

AM 

- 

35CO 

531N 

BAND 

- 

SCO 

BNDOUT 

- 

SCO 

• Ml 

• 

35CO 

53  IN 

BUI  DIM 

- 

28CO 

BU1NT 

- 

19C0 

CO 

- 

27CO 

COM 

- 

18CO 

CB  AND 

- 

2CL 

5CL 

C COUNT 

- 

8CL 

CENTEA 

- 

22C0 

A8ED 

SOLS 

CfCNST 

- 

1 1 CL 

net 

ISC  L 

CNCNST 

- 

16CL 

17CL 

CINDAT 

- 

18CL 

19CL 

CINDEE 

• 

21CL 

22CL 

23c  L 

CINMUT 

- 

26CL 

27Ct 

20C  L 

COSOMT 

- 

78C0 

COSTNA 

- 

77C0 

C05TNB 

- 

77C0 

CBRINT 

- 

33CL 

CSKS1 

- 

18C0 

CSPAED 

- 

30CL 

CTABE 

- 

35CL 

CTBLKF 

- 

36CL 

CXCNST 

- 

37CI 

DO 

- 

26C0 

DB1TH 

- 

27CO 

DCOMN 1 

- 

1 

DCOMM2 

- 

58 

DCOMN3 

- 

66 

DC0NN4 

- 

71 

DCONN5 

- 

75 

DCOMN6 

- 

90 

DE4RAB 

- 

14C0 

DEEDED 

- 

36C0 

DE1IND 

- 

36C0 

BEIT 

- 

18C0 

DELT2 

- 

28C0 

BOB 

- 

76C0 

8SEB 

DD 

- 

18C0 

END 

- 

22CO 

SOLS 

EXPS 

- 

19C0 

FO 

- 

43E0 

FI 

- 

11C0 

FIO 

• 

11C0 

F180 

• 

15C0 

F1E3 

- 

14C0 

E ININ 

- 

1SC0 

P Z 

- 

ttco 

F20 

- 

11CO 

F3 

- 

1 SCO 

f 4 

- 

11CO 

F*0 

- 

15CO 

FACTO* 

- 

36C0 

FC03 

- 

t9CO 

FCOVS 

- 

18CO 

rCSCAM 

- 

72CO 

ECAN 

- 

SCO 

MITE* 

- 

22CO 

3414 

FLOS 10 

- 

14CO 

f NBAND 

- 

2CO 

FPT5 

- 

11CO 

FEZ 

• 

19CO 

FIDO 

- 

28CO 

46EB 

FZSS 

- 

19C0 

BO 

- 

22CO 

5416 

NBAND 

- 

16CO 

REFERENCES 


ZOCL 
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1 k I E I 


>UI 


PA8E  4 


HEAPS  - 17CO  45EP 

HEP  - 17CO 

HE POP  - 16CO 

HOUTPT  - 39P1  45EP 

HSPPP  - 1600 

HSPRM  - 1600 

HT0T1  - 1700 

HT0T2  - 1700 

HTOT  3 - 1700 

HUH1T  - 1600  45EP 

IBLAHA  - ASEP 

1PATA  - A OP  I A6EP 

1 PATE  - 26(0 

1PC  - 26(0  AtEP 

109  - 26(0 

IEZP0  - A6EP 

1HENT  - 14(0  52RL 

IHPT  - 3SC0 

IPE9RT  - 29(0 

1PLT  - 35CO 

1PPT  - 35(0 

I TABLE  - 36C0 

AO  - 37C0  ASEP 

A 1 - 37(0  ASEP 

A 10  - 37(0 

A2  - 37(0 

A3  - 37CO  ASEP 

AAO  - 37(0 

A5  - 37(0 

A6  - 37(0 

A8  - 37(0 

A$1  - 27C0  52AL 

ASIP  - 18(0  52PL 

AT  - 8(0 

ATS8NP  - 22(0  54L8 

ATT  - 8C0 

LELA8S  - A1P1  A8EP  5618 

LMANO  - 2(0  91  PA 

LMN01  - 2(0  9 »9A 

LHAS  • 8(0  91PA 

LHAS2  - 8(0  91PA 

LHP T - A3EP 

LPBB  - SCO  92PA 

LHSPAP  - A8EP  92PA 

LATIN  - A7EP  92PA 

LNTRS  - 8(0  92PA 

L06API  - 1 SCO  52 AL 

L08H9  - 27(0  52AL 

L08H9I  - 19(0  S2AL 

HBANP  - 2(0 

HNBNO  - 2(0 

NAHCNT  - 21(0  47E0  A9EP 

HAHPAT  - 21C0  47E0  48EP 

NBAHP  - 2(0 

N8EAH  - 28(0 

NOOUNP  - 38P1  AAEP 

NBSPRP  - A9EP 

HATH  - 8C0 

NNAHES  - 2 ICO 

NOBTTN  - 22CO  5418 

NOPRNT  - 22CO  5418 

NOSURF  - 23(0  55L8 

NOTAPE  - 23C0  55L8 

NOVOLH  - 23CO  5518 

NPA8I  - 33CO 

NPSTftT  - 33CO 

NSPR1  - SACO 

NSPRH  - SACO 

NSPRNl  - SACO 

NtSPRP  - A9EP 
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1 N 0 E 

1 

NSURf 

SCO 

44E0 

NTIPE 

- 

47EO 

NTPA1 

- 

33CO 

ntpin 

• 

33CO 

NVSPRI 

• 

SOEt 

OPERA 

- 

mo 

OPERA A 

- 

28CO 

OPT 

- 

78C0 

PACE 

- 

33CO 

S31N 

PI 

- 

14C0 

PI23 

- 

14C0 

PINS 

- 

27CO 

PLOT 

- 

23CO 

SSLS 

PULSE 

- 

29CO 

« 

- 

76CO 

86EO 

ROM 

- 

68CO 

RRNA 

- 

RICO 

ROHR 

- 

02D1 

86E0 

ROT 

• 

68C0 

RRTA 

- 

RICO 

RRTR 

- 

8201 

86E0 

RRTN 

- 

68C0 

RRTPA 

- 

RICO 

RRTPR 

- 

8201 

ROES 

RRI 

- 

68C0 

RRIA 

- 

RICO 

RRIR 

- 

8201 

8SE0 

RECV 

- 

43ER 

RELRNR 

- 

24C0 

SOLS 

REVERR 

- 

8101 

0410 

ROES 

- 

SCO 

RV 

• 

8101 

04ER 

RVR 

- 

SOCO 

RVS 

- 

59C0 

04ft 

RVT 

• 

SOCO 

RVR 

- 

SOCO 

S 

- 

27CC 

SN1PT 

- 

uct 

SINOPT 

- 

78C0 

SNSSI 

- 

18C0 

SPREAR 

- 

23CO 

SSLS 

SPRER 

- 

3CC0 

T 

- 

TOCO 

80C0 

TPA 

- 

TOCO 

TNR 

- 

TOCO 

80E0 

TNTPAI 

• 

28C0 

TIPCPP 

• 

23CO 

SSLS 

TINE 

- 

SOCO 

TP  I N 

- 

67CO 

TOTALS 

- 

23CO 

SSLS 

TRS 

• 

11C0 

TRS2 

- 

11C0 

TASS 

- 

11C0 

TVS 

• 

23CO 

S5LR 

T WOP  I 

- 

14C0 

VPTTRN 

- 

2C0 

SOLS 

IS 

- 

28  CO 

1 

- 

TOCO 

SPIN 

- 

67C0 

■ PIT 

- 

43E8 

TORT 

- 

14C0 
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I H 1 C 

I 

P ABE 

6 

mo****  609 

1 

INCLUDE  DCONN 1 

DOP 

1 

2 

INCLUDE  DCONN 2 

DOP 

2 

3 

t 

DOP 

3 

4 

C 

DEAD  AND  PROCESS  INPUT , CARDS  AN 

TAPE# 

DOP 

4 

5 

c 

DOP 

5 

6 

9000 

CALL  IDENT 

DOP 

6 

7 

ir  (END)  65  TO  20000 

DOP 

7 

8 

c 

DOP 

a 

9 

c 

SET  UP  BAND  UNITS. 

DOP 

9 

10 

c 

DOP 

10 

11 

CALL  DCONN 

DOP 

11 

12 

c 

DOP 

12 

13 

c 

EXTRACT  SURE ACC  ANB/OR  BOTTOM  BATA  TROM 

INPUT  TAPE*  AND 

DOP 

13 

14 

c 

ARRAN6E  PROPERLY  f OR  DOP.  MR1TE 

ON  TENPORARY  FILE. 

DOP 

14 

15 

c 

DOP 

15 

16 

CALL  RBSORT 

DOP 

16 

17 

c 

DOP 

17 

18 

c 

PROCESS  TABLE  Of  REVERBERATION  TINES  AS 

NECESSARY  . 

DOP 

18 

19 

c 

DOP 

19 

20 

CALL  TCONP 

DOP 

20 

21 

c 

DOP 

21 

22 

NTNA1  * 0 

DOP 

22 

23 

c 

DOP 

23 

24 

c 

CONPUTE  REVERBERATION  TOR  EACH  TINE  IN 

TABLE. 

DOP 

24 

25 

c 

DOP 

25 

26 

10000 

NTN1N  » NTNA1  ♦ Kl 

•OP 

26 

27 

NTNAI  * NINO#  TNAX  ♦ INNT,  NTINE) 

DOP 

27 

28 

IF  t SMEAR!  MTMAI  • NTM1N 

BOP 

28 

29 

c 

DOP 

29 

30 

c 

7ERO  OUT  REVERBERATION  TABLE. 

DOP 

30 

31 

c 

DOP 

31 

32 

RO  10010  RT  • 1,  LNRB 

• OP 

32 

33 

REVERB(KT)  * 0. 

DOP 

33 

34 

10010 

CONTINUE 

• OP 

34 

35 

c 

• OP 

35 

36 

c 

CONPUTE  A PACE  FULL— THREE  TINES 

OR  ONE 

TINE  MITH  SPREADIN6. 

DOP 

36 

37 

c 

DOP 

37 

38 

BO  15000  KT  > NTM1R,  NTMAI 

DOP 

38 

39 

ITT  • IT  ♦ Kl  - NTNJN 

DOP 

39 

40 

c 

DOP 

40 

41 

c 

CONPUTE  BOUNDARY  AND  VQLUNE  REVERB • AND 

SPREAD  AS  REOU1RED. 

DOP 

41 

42 

c 

CONPUTE  TOTALS. 

DOP 

42 

43 

c 

DOP 

43 

44 

IF  (.NOT.  NOTAPE)  CALL  RBCONP 

DOP 

44 

45 

IE  (.NOT.  NOVOLN)  CALL  RVC0M9 

DOP 

45 

46 

IE  CSMRCAB)  CALL  RVSRRB 

DOP 

46 

47 

CALL  RTCONP 

DOP 

47 

48 

15000 

CONTINUE 

DOP 

48 

49 

t 

DOP 

49 

50 

c 

PRINT  REVERBERATION  BATA. 

DOP 

50 

51 

CALL  RVPRNT 

DOP 

51 

52 

c 

DOP 

52 

53 

c 

DOP 

53 

54 

IE  (NTMAI  .IT.  N TIME  > SO  TO  10000 

DOP 

54 

55 

60  TO  9000 

DOP 

55 

56 

c 

• OP 

56 

57 

c 

REWINB  TAPES  AS  REOUIREO.  ANB  EXIT. 

DOP 

57 

58 

c 

DOP 

58 

59 

20000 

IE  (NOTAPE)  SO  TO  20010 

DOP 

59 

60 

REMIND  AR1 

DOP 

60 

61 

REMIND  eM 

DOP 

61 

62 

20010 

IF  (PLOT)  REMIND  IPLT 

DOP 

62 

63 

3CC00 

STOP 

DOP 

63 

64 

END 

OOP 

64 
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inti 


JTN60L 

■ 

i * » » 

• ■ B 1 

9000 

6* 

55 

10000 

- 

26* 

54 

tooio 

- 

32 

34* 

15COO 

- 

38 

48* 

20000 

- 

7 

59* 

20010 

- 

59 

62* 

300C0 

- 

63* 

AM 

- 

60 

scour 

- 

11 

IM 

- 

61 

ICONkl 

• 

1 

DC0HN  2 

- 

2 

E H0 

• 

7 

1 01  NT 

- 

6 

IPLT 

- 

62 

(i 

- 

26 

39 

K T 

- 

32* 

33 

ATT 

- 

39* 

LNNT 

- 

27 

IPAB 

- 

32 

9IN0 

- 

27 

NOTAPE 

- 

44 

59 

NOVOLN 

- 

45 

NTINE 

- 

27 

54 

NTNAI 

- 

22* 

26 

NTNIN 

- 

26* 

28 

PLOT 

- 

62 

ABC  ONP 

- 

44 

ASSOAT 

- 

16 

Af SEAS 

- 

33* 

ATCOPP 

- 

47 

AVCOPP 

- 

45 

A VPANT 

- 

51 

RVSPAD 

- 

46 

SPREAD 

- 

26 

46 

STOP 

- 

63 

TCOPP 

- 

2u 

1PA1N 

- 

0 

umiNCis 


JO-  3* 


27*  28*  3*  54 

38  39 


PUI  7 


; 
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PA6C 

8 

1 

SUBROUTINE  IDENT 

IDENT 

1 

2 

c 

IDENT 

2 

3 

INCLUDE  D CORN  1 

IDENT 

3 

4 

COUPON  10(1? >, 1 , J ,K ,L,A 

IDENT 

4 

5 

LOGICAL.  OUT  US 

IDENT 

5 

£ 

DATA  INPf L6/0/,0UTfife/ .TRUE./ 

IDENT 

6 

7 

c 

IDENT 

7 

6 

c 

WRITE  HEAD  I N6  On  INPUT  DATA  PA4E. 

IDENT 

8 

9 

c 

IDENT 

9 

10 

WRITE  ( IPRT  f 5 01 ) HID 

IDENT 

10 

11 

501 

FORMAT  (8A6'A4,16M — INPUT  DATA  — /> 

IDENT 

11 

12 

c 

IDENT 

12 

13 

c 

READ  NEIT  VARIABLE  NAME  AND  RELATED 

DATA 

FROM  DATA 

CARD  (S). 

IDENT 

13 

14 

c 

IDENT 

14 

15 

o 

o 

CALL  INPUT(NNAMES,ID£V1NPFL6V0UTFL6> 

IDENT 

15 

It 

If  ( 60  > 60  TO  700 

IDENT 

16 

17 

IF  < END  > 60  10  900 

IDENT 

17 

18 

IF  F1NPF16  .67.  1>  END  « .TRUF. 

IDENT 

18 

19 

60  TO  600 

IDENT 

19 

20 

700 

60  * .FALSE. 

IDENT 

20 

21 

IF  (N0TAPE > 60  TO  2 000 

IDENT 

21 

22 

c 

IDENT 

22 

23 

( 

TAPE 

. 

IDENT 

23 

24 

c 

IDENT 

26 

25 

READ  «A»1)  <I0U>,  1*1,6) 

IDENT 

25 

26 

IF  <!6<6>  .NE.  0>  60  70  1000 

IDENT 

2t 

27 

c 

IDENT 

27 

28 

END  * .TRUE. 

IDENT 

28 

2» 

900 

WRITE  (IPRT ,9 01) 

IDENT 

29 

30 

901 

fOAPAT  (691, 18H***  END  OF  DUN  •*•> 

IDENT 

30 

31 

60  70  30000 

IDENT 

31 

32 

c 

IDENT 

32 

33 

c 

READ  SECOND  HEADER  RECORD.  PRESERVE 

ANT 

INPUT 

IDENT 

33 

34 

c 

TAPE  WHICH  HAS  NOT  BEEN  READ  IN  FROM 

CARDS. 

IDENT 

34 

35 

c 

IDENT 

33 

36 

1000 

READ  CARD  f 1 D( 1 ) r 1*6,  12) 

IDENT 

36 

37 

DO  1500  1 • 1.  12 

IDENT 

37 

38 

J * 1 

IDENT 

38 

39 

IF  (J  .6 T . K2)  J « J - Al 

IDENT 

39 

*0 

IF  (NAHCNT(i)  .ED.  0)  60  TO  1200 

IDENT 

40 

41 

IF  EJ  .Nf.  11)  60  70  1500 

IDENT 

41 

*2 

KS1  * KS1  • BE6RAD 

If  £AfT 

42 

*3 

60  70  1500 

IDENT 

43 

44 

1200 

IDATAfI)  * f D (I  ) 

IDENT 

44 

*5 

1500 

CONTINUE 

IDENT 

45 

46 

c 

IDENT 

46 

47 

c 

SCALE  INPUT  DATA  ABB  PRECOMPUTE  RELATED 

QUANTITIES 

IDENT 

47 

48 

c 

IDENT 

48 

49 

2000 

0ME6A  * ONE  SAD  • DE6RAD 

IDENT 

49 

50 

BELT  * BELT 2 / F2 

IDENT 

50 

51 

A SID  * KS1/DE6RAD 

IDENT 

51 

52 

BWI NT  * BVIDTH 

IDENT 

52 

53 

IF  (.NOT.  KTSBND)  BW1NT  * BW 1NT/F 1 E3 

IDENT 

53 

54 

L06MVI  * L06MV  - L064P1 

IDENT 

54 

55 

CSRSI  * COS(KSl) 

IDENT 

55 

56 

SNKS1  * SIV(RSl) 

IDENT 

56 

57 

DR  * DEL T2  • C0/F4 

IDENT 

57 

58 

COAT  * CO  • T DAT 

IDENT 

58 

5* 

FCOVS  » COAT/VS 

IDENT 

59 

60 

FC03  * C0**3  • INOfI 

IDENT 

60 

61 

IF  (FIDO  .It.  0.)  FIDO  • FI 

IDENT 

61 

62 

FISR  * F2R0**2 

IDENT 

62 

63 

FI2  * FIRO  • F2 

IDENT 

63 

64 

E IPS  - E IP (S/F10  • FLO610) 

IDENT 

64 

65 

IF  ITH7N6I  .CD.  0.)  7N7NA1  • F90 

IDENT 

65 

66 

IF  17N7N61  .67.  FIDO)  7N7RA8  • F180 

IDENT 

66 

67 

c 

IDENT 

67 

68 

NSPRNl  * A 1 

IDENT 

68 

69 

SPREAD  * SPREAD  .OR.  (NSSPRD  .N | . 0)  .OR. 

• OR. 

IDENT 

69 

70 

1 (NVSPRD.NE.O) 

IDENT 

70 
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71 

17  1.N01.  878  7*6  > 60  10  5010 

1*781  71 

72 

C 

1*781  72 

73 

C 

18773. 

1*781  73 

74 

c 

1*781  76 

75 

17  I (688***  .78.  01  .0*.  18887*6  .78.0)  .0*. 

(69871*  .78.  Oil 

1*767  75 

7* 

CALL  SPRCNP 

1*761  76 

77 

1*761  77 

7* 

c 

1*787  78 

7* 

5010 

887*8  • 887*81  - (1 

1RENT  79 

SO 

887*1  » 887*8  ♦ 887*81 

1*787  80 

SI 

17  (807*811  60  10  7010 

1*767  81 

S2 

c 

1*767  82 

85 

c 

SIT  UP  APPROPRIATE  PACE  NEARINCS. 

1*787  83 

84 

c 

1*787  86 

85 

R » 0 

1RENT  15 

86 

J • *5 

1*767  86 

87 

17  (.607.87*7**1  60  10  6020 
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SU0AOUTINE  ASSORT 
C 

INCLUDE  DCOAfc  1 
INCLUDE  DC0NN3 

COMMON  1»J,R«L,AvllvNAvlt*vNt**D»A  tD'TEST  tWRIL6 
LOGICAL  WML  4 
C 

USURP  * 0 
N8TA  * 0 

IE  (NOTAPE)  60  TO  16 
C 

TEST  * DBTTA  • SHlfT 
L • 0 
C 

C READ  DEFT M ID.  RECORD  AND  TEST  TOR  SURTACE. 

C 

1 MRFL6  * .RAISE. 

READ  CARD  I . D . N 

If  <1  .ED.  0>  60  TO  15 
If  CD  .NE.  C.)  60  TO  2 
If  CNOSURf > 60  TO  7 
NSURf  * N 
60  TO  5 
C 

C TEST  fOR  BOTTON  DEPTH, 

C 

2 If  CABSCD  - D8TTA)  ,6T.  TEST)  60  TO  7 
If  CNOBTTA)  60  TO  7 

NRTN  * N 
C 

C READ  DATA  FOR  NEXT  FATM.  SK  IF  PATHS  NOT  AT  CDESIRE»  SUkfACE 

C OR  BOTTON. 

C 

5 WRf L6  * .TRUE. 

7 DO  14  I * 1t  N 

READ  CARD  i9  CRBTHCR>fRBlCK>tRBHCK)vRDTCK)f  K * 1 f i ) 

If  C .NOT • WRFL6)  60  TO  14 
DO  9 K * I*  i 
RBHCK)  * RBHCK)  - S 
9 CONTINUE 
C 

C FUT  DATA  IN  ASCEND1N6  ORDER  Of  IF  NECESSARY. 

C 

If  CRBXC2)  . 6T.  RBI  CD)  60  TO  12 

NA  * i/2 

NB  * HLHKS 

DO  11  l * 1*  4 

NB  * NB  ♦ INKS 

NC  * NB  ♦ i 

DO  1C  N * 1«  NA 

NO  * NB  ♦ N 

A * RBICNO) 

rbiCno)  * rbicnc) 

RBICNO  * A 
NC  * NC  - 1 

10  CONTINUE 

11  CONTINUE 
C 

C SAVE  MINI  NUN  I AND  T FOR  EACH  PATH  t If  TINCNP  OPTION. 

C 

If  C .NOT.  TINCNP ) 60  TO  13 

12  L * L ♦ 1 

ININ  CL ) * RBI  Cl) 

THINCL)  * RBTC1) 

C 

C WRITE  DATA  ON  INTERMEDIATE  TAPE. 

C 

13  WRITE  IBRD  i vCRBICR)vRBTCK)fRBTHCK)fRBHCK)f  R = 1 , J) 

14  CONTINUE 
60  TO  1 


RBSORT  1 
RNSORT  2 
RBSORT  3 
RBSORT  4 
RBSORT  3 
RBSORT  6 
RBSORT  7 
RBSORT  8 
RBSORT  9 
RBSORT 10 
RBS0RT11 
RBS0RT12 
RBS0RT13 
RBS0RT14 
RBS0RT15 
RBS0RT16 
RBS0RT17 
RBS0RT18 
RBS0RT19 
RBSORT20 
RBS0RT21 
RBSORT22 
R6S0RT2? 
RBS0RT24 
RBSORT25 
RBSORT26 
ABSORT27 
RBSORT 28 
RBS0RT29 
RBSORT30 
RBS0RT31 
RBSORT  32 
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RPS0RTS8 
RNS0RT59 
RBS0RT60 
RBSCRT61 
RBS0RT62 
RBS0RT63 
RBS0RT64 
RBSORT65 
RPS0RT66 
RNS0RT67 
RBS0RT68 
RBS0RT69 
RBSORT70 
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18 

1 

SUBROUTINE  TCOHP 

TCOHP 

1 

2 

c 

TCONF 

2 

3 

INCLUDE  D CONN  1 

TCOHP 

3 

4 

INCLUDE  DCONN  3 

TCONF 

4 

5 

TIL 

TCOHP 

5 

6 
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TCOHP 

6 

7 

IF  ( • NOT  • IIIICIIP)  &0  10  1C 

TCOHP 

7 

8 

c 

TCOHP 

8 

9 
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T1NC  CONFUTATION  IS  DESIRED. 

NOTE  THAT  CONSTANT 

CHECKINS  IS 

TCONF 

9 

1C 

c 

DONE  TO  ENSURE  THAT  ALL  ADDED 

Tints  Lie  tcTnecn  i/z  mu  t 

TCONF 

1C 

11 

c 

AND  FIND  INTERVAL • NO  SUCH  CHECK  IS  NADE  ON  1NFUT 

tines  . 

TCOHP 

11 

12 

c 

ALSO  LEND TN  OE  THE  TINE  ARRAV 

0VER-E1LLIND 

•TCOHP 

12 

13 

c 

TCOHP 

13 

14 

IF  Inline  ,E«.  inn*)  tO  TO  10 

TCONF 

14 

15 

NTim  « MUM  ♦ K 1 

TCOHP 

15 

16 

IE  (NT IN£  .ED.  LNTIN)  DO  TO  9 

TCOHP 

16 

17 

TINE (NT INE)  * DELT 

TCOHP 

17 

18 

NTINE  * NT1NE  ♦ Kl 

TCOHP 

18 

19 

if  inline  .es.  LnTin)  to  to  9 

TCOHP 

19 

2C 

IE  (NOTAFE)  DO  TO  6 

TCOHP 

20 

21 

c 

TCOHP 

21 

22 
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CONPUTE  TINES  ASSOCIATED  N ITM 

OVER  EACH 

TCOHP 

22 

23 
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C0PBINAT10N  OE  TNO-NA T PATHS. 

TCOHP 

25 

24 

c 

TCONF 

24 

25 

NNAI  * 0 

TCONF 

25 

26 

DO  4 I BOUND  * Kl,  K2 

TCOHP 

26 

27 

IF  InBOOHBIlBOUHK)  .es.  0)  to  TO 

4 

TCOHP 

27 

28 

nnm  • nn«i  ♦ m 
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28 

29 

NNAI  * NNAI  ♦ NBOUNDdeOUND) 
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£9 

7C 

DO  3 I - NNIN,  NNAI 

TCONF 

30 

31 

READ  (B R 1 ) NTBL,(RBX<K>,RBT<K),RBTH(K),RBM(K),  K * 1 

• NT  DL  ) 

TCOHP 

SI 

32 

c 

TCONF 

32 

33 

DO  2 J * I,  NNAI 

TCONF 

S3 

34 

TineiNTine)  « (mini/)  « TAnunonini j>,*bi,i*t,i tnT.i >>  • fpts 

TCONF 

34 

15 

IE  (ITABLE  .EO.  0)  DO  TO  2 

TCONF 

35 

76 

IE  (TIFE(NTINE)  .BE.  FIND)  60  TO 

2 
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36 

37 

N * NTINE 
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37 

38 
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1 

TCOHP 

38 

39 
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TCONF 

3* 

4C 

IF  inTine  .tt.  LnTin)  to  to  9 

TCOHP 

40 

41 

TINE (NTINE)  « TIHE(F)  - DELT 

TCONF 

41 

42 

IE  (TINE (NT INE)  .ST.  DELT)  NTINE 

* NTINE  ♦ 11 

TCONF 

4? 

43 

IE  (NTINE  .ED.  LNTIN)  DO  TO  9 

TCOHP 

4.3 

44 

1 

T INE (NTINE)  * TINE(N>  ♦ DELT 

TCOHP 

44 

45 

if  (TineiNTine)  .it.  Pint)  NTine 

* NTINE  ♦ Kl 

TCONF 

45 

46 

if  InTine  .es.  LnTin)  so  to  9 

TCONF 

46 

47 

2 

CONTINUE 

TCOHP 

47 

48 

3 

CONTINUE 

TCOHP 

48 

49 

4 

CONTINUE 

TCONF 

49 

5C 

REWIND  0 R 1 

TCOHP 

sc 

51 

c 

TCOHP 

51 

52 

c 

ADD  PRESET  TABLE  OE  TINES. 

TCOHP 

52 

53 

c 

TCOHP 

53 

54 

6 

DO  7 1 = Kl,  LNTRS 

TCOHP 

54 

55 

IF  CTPSII)  . t E ■ PlhG)  to  TO  9 

TCOHP 

55 

56 

TINE (NTINE)  * TRS(I) 

TCOHP 

56 

57 

IE  (TINE (NTINE)  .6T . BELT)  NTINE 

* NTINE  ♦ Kl 

TCOHP 

57 

58 

IE  (NTINE  .ED.  LNTIN)  DO  TO  9 

TCOHP 

58 

59 

7 

CONTINUE 

TCOHP 

59 

6 C 

c 

TCOHP 

60 

61 

c 

CONTINUE  WITH  EVERT  1/2  SECOND 

TO  FIND.  ADD  F1N6 

INTERVAL. 

TCOHP 

61 

62 

c 

TCONF 

62 

63 

A 

TINE (NTINE)  * T INE (NT IFE-1 ) ♦ EPT5 

TCOHP 

63 

64 

IE  (TINE(NTINE)  .DE.  FIND)  60  TO 

9 

TCOHP 

64 

65 

NTINE  * NTINE  ♦ Kl 

TCOHP 

65 

66 

IE  (NTINE  .NE.  LNTIN)  DO  TO  6 

TCOHP 

66 

67 

c 

TCOHP 

6* 

68 

9 

TINE (NTINE)  * P1N6 

TCOHP 

60 

69 

c 

TCOHP 

69 

70 

c 

SORT  TINE  ARRAT  INTO  ASCENBIND 

ORDER,  EL1FINAT1N6 

DUFLI CATES. 

TCOHP 

70 
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SUBROUTINE  T C ORP 


10  CALL  S0RT(TIRE,NT1NF> 

INITIATE  PALE  COUNTER  ANR  COMPUTE  TOTAL  NUMCER  OF  PACES  Of 
OUTPUT  FOR  THIS  FILE  OF  INPUT  OATA. 

PACE  « 0 

NNAI  • ENTIRE  - ID/K3  ♦ Kl 
IF  (SPREAD  NNAI  • NT  IRE 

NPACE  * (CENNBNR  - NPS TR T > / IPEV R V > /I 4C  ♦ ID  • NNAI 
IF  (TOTALS)  NPACE  * (NTIRf  - RD/140  ♦ Kl 


TCONP  71 
TCONP  7? 
TCONP  73 
TCONP  74 
TCONP  75 
TCONP  76 
TCONP  77 
TCONP  7 f 
TCONP  79 
TCONP  60 
TCONP  81 
TCONP  62 
TCONP  83 
TCONP  84 


OD  52258 


1 N D E 


SUBROUTINE  TCO«P 
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SfPbOl 

* * * * 

1:11 

* x r x 

X X * 

REFERENCES 

1 

_ 

38 

44  • 

2 

- 

t 3 

35 

36 

67* 

3 

- 

3U 

68* 

6 

- 

26 

27 

69* 

6 

- 

20 

56* 

7 

- 

56 

59* 

6 

- 

63* 

66 

9 

- 

16 

19 

6C 

63 

66  55 

ir 

- 

7 

16 

72* 

B»1 

- 

3 1 R D 

50 

DCOPN1 

- 

3 

DCOPN3 

- 

6 

BELT 

- 

17 

38 

61 

62 

66  57 

FPT5 

- 

36 

63 

1 

- 

SCO 

30* 

33 

56* 

55  56 

1 BOUND 

- 

5C0 

26  = 

27 

29 

1 P E ¥ R T 

- 

8 0 

I TABLE 

- 

35 

J 

- 

5C0 

33  = 

36 

K 

- 

5C0 

31RD 

K 1 

- 

15 

18 

26 

26 

39  62 

76 

8C 

81 

56  66  68* 


65  56  57  65 


K? 

- 

26 

1 3 

- 

78 

K6G 

- 

80 

81 

LPT  IN 

- 

16 

16 

19 

6C 

43 

46 

56 

66 

LPT  R S 

- 

56 

- 

SCO 

37  = 

61 

66 

NNBND 

- 

8C 

NPOUND 

- 

27 

29 

NBA* 

- 

SCO 

25  = 

28 

29  = 

30 

33 

76= 

79= 

80 

NPIN 

- 

SCO 

28* 

30 

NOTAPE 

- 

20 

NPA6E 

- 

8 0* 

81  = 

NPSTRT 

- 

fu 

NTGL 

- 

SCO 

31RD 

36 

NT1PE 

- 

16 

15  = 

1e 

17 

18  = 

19 

34 

36 

37 

36 

39* 

60 

61 

62  = 

43 

44 

4 5* 

46 

56 

57* 

56 

63 

66 

65  = 

66 

68 

7 2A6 

76 

79 

81 

PAGE 

- 

77* 

PING 

- 

36 

65 

55 

64 

68 

RPN 

- 

31RD 

RPT 

- 

31RD 

36 

RBTM 

- 

31RD 

RBI 

- 

31RD 

36 

RETURN 

- 

83 

SORT 

- 

72 

SPREAD 

- 

79 

TAfcLRP 

- 

36 

TCOPP 

- 

1 

TIPCPP 

- 

7 

TINE 

- 

17* 

36* 

36 

38 

61  = 

42 

4 4= 

45 

56= 

57 

63* 

66 

68* 

72A6 

TP  IN 

- 

36 

TOTALS 

- 

81 

TRS 

- 

55 

56 

IP1N 

- 

36 
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1 c E 
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1 

SUBROUTINE  RBCOMP 

RBCOHP  1 

2 

C 

RBCOMP  2 

i 

INCLUDE  DCONM 

RBCOHP  3 

4 

INCLUDE  DCO*N  2 

RBCOHP  4 

5 

INCLUDE  DC0MN5 

RBCOHP  5 

6 

COMMON  DD(2),XD(?),TO(?),THAD(2),THBD(2),RD(2) 

RBCOHP  6 

7 

CO. "ON  CSTH»M21  ,CSTHBD<2>iO.TO<2>  .CSO.TD ( 2 >, SkO. TO ( 2) 

RBCOHP  7 

6 

COMMON  XTMU,0MTTMU,C0MTMU,S0MTHU,CTMTMU(2> 

RBCOHP  8 

9 

COMMON  RT(2),THAT(2),THBT<2) 

RBCOHP  9 

10 

RBCOHPIO 

11 

COMMON  CSALF A ,C  S A L F P , C SB  E T A , C SB  E T B , C SG  AM  A , C S G A MB 

RBCOHP 1 1 

12 

COMMON  CALFAP  fC A L F6P ,C GAMAP ,CG AM8P ,C0$ A ,CO$B 

ROCOMPI 2 

13 

COMMON  PH  I, PM MONT  ,DELR,TMMAX ,TVM1N ,T TM A X t A , DOP SKP 

RBC0HP13 

14 

C0..0.  «1(2>,>2<  2>.«3<2>,IU<2>,V1<2)>V2(2>,T3(2I,T4(2) 

15 

COMMON  I ,11 , 1 J,J A.JB.JBOUND.JC, JD,JDD,JDELT,JDMA|,JDM1N, 

JE 

RBCOMP 1 5 

16 

COMMON  JNAI .JUMP , K 8 AND ,K T 1 M , L ,L E AP ,L L ,MLMU < 2 ) , MM ,MN , Ml ,N A , NB , N E X T 

RBC0..16 

17 

EQUIVALENCE  < *L * U ( 1 > , ML > , ( ML MU < 2 > ,MU > , < DOP , 1 DOP) , < MNMI ,M N ) 

RBC0HP17 

18 

EQUIVALENCE  ( QC 0 S ,0  A ABC ) , CQS I N ,Q» ABC ) , (QC OSM VQ AB > 

RBC0HP18 

19 

DIMENSION  ID0P<1 ),MNMlf?> 

RBC0MP19 

70 

LOGICAL  TTNA  1 ,D  OPSKP 

R6COHP20 

21 

C 

RBC0MP2 1 

22 

c 

BOUNDART  LOOP  — 

RBCOHP22 

23 

c 

ONCE  THROUGH  EACH  FOR  SURFACE  AND  BOTTOM 

RBCOMP23 

24 

c 

RBC0MP24 

25 

CO  9C1C  J60UHC  * 1,  2 

R8C OMP25 

26 

J M A X = N80UND (JBOUND) 

RBCOMP 26 

27 

If  ( JM A X .|Q.  J)  GO  TO  9C10 

RBCOHP27 

28 

c 

RBCOHP28 

29 

c 

PATH  LOOP-- 

RBCOMP 29 

30 

c 

ONCE  THROUGH  FOR  EACH  SURFACE  REFLECTED  PATH  AND  EACH 

BOTTOM 

rbco..3o 

31 

c 

REFLECTED  PATH  ON  THE  INTERMEDIATE  TAPE.  (PATH  A> 

RBCOMP  3 1 

32 

c 

RMCOHP32 

33 

DO  801C  JA  * 1,  JMAI 

R6COMP33 

34 

c 

RBCOMP 34 

35 

c 

POSITION  TAPE  AT  JA'TH  RECORD. 

R8COMP35 

36 

c 

MBC0HP36 

37 

IF  (JA  .£Q.  1 > GO  TO  1G20 

RBCOMP37 

38 

CO  1011  1 * JA'J.JX 

RBCOMP38 

39 

BACKSPACE  BR1 

RBC0MP39 

40 

1010 

CONTINUE 

RBCOMP 40 

41 

c 

RBC0HP41 

42 

1 02C 

REAO  (BR1)  NA  , ( * B 1 A < 1 ) , R BT A ( 1 > , R6 T HA ( I ) , R BH A ( 1 > , I * 1, 

NA) 

RBCOMP42 

43 

BACKSPACE  BR1 

R6C0MP43 

44 

c 

RBC0MP44 

45 

c 

PATH  COMBINATION  LOOP  — 

R8COMP45 

46 

c 

ONCE  THR0U6H  FOR  EACH  COMBINATION  OF  PATHS  TO  SURFACE 

RBC0MP46 

47 

c 

OR  TO  BOTTOM.  (PATH  B> 

R BC0MP47 

4f 

£ 

RBC0MP48 

49 

DO  7080  JB  * JA,  JM  A X 

R6C0MP49 

50 

READ  (BAD  NB,(RBXB(1>,RBTB(I),RBTHB(I),RBHP(I),  1 « 1, 

NB  ) 

RBCOHP50 

51 

c 

R6C0HP51 

52 

c 

DEFINE  THE  AREA  COMMON  TO  THE  TWO  PATHS  AND  ALSO  TO  THE 

RBC0HP52 

53 

c 

CURRENT  TIME  INTERVAL,  IF  ANY  SUCH  AREA  EXISTS. 

RBCOHP  53 

54 

c 

RBC0HP54 

55 

TVMIN  * AMAX1  (TIME ( X T ) -D E L T , RBT  A ( 1 ) ,RB  TB  ( 1 > ) 

RBC0..55 

56 

TB.Itl  « CM*  1 < T I »•  E <r))*Oll  T,liBTjm»«>,*eTB(N8>> 

RBC0HP56 

57 

IF  ( T MM  A 1 .LE.  TNMIN)  GO  TO  7080 

RBCOHP57 

58 

c 

RBCO.*58 

59 

CO  1030  JC  * 1.  h* 

•BC0HP59 

60 

X(JC)  * RBIA(JC) 

RBCOHP60 

61 

1G3U 

CONTINUE 

RBC0MP61 

62 

JC  * NA 

RBCOMP62 

63 

c 

PMCOMP63 

64 

DO  1040  1 « 1 , NA 

RBC0MP64 

65 

JC  « JC  ♦ «1 

66 

I ( JC  > * RP«B 1 I) 

RBC0MP66 

67 

1 04  o 

CONTINUE 

RBC0HP67 

6> 

CALL  SORT (1 , JC) 

RBCOMP68 

69 

c 

RBC0MP69 

70 

c 

EVALUATE  TABLES  FOR  OVERLAPPED  ARIA. 

RBCOMP  70 
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SUBAOUT1PE  RBCOPP 


71 

c 

RBCONR71 

72 

4D  • 0 

RBCOPP72 

73 

74 

II  * Kl 
14  * K1 

-JIBCOPP73 

^IMKcoppt* 

75 

60  1050  1 * 1 . it 

ABCOBP75 

76 

RBCONR76 

77 

IF  (ITABLE  .ER.  0)  60  TO  1050 

RBCONR77 

70 

II  * ITABLE 

RBC0PP78 

79 

T ( JO* 1 ) - C T I JO*  1 1 « TABLRPIKD.ABRB.ABTB.IJ.BB))  * 

f NTS 

to 

IF  (ITABLE  .6*.  07  *0  TO  1050 

RBC0PP80 

81 

14  * ITABLE 

RBCOPPtl 

82 

J»  * 4B  ♦ Kl 

83 

X (4B)  * xcn 

8* 

THB(J0>  * AIPTRP(ABTL6>  • 0E6RA0 

85 

MJP)  * AINTRR(RBMB) 

RBCOPP8S 

86 

THA (JO)  * TABLKPiACD.RBIA.ABTHA.Il.RA)  • 0E6RA0 

RBCONR86 

87 

R ( JO  > • ETPCFL0610  • (AINTRP (RBMA  > * R(J0))/F20> 

RBCOPP87 

88 

COSTHACJB)  * COS  CTMACJ0)) 

RBCONR88 

89 

COSTNB(JO)  • COS (THB(JO)) 

*0 

ONTCJD)  * ANOBC0NESA  * TCJB),  TVORl) 

91 

IF  (ABs(0KT(J0>>  .6T.  PI)  0NT(J0>  m OBT(JO)  - SI6P(TVOPI,  OPT < JO > >R0C0BP91 

92 

COSOBT(JO)  * COS (ORT(JO)) 

RBCONR92 

93 

SIPORT(JO)  • SI6M(SRRT(F 1 -CO SORT (J0)**2)»0RT(J0)) 

RBCOPP93 

94 

o 

o 

CONTINUE 

95 

C 

RBCONR95 

96 

c 

transnit-receive  LOOR — 

RBC0NR96 

97 

c 

JC  * 1»  TRANSNIT  RATH  Av  RECEIVE  RATH  B 

98 

c 

JC  * 2*  TRANSPUT  RATH  B»  RECEIVE  RATH  A 

RBC0NR98 

99 

c 

RBC0NR99 

100 

00  7070  JC  * 1.2 

RBCOPIOO 

101 

c 

RBCON101 

102 

c 

TORE  “ ATT  LOOR— 

RBC0N102 

103 

c 

JE  » 1,  FORBARO  HER1SPMERE.  OR  RECEIVE  FRERUENCV 

BREA TER 

RBCON103 

10* 

c 

THAN  TRANSR1T  TR  EBVENCT • 

RBCON104 

105 

c 

JE  • 2.  AFTER  MEPISPPERE.  OR  RECEIVE  FAER.  LESS  TRAN  SPIT 

FREB.RBCOPlOS 

106 

c 

RBCOPIOO 

107 

00  7050  JE  * 1,  2 

RBCOP107 

108 

LL  * (JC  - Kl)  • LPKS2 

RBCON108 

109 

NN  * LHKS2  - LL 

RBCOP109 

110 

M • K3  - JE  * X2 

RBCOP110 

111 

NCOS  * f A 

RBCOP111 

112 

•COSH  « M 

RBCOP112 

113 

00  2020  1 * 1,  JO 

RBCOP113 

11* 

NL  • LL  ♦ I 

RBCOP11* 

115 

HU  * NN  ♦ 1 

R8COP115 

110 

IF  (0PE6A  .ER.  FO)  60  TO  2010 

RBC0N116 

117 

RA  * COSTHA(PU)  • COSBpT (I)  « COSTMA(PL) 

RBCOP117 

118 

QB  * COSTHA INU>  * SINBNT <I> 

RBCOP118 

119 

RC  * RB  • COSTHA (Pi l/FCOVS 

RBC0P119 

120 

BAB  * BA«*2  ♦ «B««2 

RBCOP120 

121 

BCOS  * (BB  • BC  * FA  • RA  • SORT (BAB  - *C**2))/RAB 

RBCON121 

122 

RSIP  • (BB  * BCOS  ' BO/BA 

ABCOP122 

123 

•COSH  * (BCOS  • COSOPT ( I > * BS1H  • SINOPT(l)) 

RBCOP123 

12* 

2010 

OOP(l)  - F2R0  * (COKE  * VS  * COSTBA(PD)  * BCOSP) 

RBCOP12* 

125 

/ (COKT  - VS  • COSTPA(PL)  • BCOS) 

RBCOP12S 

126 

2020 

CONTINUE 

RBCOP126 

127 

NN  « Kl 

RBC0P127 

128 

BORSKR  • .FALSE. 

RBCON128 

129 

SORTNU  * 0. 

RBC0P129 

130 

c 

RBCOP130 

131 

c 

LOOP  TO  PROCESS  COBSECUTIVE  STRIN6S  OF  BATA  HPICP 

ARE 

ROCOP131 

132 

c 

POPOTORIC  IP  'BOP' 

RBCOP132 

133 

€ 

RBCOP133 

13* 

3010 

JBELT  • 1 SI8P (1 . IBOP(RP)  • 1B0P(PP*1>) 

RBCOP13* 

135 

PL  * R1  ♦ (JBELT  ♦ ID/I2 

RBCOP135 

136 

NU  * K3  - NL 

RBCOP136 

137 

U * JO  - *1 

RBCOP137 

138 

00  3020  1 • PP.  IJ 

139 

IF  (JBELT  .RE.  lSISPd,  lOOP(I)  - I00P(I*1>)>  SO  TO 

3030 

ROCOP139 

1*0 

3020 

CONTINUE 

RBCOPIOO 
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in>c 

I 

SUBROUTINE  RBCOHP 

PARE  23 

141 

1 * 60 

RRC0R141 

142 

c 

FBCOM142 

143 

3030 

III  « X 

FBCOB143 

144 

60H1N  • hnni  col) 

FBCOM144 

14} 

J0NAI  • UNO! (MUX 

FBCOB145 

146 

c 

FBCOM 146 

147 

0041)  ■ 006  C J 004 1 ) 

FBCOB147 

141 

IF  (JE  .1*.  K 21  66(1)  > FIFO 

FBCOM148 

149 

600  ■ 60NAI 

FBCOB1 4* 

150 

6SS14N  3050  10  NEXT 

OBCOH1SO 

151 

• 0 3040  KBANF  > MSPFH1,  MAN* 

OOCOH1S1 

152 

IF  <te<1>  .41.  BANt((BAN*«1>>  40  70  (3070,  30«0>,  JE 

FBCOB152 

153 

o 

* 

o 

continue 

FBC0B153 

154 

40  70  7030 

FBCOM154 

155 

C 

XBCOM155 

156 

c 

100P  FO*  AIL  8ESIFE*  BANFF  AN*  FOB  ALL  FAIA  A1  C0NS7AN7  60PPL EBFBC0P 156 

157 

c 

WITHIN  4 0 AN0  * If  BAN0W10TH  tS  40EATER  THAN  OATA  SPAC1N4. 

FBC0B1S7 

151 

c 

IN  THE  fOLLONINS  T40LE S • THE  SUBSCRIPTS  1 ANO  2 REfER  TO 

FBCOB1S8 

159 

c 

OATA  POINTS  Of  THE  LOWEST  ANO  N1SHEST  fREOUENCT  — SUOSCR IPTS 

FBCOB13* 

160 

c 

Hi  ANO  NU  Off EO  TO  POINTS  OE  THE  LOWEST  ANO  HISNEST  TIRES. 

O0COH16O 

161 

c 

08COH161 

162 

3050 

ASSI4N  4020  10  NEI7 

R0COR162 

163 

c 

FBCOB 163 

164 

3060 

0012)  » 00(1) 

OBCOR164 

165 

lot?)  - 1011) 

FBCOB 165 

166 

T0(?)  * TOtl) 

OOCOH166 

167 

THAO ( 2)  » TNAOtl) 

FBCOB 167 

16* 

7BBFC  2)  - 1MBB1 1 ) 

FBCOB168 

169 

0042)  • 0011) 

O0COR169 

170 

CSTHAOt  2)  • CSTNA0C1) 

O0COP17O 

171 

CS1HBFI  2>  • CI1HBK1) 

O0COR171 

172 

ONTO!?)  • ON  TOtl) 

FBCOB172 

173 

CS0M1D(2>  • CSOMIFOI 

FBC0B173 

174 

SN0HT0C2)  • SRONT0C1) 

O0COR174 

173 

c 

XBCOB175 

176 

600  » 600  » 60ELT 

R0COH176 

177 

IF  (F0P(JF*»  .LI.  0ANB((*F***1>>  *0  70  3080 

RPC OH  177 

17* 

0011)  * 0OPC6OO) 

O0COR178 

179 

3070 

lOt 1 ) • 11600) 

FBCOB17* 

ISO 

TOtl)  • T (600 ) 

FBC0B180 

101 

TRAOff)  * r*4<6PP> 

O0COH181 

1*2 

IBB 6 1 1 1 • IBB  OFF) 

FBCOB  162 

1*3 

ROt 1 ) « 04600) 

RBCON183 

184 

CSTHAOt 1)  ■ C OSTHA 4600) 

OPC  ON  184 

185 

CSTN00I1)  - C OSTHO 4600) 

BBC  ON  1 85 

186 

ONTO 4 1 ) • ONT  1600 ) 

O0COH186 

187 

CS0N1*(1>  • COSOB!  OFF) 

0RC0N187 

it* 

R0COH188 

189 

80  TO  4010 

O0COH189 

1*0 

c 

O0COH19O 

191 

3080 

600  * 600  - 60ELT 

FBCOBIfl 

1^2 

FBCOB1V? 

1*3 

30*0 

FBCOP193 

194 

FOPSKP  - 17 ABLE  .E*.  0 

O0COH194 

1*5 

If  400PSKP)  40  TO  NEIT 

FBCOB  1*5 

196 

TOID-  AINTOPIT) 

O0COH196 

197 

THAO  1 1 ) « AINTRP ITHA) 

FBC0B1«7 

ivt 

THRO  ID  » AINTRP4THB) 

FBCOB 1 *6 

1** 

FB1 1 > > A1N7FPEF  > 

OPC  ON  199 

200 

CSTHAOt 1 ) « A INTOP 4 COS THA) 

FBCOB200 

201 

CS1MBF1 1 > « A INI FPL CO* 1MB > 

R0CON2O1 

202 

ONTO  CD  * AIN  TOP  40NT) 

R0COR2O? 

203 

CSOB1FC1)  * A INI FP (COS OBI) 

RPCOH203 

204 

SNOHTOI 1 ) - A1NTOPISINONT) 

R0COH2O4 

205 

c 

FBCOM205 
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- 
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- 
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3 39 

340 

PI 

- 

91 

• A 

- 

10C0 

117* 
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18EB 
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- 

10CO 

18EB 
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294 
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BABC 

- 

10CO 

294* 

295 
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- 

10C0 
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BB 

- 

10CO 
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289* 
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BBABC 
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311 
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293 
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56 
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- 
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93 
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- 
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123 
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- 
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- 
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race  >1 


SUBROUTINE  RVCOMP 

INCLUDE  DCONN1 
INCLUDE  D CONN  ? 

INCLUDE  DC0NN4 

CONNCN  I « J ( K v NN  , ONT  ,C SONT F .SNOMTF t FR #f R St ,F IS* .F I » C 
CONN  ON  $NPHI(?),CSPN](2)tX(2)t!(2>tZ(?)llt»lfltIIV 
CONNON  T 1 *T 2 •a1»A2«V0L1«T3 


COMPUTE  ONE-BAT  RANGE  TO  MIDPOINT  OF  TRANSMITTED  PULSE* 
PRE-COMPUTE  RELATED  VALUES  FOR  VOLUME  AND  FOR  RAN6E  LOSSES* 

R * T INf  CRT > • CO  • FPT5 

FRY  * EXPMLOtNVl  - ALPNC  • TIME(KT)>/F10  • FL0G10)/R**4 
IF  fONEGA  •€  0 • 0.)  FRV  * FRV  * 

1 (CR  ♦ Dl)««3  - AMAXlCR  - DR*  fOMO>  • PI23 

IF  ( VPTTRN)  GO  TO  9C00 

PATTERN  LOSSES  MUST  RE  COMPUTED* 

OMT  * AMODCOMEGA  • T1MEUTI*  TVOP1) 

IF  (A8S<0MT>  *ST.  PI>  OMT  « OMT  - SIGNdMOPl,  OMT) 

CSOPTF  * COSt  OMT  > 

— «W»W- •-*  S16JKSD*XIF1-C$0MTFa*2),©MT)  • FZRO 

C SOM T F * CSOMTF  • FZ2  

FCSGAMtNSPRHl ) * 0. 

IF  ( OMEGA  • ED • 0.)  FC SGAM(NSPRNl)  * FI 
T 1 ^ aMAI 1(T IMe(KT)  - DELT*  FO)  « FPT5 
12  * ((TIME(RT)  ♦ DELT)  a FPT5  - TD/F10 

COMPUTE  AVERAGE  PATTERN  LOSSES  FOR  E ACM  RAND* 


DO  8C00  1 * NSPRMl,  ND AND 
FGAM(I)  * 0. 

FR  * CBANDt  I > ♦ BANDtl  + D)  a FPT5 
FRSB  * FR  **  2 

FISB  *■  FISB  ♦ FR  SB  ♦ FR  a CSOMTF 
FX  * SORT (FXSO) 

A 1 * FCOVS  • (FR  - FZRO) 

CSA  * A 1 / FI 

IF  (ABS(CSA)  *€T*  f1>  CSA  * SI 6*<F1,  CSA) 

SNA  * SORTCF 1 - CSA**2) 

NN  * MIN0CK2  * 1 F II  (FNBAND  • SNA)  ♦ Kl , 360) 
PN  * NN 
DD  « TMOP1/PN 
D * 0* 

SNPH 1(1)  a SNORTF/FR 

SNPN1C2)  a -SNPH 1(1)  • FR/FZRO 

CSPHId)  > SORT  ( F 1 - SNPHltl)**?) 

CSPNK2)  * SORTCF1  - SNPHI(2)**2) 

IF  (ABStFZSB  - FRSB)  *LE*  FXSO)  GO  TO  2000 
CSPHId)  a S1GN(CSPNJ(D*  CSA) 

CSPNM2)  a -SIGN(CSPMKI)*  CSA) 

2000  DO  4000  J * 1 VNN 

D » D ♦ DD 

B a SNA  t S1N(D) 

G * SNA  * C0S4D) 

C 

DO  3000  K * 1 f 2 

1(1)  * CSA  • CSPHId)  - D A SNPH  1 (K ) 

Yd)  * B • CSPHId)  ♦ CSA  A SNPH  Id) 

Id)  a Xd)  • CSX  SI  A G A SNKS1 

Id)  * G A CSKSI  - X (X  > A SDKS  I 

3000  CONTINUE 

F€AM(I)  * F G AMO  ) A OIL  UNIT  v!  (1 ) *¥  (1 ) ,Z  €1  > ) 
1 a OXL(RECV»X(?)»T(2)»Z(2)) 

4000  CONTINUE 

FGAM(X)  a FGAMCD/PN  a E IPS 
IF  (OMEGA  .ED.  0.)  GO  TO  6 000 
T3  * ?1 


RVCOMP  1 
R VCOMP  2 
RVCOMP  3 
RVCOMP  4 
RVCOMP  5 
RVCOMP  6 
RVCOMP  7 
RVCOMP  8 
RVCOMP  9 
RVCOMP10 
RVCOMP11 
RVCOMP12 
RVCOMP13 
RVCOMP14 
RVCOMP IS 
RVCOMP 16 
R VC0MP17 
RVCOMP18 
RVC0MP19 
RVCOMP 20 
RVCOMP21 
RVCOMP22 
RVCOMP23 
RVCOMP24 
RVCOPP2S 
RVCOMP26 
RVC0MP27 
RVCOMP28 
RVCOMP29 
RVCOMP30 
R VCOMP31 
R VCOMP32 
RVCOMP33 
RVCOMP34 
RVCOMP3S 
RVCOMP36 

RVCOMP 37 

RVCOMP38 

RVCOMP39 

RVCOMP40 

RVCOMP41 

RVCOMP4? 

RVCOMP43 

RVCOMP44 

RVCOMP45 

RVC0MP46 

RVCOMP47 

RVCOMP48 

RVCOMPSO 

RVCOMPS1 

RVCOMP52 

RVCOMPS3 

RVC0MPS4 

RVCOMPS5 

RVCOMPS6 

RVCOMPS7 

RVCOMP38 

RVCOMP39 

RVCOMP60 

RVCOMP61 

RVCOMP62 

RVCOMP63 

RVC0MP64 

RVCOMP65 

RVCOMP66 

RVCOMP67 

RVCOMP68 

RVC0MP69 

RVCOMP70 
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SUBROUT  ICE  RVCONP 
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A2  » ri2  • mi 

RVCONP71 
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A1  « Al/SORT < A2) 

RVCONP72 
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RVCONP73 
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VOLl  * 0. 

RVCONP74 
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c 

RVCONP75 

7* 
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RVCONP76 
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71 
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T3>)*FCi>RVCONP78 

7* 

5000 

CONTINUE 

RVC0NP79 

80 
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RVCONP80 

81 

COTO  70C0 

RVCONP81 

82 
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FCSCANClO)  * f C CVS  * <BAND<I«1>  - El  RO)  / <B  AN  D(  1 ♦ 1 ) 

♦ F 2 R 0)  R VCONP82 

83 

7000 

F6AN  ( I > * FGANCI)  • ( F CSC AN  < 1 ) - FCSCAMI«1>> 

R VCONP8I 

84 

RVV(KTT*I>  * FRV  * ECAN(I) 

RVCONP84 

85 

8000 

CONTINUE 

RVCONP85 

86 

c 

RVCONP86 

87 

c 

FOR  STRAICHT-RUNN1N6  CASE*  PATTERN  LOSSES  FOR  ANT 

BAND 

BILL  RVC0NP87 

88 

c 

BE  THE  SANE  AT  ALL  TINES. 

RVCONP88 

89 

c 

RVC0NP89 

90 

IF  (ONEGA  • E 0 • 0.>  VPTTRN  * .TRUE. 

RVC0NP90 

91 

CO  TO  30000 

RVC0NP91 

92 

9000 

DO  1CC00  I = NSPRHl*  NBAND 

R VCONP92 

93 

R VV (KTT  * I ) * FRV  * F6AN ( 1 ) 

RVCONP93 

94 

1GC00 

CONTINUE 

RVCONP94 

■95  " 

3'COGQ 

RETURN 

RVCONP95 

96 

END 
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23 
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62 

52 
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- 

62 

63 
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- 
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24 
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7C0 
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B 

- 
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46* 
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56 

57 

BC0NN 1 

- 

3 
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- 

4 

BC0HN4 

- 

5 

BB 

- 

7CO 

45* 

55 

BELT 

- 

28 

29 

B* 

- 

16 

EIP 

- 

14 

EIPS 

- 

68 

ro 

- 

16 

28 

78 

' 

M 

- 

24 

27 

41 

62 

69 

1 

MC 

- 

14 

29 

f€03 

- 

80 

MOTS 

- 

39 

82 

fCSBAN 

- 

26* 

27* 

80* 

82* 

83 

MAR 

- 

34* 

65* 

68* 

83* 

86 

FL0B10 

- 

14 

- 

43 

MT5 

- 

13 

28 

29 

35 

M 

- 

6C0 

35* 

36 

37 

39 

MSA 

- 

6C0 

36* 

37 

51 

73 

MV 

- 

7CO 

14* 

15* 

86 

93 

M 

- 

6C0 

38* 

40 

67 

MSB 

- 

6CO 

37* 

38 

51 

M2 

- 

25 

71 

f 1*0 

- 

24 

39 

48 

82 

MS* 

- 

37 

51 

73 

C 

- 

6CO 

57* 

62 

63 

1 

- 

6C0 

33* 

34 

35 

65 

i 

92* 

93 

If  11 

- 

43 

f 

2 

- 

6C0 

54* 

76* 

K 

- 

6CO 

59* 

60 

61 

62 

* 1 

- 

43 

*2 
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43 

AT 
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13 

14 

21 

28 

29 

ATT 

- 

84 

93 

«* 

LOBNV 1 

- 

14 

•• 

NINO 

- 

43 

* •> 

NBANB 

- 

33 

92 

$ 

NN 

• 

6C0 

43* 

44 

56 

i # 

NSPRN 1 

• 

26 

27 

33 

92 
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• 

15 

21 

27 

69 

78 
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- 

6C0 
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26 

60  61 


>0  12  S3  84 


APPENDIX  G 


ism 


OD  52258 


1 N » C * 


SUBROUTINE  RVCONP 


OIL 

- 

65 

66 

PI 

- 

22 

PI23 

- 

16 

PN 

- 

6C0 

44* 

«5 

68 

R 

- 

7CO 

13* 

14 

16 

RECV 

- 

66 

RETURN 

- 

95 

RVCONP 

- 

1 

RVV 

- 

84* 

93* 

SION 

- 

22 

24 

41 

52 

S3 

SIN 

- 

56 

SNA 

- 

6C0 

42* 

43 

56 

57 

SNR  S 1 

- 

62 

63 

SNONTF 

- 

6C0 

24* 

47 

SNPHI 

- 

7CO 

47* 

48* 

49 

50 

SORT 

- 

24 

38 

42 

49 

50 

T1 

- 

8CO 

28* 

29 

70 

T2 

- 

8CO 

29* 

77 

80 

T3 

- 

SCO 

70* 

77* 

78 

TINE 

- 

13 

14 

21 

2« 

29 

T WOP  1 

- 

21 

22 

45 

VOL  1 

- 

SCO 

74* 

78* 

80 

• 

17 

90* 

X 

- 

7CO 

60* 

62* 

63 

65 

XN1T 

- 

65 
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- 

7C0 

61* 

65 

66 

Z 

- 

7CO 

63* 

65 

66 
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INCLUIE  iCONMI 
1NCIUIE  RCORN  2 
COMMON  IMAI'1,1 


10  13  i ■ Kl • K3 

•0  11  I • NSPRNl • MfAMR 

11  * I-OIMH 

00  10  I ■ I1»  MSM1 
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RV<?,UvJ>  • »V(2,IM)  ♦ RV<1,IvJ>  • SME»CKKVJ> 
11  * II  ♦ 1 

If  CI1  .Cl.  UIIANO)  CO  TO  11 

10  CONTINUE 

11  CONTINUE 
13  COMTlNUf 

RETURN 
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RVSPRR  1 
RVSPRR  2 
OVSMO  3 
IVSPRt  4 
IVSfRD  5 
IVSMO  6 
mom  7 
RVSNO  8 

msom  o 
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RVSPRR11 
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RVSPRR13 

RVSPR014 

RVSPRR15 

RVSPRR16 

RVSPR017 
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14 
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11 

©COMBI 

- 

3 

0COHN? 

- 

4 

1 
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SCO 

8* 

11 

- 

SCO 
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- 
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• 

SCO 

7* 

K 

- 

5C0 

10* 

K 1 

• 

7 

10 

K3 

- 

7 
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- 

SCO 

11* 

- 

14 

NBANB 

- 

8 

NSPR1 

- 

10 

NSPRN 

- 

9 
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- 

8 

11 

RETURN 
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18 
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- 
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1 
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12 
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NOE 

1 

RA6E  37 

i 

SUBROUTINE  RTCONR 

RTCONR  1 

2 

€ 

RTCONR  2 

3 

INCLU6E  6C0NN1 

RTCONR  5 

4 

INCLURE  RC0NN2 

RTCONR  4 

S 

RTCONR  5 
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C 

RTCONR  6 

T 

IRAK  » KTT 

RTCONR  7 

• 

IE  (.NOT.  SRRER9I  60  TO  2 

RTCONR  6 

t 
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RTCONR  9 
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TO  2 

RTCONRIO 

11 

1NAX  * 95 

RTC0NR11 

12 
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RTCONRl? 

1) 

IE  (TV6»  T • TV6E(TINE(9T» 

RTC0NR15 

1* 

• 011-1,  NNBNR 

RTC0NR16 

1$ 

TR  • T 

RTC0NR15 

1» 

IE  (EILTERI  TR  » T * RRESd) 

RTCONR 16 

17 
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11 
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RTCONRIO 

1* 
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Rvad.ji 

RTC0NR27 

26 

RVV(1,4) 

RTCONR26 

29 

c 

ft  TC0NR29 

50 

4 R VT ( I ,4 ) « RVSCI  ,4)  ♦ RV6<1,4> 
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IE  (J  .LI.  NN6NR)  60  TO  5 
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5 CONTINUE 
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1 

SUBROUT  INF  RVPRNT 

RVPRNT  1 

2 

c 

RVPRNT  2 

X 

INCLUDE  D CON  N 1 

RVPRNT  3 

4 

INCLUDE  D CONN  ? 

RVPRNT  4 

5 

COMMON  ItJvKvL,II 

RVPRNT  5 

6 

0 IRENS  ION  HRVBU) 

RVPRNT  6 

7 

DATA  HT1NE/6HTINE  " / , ( H R V P < I > , I * 1,6» 

RVPRNT  7 

8 

1 /6H  SUR,  6MFACE  • 6MB0TT0P,  6M  VOL,  6hUPE  , 6H  TOTAL/ 

RVPRNT  8 

9 

DATA  ITOTAL/1/ 

RVPRNT  9 

1C 
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RVPRNT 1 0 

11 

IF  (PLOT)  WRITE  < 1 PL  T ) T 1 ME (NTN I M ) ,HNBN D , (UNDO UT < K ) ,BNDOUT (R ♦ 1 ) , 

RVPRNT 1 1 

12 

1 RVS  (2vK),RVe(?(K),RVV(2 vK),RVT(2fK) , K * 1,  HNBND) 

RVPRNT 1 2 

13 

IE  (NOPRNT)  00  TO  3C000 

RVPRNT13 

14 

IF  (TOTALS)  60  TO  2000 
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15 
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17 
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RVPRNT 17 
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1 VS,C0,FIROfDC,S,KSIDfOHE6AD,P3NW,DElT2,DeTTH 

RVPRNT18 

19 

If  (lOTACSI  to  TO  2920 

RVPRNT 19 

2C 

( 

RVPRNT20 

21 

IF  (SPREAD)  WO  TO  5 

RVPRNT  21 

22 

WRITE  ( I PRT  , 1 1 ) MPAND,  (MT1ME  VT1**E  (K>  , K * NTR1N,  NTPAX) 

RVPRNT  2 2 

23 

WO  TO  6 

24 

c 

23 

5 

UNITE  (IPRT, 11)  HEANDt<HTIHE.TI«E<NTNIN>,  a * 1,  I> 

R VPRNT2  5 

26 

c 

27 

6 

WRITE  (IPRT, 12)  MOUTPT,  (HR V6  , 1 * 1,  I) 

RVPRNT27 

2£ 

If  (TOTALS)  tO  TO  2100 

R VPRNT28 

29 

DO  1 A - A 1 * At 

RVPRNT29 

3C 

UNITE  (IPPT.20) 

RVPRNT30 

31 

c 

RVPRNT31 

32 

DO  2 L * Kit  R5 

RVPRNT32 

33 

WRITE  (IPRT,  2t>  BNDOUTW),  ONDOUT (J • 1 > . 

RVPRNT33 

34 

1 (RVSOl.J),  DuPEII.i).  UVVdliJ),  RVTUI.J),  II  * 1,  I) 

RVPRNT34 

35 

J * J ♦ IPEVRY 

RVPRNT  35 

36 

IF  ( J • 6 T . AN  PSD ) WO  TO  4 

RVPRNT36 

37 

c 

RVPRNT  37 

3f 

c 

CONTINUE 

RVPRNT38 
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3 

CONTINUE 
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60  TO  1 

RVPRNT40 

41 
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RVPRNT41 
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4 
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RVPRNT42 

43 
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R VPRNT43 

44 

WO  TO  3COOO 

RVPRNT44 

45 

c 

RVPRNT45 

4e 
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If  (.NOT.  SPREAD)  1 * R1 

RVPRNT46 

47 
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RVPRNT47 

4t 

II  * N T*  I N 

RVPRNT48 

49 

2 Cl  u 

ITOTAL  * 1T0TAL  - Hi 

RVPRNT49 

5 C 

If  (ITOTAL  .WE.  C>  tO  TO  2100 

R VPRNT  50 

51 

ITOTAL  * K43 

RVPRNTS1 

52 

WO  TO  1 

RVPRNT52 

53 

c 

RVPRNT53 

54 

2020 

WRITE  < IPRT  ,13)  MTOT1,MNbNd«HTOT2,BW1DTH,HTOT3 

RVPRNT  54 

55 

WO  TO  6 

RVPRNTS5 

56 

c 

RVPRNT56 

57 
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If  (NOD (ITOTAL.AS ) .ED.  0)  UNITE  (1PNT.20) 

RVPRNT57 

58 

UNITE  (IPRT, 21)  TINE(II),  (RVS(A,LN8ND1),NVB(A,LNBND1), 

RVPRNT  58 

59 

1 RVV(R,LHBND1> ,RVT(K  ,LMBND 1) , K * J,  I) 

R VPRNT  59 

60 

1 * 1 ♦ A1 

RVPRNT60 

61 

J * J ♦ A 1 

RVPRNT61 

62 
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RVPRNT62 

63 

If  (11  .IE.  NTNAX)  tO  TO  2010 

R VPRNT63 

64 

c 

RVPRNT64 
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RETURN 
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1 • 1,  31/1(0  040401470  / 


4474  <07071(11, 
4474  <07072(11, 
4474  <07073(11, 


CLOCK 
CLOCK 
CLOCK 
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AIM  - 1 SIA 

•»1  - 1 5BA 
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ICORNA  - 4 

if AAAI  • 12IA 

>1  - SIA 

>10  - 9IA 
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fU3  - 13IA 
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>3  - 13IA 

>4  - 9IA 

>90  - 14IA 

>40*10  • 1 OIA 
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HEAIS  - 17IA 
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HT0T1  - 2 SIA 

MT0T2  - 26IA 

NT0T3  - 27IA 

I - 7IA  OIA  9IA  1 7|A  2SIA  260*  27IA  31IA 

HINT  - 101 A 

HIT  - 1 SIA 

mviT  - 2siA 

mt  - i sia 

1MT  - 1 SIA 

ITAINN  - 2SIA 

J - 31IA 
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(1  - 6IA 
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A 3 - 6IA 

(40  - AIA 
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(A  - AIA 

(0  - AIA 
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TIS2  - ABA 
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1 

SUBROUTINE  SORT  f AM  AT  «L  ENCTH  > 

SORT 

1 

2 

DIMENSION  ARRAV(I) 

SORT 

2 

3 

C 

SORT 

3 

4 

c 

SORT  IN  ASCENDING  ORDER  (SIMPLE  REPLACEMENT 

SORT) 

SORT 

4 

5 

c 

SORT 

5 

6 

If  (LENGTH  .LE.1I  GO  TO  4 

SORT 

6 

7 

DO  2 4 * 2, LENGTH 

SORT 

7 
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9 

TEMP  * ARRAV(i) 
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SORT 
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16 

17 
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SORT 
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SORT 
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SORT 
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CONTINUE 

SORT 
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SORT 
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33 
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SORT 

33 
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FUNCTION  TABLCPC  ARB , 1NDEP VDE PNO  , 1 TARNN , I TABPI > 

COPPfN  /CTBLKPS  PEL INP , DEI BE P , f AC TOA , I T ABLE 
omr.NSiON  iNOCP(l)vt£rNo(1) 

BEAL  INBEP 

OELIND  * IND  E P(  2 ) - JNDEP  (1) 

BO  3 ITABLE  * 1TABNN,  ITABM1 
FACTOR  * ABB  - I NBE P < 1 T ABLE > 

IF  (FACTOR  • OELIND)  5,  6,  3 
3 CONTINUE 


12 


U 

17 

18 

19 

20 
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4 ITABLE  * 0 
BO  TO  7 

5 IF  (ITABLE  . E 9.  1>  BO  TO  4 

OELIND  * INDEPd TABLET  - I NO E P ( I T ABLE -1 > 
FACTOR  * FACTOR/DEL1NO 

6 TABLKP  * A1 NT BP( OEPND ) 

7 RETURN 
END 


TABLKP  1 
TABLKP  2 
TABLKP  3 
TABLKP  4 
TABLKP  5 
TABLKP  6 
TABLKP  7 
TABLKP  8 
TABLKP  9 
TABLKPIO 
TABLKP11 
TABLKP12 
TABLKP13 
TABLKP14 
TABLKP15 
TABLRPU 
TABLKP17 
TABLBP1S 
TABLKP19 
TABLKP20 
T ABLKP21 
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rui 


SYMBOL  Ilf CRINCIS 


1 

• 

* 

10* 

— . 

2 

• 

i 

12* 

AIMTRM 

- 

i 

7* 

11* 

CTBLRM 

- 

3(t 

BCLBEM 

- 

3(0 

10* 

11 

BEilMB 

3(0 

BlfMB 

• 

uc 

*01 

7 

10 

11 

EACTCR 

• 

3(0 

* 

11 

STABLE 

- 

3(0 

7 

10 

11 

RETURN 

- 

12 

Mil  53 


SUBROUT INE  INPUT (NAPTAO. RATA. 1NPFL6.0UTFL6! 

INPUT 

1 

: 

INPUT 

2 

COHHON  /C  *NCOO/  lNAP.IlNA6CtlPACE<60!,lBLANK.NUNBER(21) 

INPUT 

3 

EGUl VALENCE  ( IGUOT E , NUHB C R ( 1 7) ! v f 1 CONN* .NUHBE R (20 > ! 

INPUT 

4 

EQUIVALENCE  ( IE GUAL .HUNGER 4 21 !! 

INPUT 

5 

COPPON  /CNACMN/  NVORR #NCNAR .PAXCOL 

INPUT 

6 

: 

INPUT 

7 

LOGICAL  BCRFL6*RECPT.E IPNT.C XB0TM.0UTFL6. LVALUE. L TRUE 

INPUT 

8 

ROUBLE  PRECISION  VALUE*R10 

INPUT 

9 

RIPEN  SION  NABTARCll.GATAd!  .AVALUEd! 

INPUT 

10 

E GUI  VALENCE  ( VALUE • AVALUE . 1 VALUE .LVALUE VNANE  > * (1 . 1 IP AGE! . 

INPUT 

11 

1 LNTBUE.LT RUE  > 

INPUT 

12 

RATA  ITTPE/0/ .4/ 3 / . I TRUE / 5NTRUE •/.R10/1 »R1/ .FRNAT1/6H(80A1 >/ . 

INPUT 

13 

1 LTRUE/.TRUE./ 

INPUT 

14 

INPUT 

15 

1TTPE  * 1*  INTEGER  (NO  REGIHAL  POINT! 

INPUT 

16 

ITYPE  ■ 2.  REAL  CNITH  OR  UITNOUT  EXPONENT! 

INPUT 

17 

1TVPE  * 3.  ROUBLE  PRECISION  CO-TYPE  EXPONENT! 

INPUT 

18 

1TVPE  * 4.  CONPLE1  LTVO  REAL  NUNOERS  IN  PARENTHESES! 

INPUT 

19 

1 TYPE  * S.  LOGICAL  (.TRUE.  OR  .FALSE*  ONLY! 

INPUT 

20 

I TYPE  * 6.  ALPHABETIC  (ENCLOSER  IN  APOSTROPHES.  THUS  *ABC"! 

INPUT 

21 

: 

INPUT 

22 

GO  TO  (1260.  1600,  1030.  2000!#  4 

INPUT 

23 

1C30 

REAR  FRPAT1,  IHAGE 

INPUT 

24 

IF  COUTFLG!  PRINT  2001.  IHAGE 

INPUT 

25 

1*0 

INPUT 

26 

GO  TO  1600 

INPUT 

27 

1100 

1 * I ♦ 1 

INPUT 

28 

IF  (ITYPE  .EG.  6!  GO  TO  1160 

INPUT 

29 

IF  ( IHAGE  ( I ) .EG.  1BLANR ! 60  TO  19C0 

INPUT 

30 

IF  (IHA6E(1)  .NE.  1CONHA!  60  TO  (1130#  1110!*  4 

INPUT 

31 

lies 

60  TO  ( 1 2 9C  . HOC),  4 

INPUT 

32 

1110 

RO  1120  X * 1.  19 

INPUT 

33 

IF  (1HA6E(I)  .NE.  NUHBERCK!!  GO  TO  1120 

INPUT 

34 

IF  (K  .LT.  18!  60  TO  1300 

INPUT 

35 

IF  (RECPT!  60  TO  1300 

INPUT 

36 

1120 

CONTINUE 

INPUT 

37 

IF  <L  ♦ NAPE  .NE.  0!  60  TO  1995 

INPUT 

38 

4 * 1 

INPUT 

39 

IF  (RECPT)  ITYPE  * 5 

INPUT 

40 

1130 

IF  <IPA6E(1!  .EG.  I EQUAL ) 60  TO  1200 

INPUT 

41 

1140 

IF  (L  .LT.  NVORR)  60  TO  1150 

INPUT 

42 

IF  (ITYPE  .NE.  6)  60  TO  (1900.  1210).  4 

INPUT 

43 

IF  (BCRFL6!  1*1-1 

INPUT 

44 

60  TO  (1500*  1210)#  4 

INPUT 

45 

1150 

FLR(L.NCHAR.NANE)  * FL R <0 .NCHAR . IPA6E ( 1 ! ) 

INPUT 

46 

L * L ♦ NCHAR 

INPUT 

47 

60  TO  (1900.  1140).  4 

INPUT 

48 

1160 

IF  (BCRFL6!  60  TO  1170 

INPUT 

49 

IF  ( IHAGE (1 ! .EG.  ICOPHA!  60  TO  1200 

INPUT 

5C 

60  TO  1900 

INPUT 

51 

1170 

IF  ( IHA6E (I ) .NE.  1 QUOTE ) 60  TO  1140 

INPUT 

52 

NGT  * 1 - NQT 

INPUT 

53 

IF  (NGT  .EG#  0)  60  TO  1140 

INPUT 

54 

IF  (IPA6E (I ♦ 1 ) .NE.  IGUOTE!  BCRFL6  * .FALSE. 

INPUT 

55 

60  TO  1900 

INPUT 

56 

1200 

4*2 

INPUT 

57 

K * I 

INPUT 

58 

1 * 81 

INPUT 

59 

60  TO  1140 

INPUT 

60 

1210 

I « X 

INPUT 

61 

IF  (ITYPE  .GE.  5!  GO  TO  1500 

INPUT 

62 

NNAHES  * NANTABd! 

INPUT 

63 

1NPFLG  « 0 

INPUT 

64 

NAPS A V * NANE 

INPUT 

65 

RO  1220  X « 1,  NNAHES 

INPUT 

66 

IF  (NAPE  .NE.  NAHTAP (K+1 ! ) GO  TO  1220 

INPUT 

67 

IF  (ITYPE  .EG.  0!  60  TO  1260 

INPUT 

68 

4*1 

INPUT 

69 

60  TO  1999 

INPUT 

7C 
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INftEI 

SUBROUTINE  INPUT CNAHTAB 

PAST  54 

71 

1220 

CONTINUE 

INPUT  71 

72 

E0  TO  1996 

INPUT  72 

73 

1260 

NL  * NNANE5  ♦ I « 1 

INPUT  73 

74 

NH  * NN ARES  ♦ NL 

INPUT  74 

75 

N C - NNARES  ♦ NR 

INPUT  75 

76 

INPUT  76 

77 

17  EIPAtEEI)  ,E0  • I E9UAL > *0  TO  1600 

INPUT  77 

78 

J • 2 

INPUT  78 

79 

NANTABCNO  * NTRUE 

INPUT  79 

eo 

to  TO  1998 

INPUT  80 

81 

1300 

If  <ft  .41.  1C ) 90  TO  1320 

INPUT  81 

82 

If  f .NOT • E ■ PNT > SO  TO  1310 

INPUT  82 

83 

I EIP  « I|*P  • 10  ♦ S 14m (K  - 1, 

1ES1SN) 

INPUT  S3 

84 

SO  TO  1900 

INPUT  84 

85 

131u 

If  (6ECPT  > L * L ♦ 1 

INPUT  85 

86 

VALUE  • VALUE  * HO  ♦ DftLE  ( fLOA  T (SISNCK  - 1,  IVS1SN))) 

INPUT  86 

87 

SO  TO  1900 

INPUT  67 

88 

1320 

I * K - 10 

INPUT  88 

89 

to  TO  C1330,  1330,  1350,  1360, 

1380,  1900,  1390,  1325,  1335),* 

INPUT  89 

90 

1325 

ITTPE  « 3 

INPUT  90 

91 

1330 

If  (.NOT.  OECPT)  60  TO  1340 

INPUT  91 

92 

If  <K  .EO.  2)  1ESI6N  * - 1 

INPUT  92 

93 

1335 

EIPNT  * .TRUE. 

INPUT  93 

94 

SO  TO  1900 

INPUT  94 

95 

1340 

If  Cl  .Eft.  2 > 1VSI6N  * - 1 

INPUT  95 

96 

to  TO  1900 

INPUT  96 

97 

1350 

ftECPT  * .TRUE. 

INPUT  97 

98 

ITTPE  * 2 

INPUT  98 

99 

SO  TO  1900 

INPUT  99 

100 

1360 

IREPT  * VALUE 

1NPUT1CC 

101 

to  TO  1620 

1NPUT101 

102 

1380 

ITTPE  * 4 

INPUT  102 

103 

SO  TO  1900 

1NPUT103 

1C* 

1390 

ITTPE  « 6 

INPUT  1 04 

105 

BCftfLS  * .TRUE. 

INPUT  1 35 

106 

* * i 

1NPUT1G6 

107 

MOT  * 0 

INPUT 1C7 

108 

60  10  1900 

INPUT  1 0F 

109 

1*00 

IF  EINPFLt  .ME.  0)  fO  TO  16C0 

INPUT10« 

110 

L * 1 EIP  - L 

INPUT  110 

111 

I E IP  « ABS<L> 

INPUT111 

112 

00  1*30  I « 1 , I EIP 

INPUT  112 

113 

IF  EL)  1*10,  1**0,  1*20 

INPUT113 

114 

1410 

VALUE  * VALUE/ftIC 

1NPUT1 14 

115 

to  TO  1*30 

1 NPUT 115 

116 

1420 

VALUE  * VALUE  • HO 

1NPUT116 

117 

1430 

CONTINUE 

INPUT117 

118 

1440 

If  CITTPE  .NE.  4)  SO  TO  1500 

INPUT  118 

119 

IF  ECIBOTH)  fO  TO  1500 

INPUT  119 

120 

CIVAL  * VALUE 

INPUT  1 20 

121 

ClftOTN  * .TRUE. 

INPUT  121 

122 

SO  TO  1610 

INPUT122 

123 

u» 

o 

c 

IF  EINPFLf  .ME.  0)  tO  TO  1600 

1 NPUT 1 2 3 

124 

ftO  1590  1*1,  IREPT 

I NPUT 124 

125 

1505 

INftEI  * NANTABCNC)  ♦ 1 

INPUT12S 

126 

If  C INftEI  .ST.  NART  Aft f NR>)  SO 

TO  1997 

INPUT  126 

127 

NANTABCNO  « INftEI 

INPUT127 

120 

INftEI  * INftEI  ♦ NARTABCNL) 

INPUT  128 

129 

to  TO  <1520,  1530,  15*0,  15*0, 

1560,  1580*  15  7 0,  ITTPE 

INPUT129 

130 

1520 

I VALUE  * VALUE 

1 NPUT 1 30 

131 

1525 

0ATAE1N0EI)  « AVALUE  <1 > 

INPUT  131 

132 

to  TO  1590 

1NPUT132 

133 

1530 

0A7AC INftEI)  * VALUE 

INPUT  133 

13* 

to  TO  1590 

1 NPUT 134 

135 

15*0 

INftEI  * INftEI  ♦ NANTABCNO  - 1 

INPUT135 

136 

If  CITTPE  .Eft.  A)  SO  TO  1550 

1NPUT136 

137 

• ATA  C IftftE IM  > * AVALUE  C2> 

INPUT137 

138 

to  TO  1525 

I NPUT 1 38 

139 

1550 

OATAEIN9EI)  * CIVAL 

1NPUT139 

1*0 

0ATAEIN6EIA1 ) • VALUE 

1 NPUT 140 
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SUBROUTINE  INPUT (NAMTA8 ,DAT A , 1NPFL6 ,OUT F L6) 


MCE  55 


41  (0  TO  1590 

42  1560  LVALUE  « •FALSE. 

43  If  (NAPE  .£6.  I 7 RUE  ) LW.-i.Uc  * .TSUE. 

44  60  TO  1*25 

45  1570  N > INDEX  - 1 EXP 

46  DATA (INDEX)  * DATA(N) 

47  IF  in  .EO.  L > 60  TO  1585 

48  60  TO  1 5C5 

49  1580  OATA(INDEX)  * AV  ALU E ( 1 1 

5C  I EXP  * I EXP  ♦ 1 

51  IF  C6CDFL6)  60  TO  1610 

52  1TTPE  - 7 

53  1585  L • INDEX 

54  1590  CONTINUE 

55  1600  CXB01M  * .FALSE* 

56  INEPT  * 1 

57  J » 2 

58  X * 0 

59  ITTPE  * 1 

60  1610  L * 0 

61  IEIP  * 0 

62  1ESI6N  * 0 

63  1VSI6N  * 0 

64  DECPT  * .FALSE. 

65  EXPNT  * .FALSE. 

66  1620  VALUE  * 0. 

67  1900  IF  (I  .LT.  NAXCOL)  60  TO  1100 

68  IF  (K  .NE.  0)  60  TO  1105 

69  60  TO  1030 

70  1995  IF  (1NPFL6  .NE.  0>  60  TO  1900 

71  INPFL6  * 1 

72  1996  INPFL6  * 1 ♦ INPFL6 

73  1997  1NPFL6  • 1 ♦ 1NPFL6 

74  J « 4 

75  1998  ITTPE  * 0 

76  1999  NETUNN 

77  2000  IF  (1NPFL6  .Et.  0)  60  TO  16C0 

78  IF  (.NOT.  0UTTL6)  PRINT  2001,  1MA6E 

79  2001  FORMAT  (20l,80Al> 

8C  00  2010  J « 1,  90 

81  IMA6E  (J  > * I BLANK 

82  2010  CONTINUE 

83  1NA6E  (I INA6E  > ■ NUNBER ( 14  ) 

84  PRINT  2001,  INA6E 

85  J * 2 

86  60  TO  (2020,  2030,  2040),  INPFL6 

87  2020  PRINT  2021,  NANSAV 

88  2021  FORMAT  (20X,  17NT00  MUCH  DATA  IN  , A6'/> 

89  60  TO  1999 

90  2030  PRINT  2031,  NANSAV 

91  2031  FORMAT  <20X,  A6,  17N  NOT  IN  NAME  LIST/) 

92  60  TO  1999 

93  2040  PRINT  2C41 

94  2041  FORMAT  (20K,  33MSTNTAR  ERROR  OR  1LLE6AL  CHARACTER/) 

95  60  TO  1999 

96  END 


1NPUT141 
JNPUT142 
INPUT  143 
INPUT  1 44 
INPUT  145 
INPUT  146 
1NPUT147 
INPUTU6 
1NPUT149 
1NPUT150 
1NPUT151 
INPUT  152 
1NPUT15I 
INPUT154 
INPUT  155 
INPUT  156 
INPUT157 
INPUT  158 
INPUT  1 59 
1NPUT16G 
INPUT  161 
1NPUT162 
INPUT163 
INPUT  164 
INPUT  165 
1NPUT166 
1NPUT167 
INPUT  168 
INPUT169 
INPUT  170 
1NPUT171 
INPUT172 
INPUT  1 73 
1NPUT174 
INPUT  1 75 
1NPUT176 
INPUT  177 
INPUT  178 
INPUT179 
INPUT180 
INPUT  181 
INPUT  1 82 
INPUT  183 
INPUT  184 
INPUT  165 
INPUT  186 
INPUT187 
1 NPUT 188 
INPUT189 
JNPUT190 
INPUT  191 
INPUT  192 
INPUT  193 
1 NPUT 1 94 
INPUT  195 
I NPUT 196 
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SUeROUTim  INPUT<NA«TA8,DATA,lNPri4,0UmG> 


SfMOt  .»•««...  r .«  f ..  t KEMIEKCiS 


103C 

. 

23 

24* 

169 

1110 

- 

26* 

167 

11C5 

• 

32* 

Uf 

1110 

- 

31 

*3* 

lire 

• 

33 

34 

37* 

me 

• 

31 

41* 

me 

- 

42* 

48 

52 

54 

60 

IISt- 

- 

42 

46* 

11*? 

- 

29 

49* 

me 

- 

49 

52* 

me 

- 

32 

41 

50 

57* 

1210 

- 

43 

45 

61* 

me 

- 

66 

67 

71* 

1260 

- 

23 

68 

73* 

1500 

- 

35 

36 

•1* 

131C 

- 

(2 

85* 

m: 

- 

(1 

88* 

ms 

- 

89 

9C* 

1T3C 

- 

89 

91* 

ms 

- 

89 

93* 

1340 

- 

91 

95* 

use 

- 

89 

9 7* 

1 !fcC 

- 

(9 

ICO* 

1J6C 

- 

89 

1 02* 

13*0 

- 

89 

104* 

HOC 

• 

?2 

1C9* 

1410 

- 

111 

114* 

1420 

- 

113 

116* 

143C 

- 

112 

115 

117* 

1440 

- 

113 

118* 

1300 

- 

45 

62 

118 

119 

123* 

1503 

- 

125* 

148 

1520 

• 

129 

130* 

1323 

- 

131* 

138 

144 

153C 

- 

129 

133* 

134: 

- 

129 

135* 

15  30 

• 

136 

139* 

156C 

- 

129 

142* 

1520 

- 

129 

145* 

1520 

- 

129 

149* 

15(3 

• 

147 

153* 

15*0 

- 

124 

132 

134 

141 

154* 

1(00 

• 

23 

27 

77 

109 

123 

133* 

177 

1610 

- 

122 

131 

160* 

1620 

- 

101 

166* 

1900 

- 

30 

43 

48 

51 

36 

86 

(7 

99 

103 

108 

167* 

170 

1995 

• 

38 

170* 

1996 

- 

72 

172* 

1997 

• 

126 

173* 

1991 

- 

83 

175* 

1999 

- 

73 

176* 

189 

192  . 

195 

2000 

- 

23 

177* 

2001 

• 

259* 

178M 

179* 

1849* 

2010 

- 

180 

182* 

2020 

- 

186 

187* 

2021 

- 

1879* 

188* 

2030 

• 

186 

190* 

2031 

- 

1*39* 

191* 

2040 

• 

186 

193* 

204 1 

- 

1939* 

194* 

493 

• 

111 

4V4LII2 

- 

1091 

1119 

131 

137 

169 

acm* 

- 

816 

44 

49 

55* 

103* 

151 

CIHC06 

- 

3(L 

C94CMN 

• 

6(1 

CIIOTH 

• 

8L6 

119 

121* 

155* 

ca«*l 

- 

120* 

139 
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SUMOUTJNt  lNPUTO.««T»B,D»T*  ,lNMl4,0ljTri4> 


rut  57 


DIO 

- 

voe 

130A 

86 

114 

116 

DATA 

- 

1 AG 

1061 

in- 

133- 

137* 

139* 

1*0- 

146* 

149- 

• PIE 

- 

66 

OECPT 

- 

bl  G 

36 

40 

85 

91 

97* 

164* 

E*PNT 

- 

8LG 

62 

93* 

165* 

no 

- 

46* 

hoai 

- 

86 

mmti 

- 

1366 

266* 

1 

- 

TIER 

26* 

28* 

30 

31 

34 

41 

44* 

50 

52 

55 

58 

59- 

61* 

77 

167 

161 ANK 

- 

SCO 

30 

181 

X COMMA 

- 

4ER 

31 

50 

1 EQUAL 

- 

5EO 

41 

77 

1ES1GN 

- 

63 

92* 

162* 

ll»f 

- 

63- 

lie 

in* 

112 

145 

150- 

161* 

1 IMAGE 

- 

SCO 

11f« 

183 

IMAGE 

- 

SCO 

2400 

25PR 

30 

31 

3* 

41 

5 0 

52 

55 

77 

178PR 

181* 

183* 

184PR 

I NAM 

- 

SCO 

1NOE* 

- 

125- 

126 

127 

128* 

131 

133 

135- 

137 

139 

140 

U5 

tit 

149 

153 

INPELG 

- 

TAG 

64* 

109 

123 

170 

171- 

172- 

173- 

177 

186 

INPUT 

- 

1 

IfiUOTE 

- 

4ER 

52 

55 

10EPT 

- 

1C3- 

124 

1 56* 

1 TRUE 

- 

1 30  A 

113 

1 TYPE 

- 

1 SO  A 

29 

40* 

43 

62 

68 

9 j* 

98* 

102- 

104- 

118 

129 

1*6 

152- 

159* 

175- 

1 VALUE 

- 

TIER 

13C- 

1 VS  ION 

- 

86 

95* 

163* 

4 

- 

1S0A 

23 

31 

32 

39* 

43 

45 

48 

57* 

69* 

7i- 

106* 

157- 

17*. 

180* 

181 

185* 

K 

- 

33- 

34 

35 

58* 

61 

66* 

67 

72 

81 

83 

66 

88* 

89 

92 

95 

112- 

124* 

158* 

160 

L 

- 

36 

62 

46 

4 7* 

85- 

110- 

in 

111 

147 

153- 

160* 

LTPUE 

- 

8LC 

1286 

140  A 

LVALUE 

- 

316 

her 

142- 

1*3- 

M 

- 

145- 

146 

147 

MAICOL 

- 

6C0 

167 

NAME 

- 

TIER 

30 

46 

65 

67 

143 

NAMSAV 

- 

t 5* 

107PO 

190P* 

- 

TAG 

1001 

63 

67 

76- 

79* 

125 

126 

127* 

128 

135 

NC 

- 

75* 

76 

79 

125 

127 

135 

NCNAA 

- 

6C0 

46 

47 

NL 

- 

73* 

74 

128 

NM 

- 

76* 

75 

1?6 

NNAMES 

- 

63* 

66 

73 

74 

75 

NOT 

- 

53* 

54 

107* 

NTAUE 

- 

1 1E  0 

79 

NUMBER 

- 

3C0 

4ER 

5ER 

34 

183 

NMORO 

- 

6C0 

42 

OUTfLG 

- 

1 AG 

6LG 

25 

178 

Of TURN 

- 

176 

SIGN 

- 

33 

86 

VALUE 

- 

90S 

TIER 

86* 

100 

114- 

116- 

120 

130 

133 

140 

166- 
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•LOCK  BATA  COB  INK U1  tUOOOUTINC 

t BLOCK  BATA  BLOCK  1 

2 COMMON  /CINCOB/  IRAN .1 |M ABE  t|N ABE  CB01 .1BLANR .KUMOF A ( 21 2 BLOCK  2 

3 BATA  IBLANK/ IN  / BLOCK  3 

A BATA  (NUNBEBCl>tl>1.21>/1N0.1N1,1N2,1N3.1N4tlN3.1N«.1N7.1Nt,1N*.  BLOCK  4 

3 1 1N«,lN-,1N..1N»t1NC,1N>,1N*,1NB,lNC,1N,,lN»/  BLOCK  3 

« C BLOCK  A 

7 C MOTE  THAT  TNC  EOLLOM1NC  BATA  It  NACNINE  BEMENBENT.  BLOCK  7 

t C NMOBB  • NO.OE  BITS  IM  NACNINE  MOBB.  BLOCK  ( 

B C NCNAB  • NO.OE  BITS  IN  ALBNA "BUNCO 1C  CNABACTCB . BLOCK  * 
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